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Rapid Prototyping, RP in short, is a layered manufacturing process to make pro- 
totypes fast from CAD data. RP is emerging as a key prototyping technology with its 
ability to produce even complicated parts within hours. RP is a term which embraces 
a range of new technologies for producing accurate parts directly from CAD models, 
in a few hours, with little need for human intervention. This means that designers 
have the choice of making physical models of the conceptual design more frequently, 
allowing them to assess the form, fit and functionality of the design. This also allows 
them to discuss manufacturing issues with an easy-to-interpret prototype which can 
be held, viewed, studied, tested and compared. As a result, a product is improved 
and the product development cycles and costs substantially reduced. 

RP systems take the solid model of a part as an input and make a physical model 
or prototype without using tools or fixtures. Unlike conventional machining processes, 
which remove material to obtain a desired shape, RP systems are additive in nature. 
A part is built progressively, layer by layer. This technology has also been referred to 
as layered manufacturing technology, free-form fabrication, model making, desktop 
manufacturing, etc. 



The layer by layer manufacturing technologies poses several accuracy issues. Some 
of these issues are technical and may be as a result of software and hardware lim- 
itations, while others may be due to finishing operations. There are limitations in 
CAD to RP translation which generally take place in the CAD system itself, for e.g., 
tolerance for tessellation, convex boundary error, flipped normal, mid-line node, etc. 
The dimensions of the delivery systems of RP determine the minimum feature size 
that can be built in the RP system. RP parts usually exhibit a staircase effect on 
slanted and curved surfaces because of the layering process. Areas that overhang the 
prior layer must be supported to prevent their ciurl deformation. Orientation of the 
build is to be determined considering several factors, such as, surface flnish, build 
time, distortion, etc. 

The RP systems consist of mechanisms to move optics, build head, elevator plat- 
form, Z-stage platen, etc. The links of a mechanism are manufactured with some tol- 
erances on the link lengths and clearances at the joints. The tolerances and clearances 
cause mechanical error in the desired position of the tool or the platform. Mechanical 
Error is a significant accuracy issue in RP processes. A search of the literature has yet 
to disclose any work on analysis and synthesis of mechanical error in RP processes. 
The objective of the present work is the analysis and synthesis of mechanical error in 
RP processes. Three RP processes. Laminated Object Manufacturing (LOM), Fused 
Deposition Modeling (FDM) and StereoLithography (SL) are considered. Using the 
stochastic approach, a methodology for analysis as well synthesis for mechanical er- 
rors such as tolerances and clearances of salient feature dimensions of these three RP 
processes have been studied. In order to achieve a desirable performance accuracy 
of an RP process, it is necessary to allocate optimally the levels of tolerances and 
clearances. This aspect is also studied in the present work. 

While designing a mechanism either for function generation or for path generation 
or for rigid-body guidance, it is necessary to take into account not only the struc- 
tural error but also the mechanical error. Mechanical error in mechanisms has been 
dealt considerably in the past literature. Indiscriminate allocation of tolerances and 
clearances may lead to unacceptable deviation of the output. On the other hand, 
the cost of manufacturing and assembling such a mechanism with stringent values 
of tolerances and. clearances may be prohibitive. A designer would, therefore, hke to 
assign tolerances and clearances such that the mechanical error is within the specified 
limi ts and the levels of tolerances and clearances are as high as possible.Tolerances 
and clearances are random in nature. The stochastic approaches have been found to 
be more suitable for both analysis and synthesis of mechanical error. In the present 
work mechanical error in RP processes is studied using the stochastic approach. 

Chapter 1 of the thesis gives an overview of RP processes. It describes three 



RP processes LOM, FDM and SL in brief. The accuracy issues in RP processes are 
discussed. A review of the hterature pertaining to analysis and synthesis of mechanical 
error in mechanisms is presented. 

Chapter 2 of the thesis presents a unified approach for studying the effects of 
manufactinring errors in RP processes. It gives random errors associated with link 
lengths, joints in hnkages, coordinates of a point and orientation of a link. Stochastic 
model of an RP process is discussed. The coordinates of a point on the work surface 
traced by the laser beam or the tip of the extruder head is expressed as a function 
of random variables involved in the process. If there are n random variables Vj, V 2 , 

. . . , Vn involved in the process then the coordinates of the point are given by x(Vi, 
V 2 , , Vn), y{Vi, V 2 , ... ,'Vn) and z{Vi, V 2 , ... , Vn). The methodology to associate 

any given range of a dependent variable x, y or z with its corresponding probabihty 
is given. For normally distributed random variables, the dependent variable may as 
well be taken as normal for more than five random variables. For a given probability, 
the range of a dependent variable is evaluated from the variance of the dependent 
variable. The variances of the dependent variables x, y and z are expressed in terms 
of the influence coeflSicients (0:^ , (^)^‘ 

Chapter 3 performs analysis of mechanical error in LOM. The LOM process cuts 
and binds foil . Geometric model of the LOM process consists of rollers, work surface, 
location of the laser source, etc. The geometric model allows to express a point on 
the contour of a shce on the LOM sheet. The work surface and the laser beam is 
expressed parametrically and the intersection of the laser beam and the work surface 
is found. Newton-Raphson method is used to solve the transcendental equations. 
Expressions for the influence coeflicients and the partial derivatives needed for them 
are derived. A numerical example is given for a set of input values. The variances 
D[x], D[y] and D[z], and their sum are hsted for a grid of points on the work surface. 
The three-sigma bands of error in tracing several hues on the work surface is plotted. 
This is the band with probability 0.9973. 

Chapter 4 performs analysis of mechanical error in FDM. The FDM process com- 
prises of an extruder head moving in XY-plane. The tip of the extruder head deposits 
material on the Z-stage platen. It draws the contour of a slice on the platen and fills 
the contour with roads. Geometric model of the FDM process consists of links, hinge 
joints, location and orientation of Z-stage platen, etc. The expression for a point on 
the contour of a slice traced by the nozzle tip is found. The influence coefficients and 
the partial derivatives needed for them are derived. A numerical example is presented 
for given tolerances and clearances, and other input values. The variances and their 
sum are listed for a grid of a points traced by the nozzle tip. The three-sigma bands 
of error in tracing few example curves is plotted. 


VI 



Analysis of mechanical error in SL is presented in Chapter 5. In the SL process, 
the laser beam moves across the surface of a photocurable liquid acrylate resin in a vat 
and traces the contour of a shce. Geometric model of the SL process consists of the 
location of the laser source, orientation of the mirrors, location and orientation of the 
elevator platform, etc. Parametric representation of the ray incident on the mirror is 
found. For a desired reflection from the galvanometer-driven mirror, the orientation 
of the mirror and the intersection of the incident ray and the mirror is found. Newton- 
Raphson method is used to solve the non-linear equations. Expression for a point on 
the contour of a slice is drawn by the laser on the resin surface is derived. A numerical 
example for a set of input values gives the list of variances. The three-sigma bands 
of error for few example curves are also plotted. 

Chapters 2, 3, 4 and 5 deal with the analysis of mechanical error in RP processes. 
The synthesis part is the inverse of the above problem. In synthesis, the designer 
has to decide the levels of tolerances on the random variables Vi for certain allowable 
tolerance hmits on the dependents variables x, y and z. Chapter 6 deals with the 
optimal allocation of tolerances and clearances in RP processes. Since the constraints 
axe linear functions of design variables, therefore, the optimization is done using 
Lagrange Multipher Technique. Optimal allocation of tolerances and clearances is 
done in the LOM, FDM and SL processes. Optimization is also done using Genetic 
Algorithm, GA in short, for FDM and SL. The optimum values of variances of random 
variables obtained using Lagrange method and GA are compared. 

The conclusions of this work are presented in Chapter 7. The scope of further 
work has also been discussed in this chapter. 



ACKNOWLEDGMENTS 


I would like to express my deep sense of gratitude to my thesis supervisor Prof 
Sanjay Govind Dhande for his guidance over the past few years and for his patience 
during the progress of the work. I would also like to thank him for his interest and 
invaluable suggestion during the preparation of this document. 

I express my sincere gratitude to Prof Kalyanmoy Deb for his suggestions dming 
the synthesis part of the problem and for his encouragement during the progress of 
the work. His suggestions and explanations in using Lagrange multipher technique 
and real coded genetic algorithm were of great help. I gratefully acknowledge him for 
providing me the code for GA in C language. 

I am grateful to Prof A K Mallik, Prof B Sahay, Prof Prashant Kumar, Prof K 
K Saxena, Prof N N Kishore, Prof Vijay Gupta, Prof Sudipto Mukherjee, Prof S K 
Choudhry, Prof P S Ghoshdastidar, Prof Gautama Biswas, Dr Amitabha Mukerjee 
and all the faculty members of the institute for the courses they taught me and for 
their moral support during this work. I thank Dr K C Goel and Dr C M Singhal for 
the medical care provided to me. 

I express my sincere thanks to Mr A D Bhatt for the support and encouragement 
during our dissertations. I thank Shiva Prasad who was an instant help during my 
computation. I thank Saurabh Vishal for the moments shared together during the 
final phases of our respective theses. 

I would be filled with nostalgia when I remember the light moments shared with 
Atulji, Charu bhabhi and children. Their company and their constant encouragement 
kept us cheerful. I enjoyed the company of Mahesh Gupta and Kushal Qanungo during 
sports and leisure time in hall V. I gratefully acknowledge their support during my 
dissertation. 

I wish to thank all CAD lab members KVND Ramesh, Subhash Mishra, PVM 
Rao, Jakka Venkatesh, Palanisami, Kumaraswami, Rahul, Bansal, M Anil, Manoj, 
Siddharth, Anil, Puneet, Mukul, Sarvanakumar and M L Jain for their help in various 
ways during the course of my thesis work. My sincere thanks to the staff of CAD 
Project Laboratory for their help and cooperation during the thesis work. 



I thank my friends in the campus and outside for their moral support. I thank 
Sudhir, Satya Prakash, Dwarika, Ashutosh Mishra, Jugal, Mini Ghosh, Sanjay Ku- 
mar, B N Singh, Jilaniji, Anupam Gupta and Rajiv Awasthi. 

I express my deep sense of gratitude to my wife Ranjana for her cooperation and 
endurance during the thesis work. Her constant support and encouragement while 
she herself was completing her dissertation will always be cherished. 

Finally, I wish to express my heartfelt thanks and gratitude to my parents, brother, 
sister and all family members for their constant support and encouragement during 
the work. Ay, 

(Sanat Agrawal) 



Contents 


List of Figures xiii 

List of Tables xv 

List of Symbols xvi 

1 INTRODUCTION 1 

1.1 Introduction to Rapid Prototyping 1 

1.1.1 Process sequence of Rapid Prototyping 2 

1.1.2 Overview of RP processes 3 

1.1. 2.1 Laminated Object Manufacturing Process 3 

1.1. 2. 2 Fused Deposition Modeling Process 5 

1. 1.2.3 Stereolithography Process 6 

1.1.3 Performance Parameters of RP Processes 8 

1.2 Analysis of Mechanical Error in Mechanisms: A Brief Review 13 

1.3 Synthesis of Mechanical Error in Mechanisms: A Brief Review .... 16 

1.4 Objective and Scope of the Present Work 17 

2 STOCHASTIC MODELING OF MECHANICAL ERROR 20 

2.1 Introduction 20 

2.2 Random Errors in the Elements of the Geometric Model of an RP Process 21 

2.2.1 Random Error in a Revolute Joint and a Revolute Dyad ... 22 

2.2.2 Random Error in a Spatial Revolute Joint 24 

2.2.3 Random Errors in the Coordinates of a Point 26 

2.2.4 Random Errors in the Orientation of a Link or an Axis .... 26 

2.2.5 Stochastic Model of an RP Process 27 

3 ANALYSIS OF MECHANICAL ERROR IN LOM 31 

3.1 Geometric Modeling for Description of the LOM Process 31 

3.1.1 The Transformation Matrices [Tq], [Tq] and [T^] ........ 34 

3.1.2 Mathematical Representation of the Work Surface ....... 36 



3.1.3 The End Points of the Edge AC of the Work Surface 37 

3.1.4 Solving for u and r by Newton-Raphson Method 40 

3.1.5 The Four Corner Points of the Work Surface in the Base Erame 

XYZ 41 

3.1.6 Parametric Representation of the Laser Beam 42 

3.1.7 Intersection of the Laser Beam and the Work Surface 42 

3.2 Stochastic Error Analysis of the LOM Process 46 

3.2.1 Partial Derivatives of the Elements of [Tf] 50 

3.2.2 Partial Derivatives of xi, yi, X 2 , y 2 and Z 2 54 

3.2.3 Partial Derivatives of the Elements fij 56 

3.2.4 Partial Derivatives of the Parameters u and r 58 

3.2.5 Partial Derivatives of the Elements of [Tq] 59 

3.2.6 Partial Derivatives of the Elements of [Tq ] 61 

3.2.7 Partial Derivatives of the Elements of [p] 62 

3.2.8 Partial Derivatives of u and w 63 

3.2.9 Partial Derivatives of the Elements of [T] . , 65 

3.2.10 The Influence Coefficients 

3.3 Numerical Results 66 

4 ANALYSIS OF MECHANICAL ERROR IN FDM 72 

4.1 Geometric Modehng for Description of the FDM Process 72 

4.2 Stochastic Model of the FDM Process 75 

4.2.1 Mean Values of the Random Variables 78 

4.2.2 Variances of the Random Variables 78 

4.2.3 The Angle ja 79 

4.2.4 The Angle 7J, 81 

4.3 The Transformation Matrix [Tq] 82 

4.4 The Transformation Matrix [Tq ] 82 

4.5 The Point Q on the Contour of a Slice Traced by the Nozzle Tip ... 85 

4.6 The Means and the Variances of the Dependent Variables x, y and z . 86 

4.7 Partial Derivatives of the Elements of [X4] 86 

4.8 Partial Derivatives of the Elements of [Tjg] 87 

4.9 The Influence Coefficients j (•^) and 88 

4.10 Numerical Results 89 

5 ANALYSIS OF MECHANICAL ERROR IN SL 94 

5.1 Geometric Modehng for Description of the SL Process 94 

5.2 Stochastic Model of the SL Process 96 



5.3 The Incident Ray EF 97 

5.4 The Transformation Matrices [Tq] and [T® ] 98 

5.5 Solving for the Mirror Plane and the Intersection Point F 99 

5.6 Orientation of the Frame X 2 Y 2 Z 2 Using Direction Cosines of X 2 and 

Z 2 axes 106 

5.7 Ejqiression for tf , the Value of the Parameter t at the Point F . . . . 110 

5.8 The Reflected Ray FQ Ill 

5.9 The Point Q on the Contour of a slice Drawn by the Laser on the Resin 

Surface 112 

5.10 The Means and the Variances of the Dependent Variables x, y and z . 113 

5.11 The influence coefficients 114 

5.11.1 Partial Derivatives of [wl_sq] 114 

5.11.2 Partial Derivatives of sg 115 

5.11.3 Partial Derivatives of [fp2], [f 2 ], [r2_t/j, [rjtf\ and tf . . . 115 

5.11.4 Partial Derivatives of [ef, [62] and [e] 115 

5.11.5 Partial Derivatives of the Elements of [Tf] and [t21_3] 116 

5.11.6 Partial Derivatives of the Elements of [Tf ] and [t02_3] 118 

5.12 Numerical Results . 119 

6 OPTIMAL ALLOCATION OF TOLERANCES AND CLEARANCES 123 

6.1 Optimization Using Lagrange Multiplier Technique 124 

6.2 Optimization Using Genetic Algorithms 128 

6.3 Optimum Tolerances in LOM 129 

6.4 Optimum Tolerances and Clearances in FDM 131 

6.5 Optimum Tolerances in SL 133 

7 CONCLUSIONS 137 

7.1 Conclusions 137 

7.2 Scope for Further Work 138 


References 


140 



List of Figures 


1.1 Process Sequence of RP Processes 2 

1.2 A Schematic Diagram of the LOM Process 4 

1.3 A Schematic Diagram of the FDM Process 5 

1.4 The Schematic Diagram of Stereolithography Process 7 

1.5 Tolerances in Faceted Formatting 8 

1.6 Convex Boundary Error Issue 9 

1.7 Mid-Line Node 10 

1.8 Closure Error Causing Hole in Model 10 

2.1 An RP Process 21 

2.2 A Revolute Dyad 22 

2.3 Stochastic Model of a Revolute Joint 23 

2.4 Stochastic Model of a Revolute Dyad 24 

2.5 Clearance Model for a Revolute Joint in a Spatial Mechanism .... 25 

2.6 The Link AB in XYZ Frame 27 

3.1 Geometric Model of the LOM Process 32 

3.2 The Cyhndrical Surfaces Pi (u,u) and q 2 (r, s) 33 

3.3 Transformations [Tq] and \r^] 34 

3.4 The Coordinate Frames AF.Z and A'lliZ'i 35 

3.5 Mathematical Model of the Work Surface 36 

3.6 The Bihnear Surface ABDC in XYZ Object Space 37 

3.7 Tangent Points A and C in the frame X\YiZi 38 

3.8 The Incident Laser Beam 42 

3.9 Intersection of the Laser Beam and the Work Surface 43 

3.10 Geometric Model of the LOM Process Showing Errors in Dimensions 47 

3.11 The 3cr band for a Horizontal Line in uw Parametric Space from (.l,.l) 

to (.9,.l) 70 

3.12 The 3(7 band for a Vertical Line in uw Parametric Space from (.9,.l) 

to (.9, .9) 70 



3.13 The 3(7 band for an Inclined Line in uw Parametric Space from (.l,.l) 

to (.9, .9) 71 

4.1 Motion of Extruder Head in FDM 73 

4.2 Geometric Model of the FDM Process 74 

4.3 Equivalent Kinematic Chain 76 

4.4 The Angles 7^ and 75 80 

4.5 Angle 7 Made by a hne OP With the XZ-plane 81 

4.6 The Transformation [Tq] 83 

4.7 The Transformation \Tq] 84 

4.8 The 3cr band for Tracing a Curve by the Nozzle Tip whose Projection 

in AoTo“P1^3Qe is an Inclined Line from (.05, .05) to (.20, .20) 92 

4.9 The 3cr Band for Tracing a Curve by the Nozzle Tip whose Projection 

in AoFb-Plane is a Circle With Center at (.15, .10) and Radius .05 m 93 

5.1 Geometric Model of the SL Process 95 

5.2 Geometric Model of the SL Process Showing Error in Dimensions . . 96 

5.3 The vector EF • 98 

5.4 The Transformations and [T®] 99 

5.5 The Mirror Plane and the Plane of Reflection 101 

5.6 Four Possible Solutions to Unit Normal Vector on the Plane of Reflectionl04 

5.7 Direction of X 2 Axis 107 

5.8 Directions of X 2 and I2 Axes 108 

5.9 The Reflected Ray FQ Ill 

5.10 The 3cr band for an Inclined Line in Object Space from (.01, .28) to 

(.28, .01) 122 

5.11 The 3cr Band for an Ellipse With Center at XiYx = (.15, .15), Semi- 
major Axis = .14, Semiminor Axis = .10 and Inclination Angle = 0° 122 

6.1 The region Uj + aibi + a^Ci >0 127 

6.2 Feasible Region of ax and 128 



List of Tables 


3.1 Variances and Their Sum at Several Points in uw Parametric Space on 


the Bilinear Surface ABDC in LOM 68 

3.2 The Influence Coefficients at (u = 1.0, w = 1.0) in uw Parametric 

Space on the Bilinear Surface ABDC in LOM 69 

4.1 Variances and Their Sum at the Nozzle Tip Q in FDM Obtained by 

Varying h and Iz 91 

4.2 The Influence Coefficients for li = .25 and Iz = .25 in FDM 92 

5.1 Variances and Their Sum at Several Points on the Resin Surface ob- 
tained by varying Qx and qy 120 

5.2 The Influence Coefficients at qx = 0.0 m and qy = 0.3 m in SL . . . . 121 

6.1 Optimal Allocation of Tolerances in LOM . 130 

6.2 Optimum Values of Variances pi in FDM from Lagrange Multiplier 

Technique and from Real Coded GA 132 

6.3 Optimal Allocation of Tolerances in FDM 133 

6.4 Optimum Values of Variances pi in SL from Lagrange Multiplier Tech- 
nique and from Real Coded GA 135 

6.5 Optimal Allocation of Tolerances in SL 135 



List of Symbols 


influence coefficient 

bi influence coefficient 

a influence coefficient 

dj radial clearance in the pair ij 

Cjk covariance of the dependent variable Vj and 

DSUM sum of the variances of the dependent variable x, y and z 

D[Vi] variance of the random variable Vi 

D[x] variance of the dependent variable x 

E[w] expected value of the dependent variable w 

FNi unit vector normal to the mirror plane GHIJ in SL 

li nominal length of the i*^ link 

l(t) parametric representation of the line AC in LOM in XiYiZi frame 
hfixed length of the portion of link 1 in FDM where there is no motion 
hmax maximum variation of the length of fink 1 

hfixed length of the portion of fink 3 in FDM where there is no motion 
hmax maximum variation of the length of link 3 
m[Vi] mean of the random variable V 
m[x] mean of the dependent variable x 

M modified cost function 

Pi(u, v) parametric representation of a cylindrical surface in XiY^Zi frame 
q(u,w) parametric representation of the bihnear surface ABDC in LOM 


in XYZ frame 



<i2{r, s) 
Q(t) 
r(t) 
[^2-t/] 

Ri 

Rii 

sq 

tf 

[i02_3] 

[i21_3] 

[t;] 

Vi 

( 4 ’ Vij) 

Greek 

7 

€i 

e 

Ai, A2, A3 

pi 

a 

4 > 


parametric representation of a cylindrical surface in XiYiZi frame 

parametric representation of a cylindrical surface in X2Y2^2 frame 

parametric representation of the ray EG in LOM in XYZ frame 

ray EF in the base frame XoVoZo in SL 

position vector [r2(t)] at t =tf in SL 

actual length of the link 

equivalent length of the Hnk 

value of the parameter s at the point Q in SL 

value of the parameter t at the point F in SL 

third column of [T^] in SL 

third column of [Tf] in SL 

homogeneous transformtion matrix of the frame XiYiZi with 
respect to the frame XjYjZj 
random Variable 

postion vector [loi] at s = sg in SL 
coordinates of the pin axis in frame 

Symbols 

angle made by a hnk with a plane perpendicular to its socket 

error on a random variable 

azimuth angle of the laser beam 

lagrange multipliers 

variance of the random variable V 

standard deviation of error at a point on the work surface 
zenith angle of the laser beam 



Chapter 1 


INTRODUCTION 


1.1 Introduction to Rapid Prototyping 

Rapid prototyping is emerging as a key prototyping technology with its ability to 
produce even complicated parts within hours. Prototyping is an essential part of the 
product development cycle required for assessing the form, fit and functionality of a 
design before a significant investment in tooling is made. Untfi. recently, prototypes 
were largely hand-made by skilled craftsmen, adding weeks or months to product 
development time. 

Rapid Prototyping (RP, in short) is a term which embraces a range of new tech- 
nologies for producing accurate parts directly from CAD models, in hours, with little 
need for human intervention. This means that designers have the choice of making 
physical models of the conceptual design more frequently, allowing them to check the 
form and assembly of the design. This also allows them to discuss manufacturing 
issues with an easy-to-interpret prototype which can be held, viewed, studied, tested 
and compared. Consequently, a product is improved and product development cycles 
and costs substantially reduced. It has been claimed that RP can cut new product 
costs by upto 70 % and product development time by upto 90% [95]. 

RP systems take the solid model of a part as an input and fabricate a physical 
model or prototype without using tools or fixtures. Unlike CNC machines, which 
remove material to obtain a desired shape, RP systems are additive in nature. A 
part is built progressively, layer by layer. This technology has also been referred to 
as layered manufacturing technology, free-form fabrication, desktop manufacturing, 
model making, low- volume manufacturing, material incress manufacturing, etc [2]. 
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1.1.1 Process sequence of Rapid Prototyping 

The RP process starts with the creation of a sohd model and ends with a finished 
prototype. The process sequence of rapid prototyping consists of the following four 
tasks: 


• Geometric Modeling 

• Data conversion and checking 

• Part building 

• Postprocessing 

The process sequence is shown in Fig. 1.1. 



Figure 1.1: Process Sequence of RP Processes 


Geometric Modeling 

Creating an accurate and complete solid model is a general prerequisite of RP processes. 
The solid model must be a closed volume. There are several sohd modelers available, 
for e.g., Pro/ENGINEER, I-DEAS, Unigraphics, CADDS 5i and CATIA. Geometric 
models can also be constructed using data from Computerized Axial Tomography 
(CAT), Magnetic Resonance Imaging (MRI), reverse engineering, etc., depending 
upon the need. 
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Data Conversion and Checking 

Geometric data of a model is transferred from a CAD system to an RP system 
primarily through the STL (STereoLithography) file format. This format has become 
a de facto standard for interfacing with RP systems. The solid model is converted to 
STL format. The STL is a faceted format and consists of connected three dimensional 
triangles representing the part shape. The vertices of the triangle are ordered to 
indicate which side of the triangle contains the part mass [62]. Prom a mesh generation 
point of view, triangular meshes are easy to generate than quadrilateral meshes. This 
is because a triangle is a simplex while a quadrilateral is not [133, pp. 930]. Some 
RP processes cannot effectively handle sliced models, (or “2|D” models, as referred 
to in hterature by Dolenc and Makela [32] ), and so a faceted model is taken and 
sliced again. This is because the sliced model cannot be positioned arbitrarily in the 
workspace of an RP system [32]. 

The translation from CAD to RP resxilts in some errors. The typical errors are 
flipped normals, mid-line node, closure errors (holes) and truncation errors [36]. The 
validity of the STL file is verified using a correction software and the file is repaired. 

Part Building 

This process starts with placing the STL file for the optimum part orientation, and 
extracting the shce information from it, checking the slices for any open curves and 
other defects and repairing them, adding the part-building parameters, if any, and 
sending the file thus created to the RP machine. The machine then builds the part 
described by the file, and this process requires no or minimal human intervention. 
The part is thus built by the machine. 

Postprocessing 

In the postprocessing stage, some manual operations are done and these operations 
require expertise of the operator. The manual operations includes removal of supports 
and excess material, if any, and postcuring, if required. 

1.1.2 Overview of RP processes 

1. 1.2.1 Laminated Object Manufacturing Process 

The Laminated Object Manufacturing (LOM) process from Hehsys Inc. cuts and 
binds foils as illustrated in Fig. 1.2. The undersurface of the foil has a thermoplastic 
binder that, when pressed by a heated roller, causes it to glue to the previous foil. The 
foil is cut by a GO 2 laser following the contour of the slice. A computer controlled 



Chapter 1. Introduction 


XY positioning device and a set of mirrors are used to direct the laser beam. To help 
the removal of excess material once the parts have been built, the exterior of the shce 
is hatched, as opposed to the stereolithography process where the interior is hatched. 
This excess material remains in place during the build process to act as support. The 
thickness of the foil is not constant. Therefore, a sensor measures the current foil 
thickness, and the foil is shced accordingly. 



The Helisys product line consists of LOM-1015Plus and LOM-2030H machines. 
Material thickness in LOM-1015Plus is from 0.08 to 0.25 mm (0.003 to 0.008 in) 
and the maximum part envelope size is 381x254x356 mm (15x10x14 in) [55]. Part 
accuracy is ± 0.125 nun (± 0.005 in) anywhere in the work envelope [2, 95, 31]. 
The LOM process can produce parts from paper, plastics, metal foil, composites and 
ceramics in sheet form [112, 63, 2]. 

Advantages of LGM include a wide variety of organic and inorganic materials 
which can be used. Speed is another strong point of LOM. As only the outhnes of 
the parts need to be traced by the laser, this method is very fast. This method is 
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best suited for building large parts. The process produces virtually no internal stress 
and associated undesirable shrinkage and warping [55, 95]. The procedure is not well 
suited for manufacturing parts with thin walls in z direction and hollow parts. The 
parts are strong in x and y directions and have a tendency to delaminate in z direction 
in intricate parts [62, 112]. Other disadvantages of this technology are that there is 
a large amount of scrap, the machine must be constantly manned, parts need to be 
hand-finished and the shear strength of the part is limited by the bonding strength 
of glued layers [95]. 

1.1. 2. 2 Fused Deposition Modeling Process 

In the Fused Deposition Modeling (FDM) process from Stratasys Inc., a spool of 
1.27 mm (0.05 in) diameter filament feeds into the unit’s extruding head. The extruder 
head moves in the XY-plane (See Fig. 1.3). The filament is melted to liquid inside 
the extruder head at 82° C (180°F) by a resistance heater. FDM produces parts from 
sohd CAD models. The model is then converted into STL format in the CAD system 
and sent to the FDM slicing software, called QuickSfice. There the .STL file is sliced 
into thin cross sections of desired thickness, creating a .SLC file. 



Supports are created for overhanging parts and sliced as well. The sliced model 
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and support are converted into a .SML file (Stratasys Machine Language file) that 
contains actual instructions for the FDM machine. A tool path is generated which 
is followed by the numerically controlled extruder head. As the head moves in X 
and Y directions following the tool path, the thermoplastic material is extruded out 
of a nozzle and then deposited in ultra thin layers, one layer at a time. Since the 
envelope surrounding the head is maintained at a temperature below the melting 
point of material, the extruded material quickly solidifies [63, 95]. The extruder head 
has two nozzles, one for the part material and the other for the support material. 
The support can be easily removed by breaking away. The part is built on a foam 
foundation attached on a Z-stage platen. The Z-stage platen moves downwards as 
the part is built progressively. 

The system uses ABS — a tough nylon-like thermoplastic, ABSi — a high im- 
pact ABS, machinable wax, investment casting wax, elastomers, etc., as modeling 
materials. 

The specification of FDM2000 are as follows. The work envelope is 254x 254x 254 mm 
(10x10x10 in). The operator may select a road width between 0.254 and 2.54 mm 
(.010 and .100 in) and a layer thickness between .05 and .762 mm (.002 to .030 in). 
Models can be produced within an accuracy of ± 0.127 mm (± 0.005 in) [114]. 

The FDM system has two position sensors. Theses sensors are to place the ex- 
truder head in the home position before starting the part build. They do not sense 
the position of the head during part build. 

This system may be viewed as a desktop prototyping facility in a design office since 
the modehng materials are cheap, non-toxic, odorless and environmentally safe. A 
wide range of material is available in various colors. Parts made by this method have 
a high stability since they are not hygroscopic. The system doesn’t waste material 
during or after producing the model. No postcuring of model is required [13]. A 
drawback of the FDM process is that the surface finish of a part is poor relative to 
other processes. 

1. 1.2,3 Stereolithography Process 

The StereoLithography Apparatus (SLA) system from 3D Systems Inc. is the most 
representative system running on StereoLithography (SL) process. The SL process 
starts with a solid CAD model of an object that is downloaded to a sficing algorithm, 
which cuts the CAD model into many thin layers. The part itself is built in a vat of 
liquid acrylate or epoxy resin that sohdifies when the laser light hits it (See Fig. 1.4). 

An elevator platform rests below the liquid surface at a distance equal to the thickness 
of the first bottom-most layer of the part to be built. An ultraviolet laser beam from 
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the HeCd-laser is deflected by the galvanometer-driven mirrors and scribes the contour 
of a shce of the part on the resin surface. A second, HeNe-laser is used to ensure that 
the surface of the liquid is in correct location [60, 31]. 



Figure 1.4: The Schematic Diagram of Stereolithography Process 

The interior of the contour is then hatched using a hatch pattern. After a layer is 
built, the elevator platform is lowered a programmed distance so that a new coating 
of hquid resin covers the solidified layer. A recoater blade then sizes the coating 
depth. The laser draws a new layer on top of the first layer and the part is built 
layer by layer from the bottom up. Finally the part is removed from the vat and 
placed in a curing oven, which floods the object with ultraviolet fight to complete the 
polymer solidification process [2, 95]. For parts with overhanging geometry, support 
structures are created that attach to the elevator platform to support the part while 
it is being built [63, 62]. Postprocessing for the SL process include removal of excess 
resin, postcuring, and surface finishing as required. 

The SL product line consists of SLA 190, SLA 250, SLA 500, SLA 3500, SLA 5000 
and SLA 7000. The SLA 250/50 specifications are as follows. The layer thickness 
can be selected between .063 and .76 mm (.0025 and .03 in) [60, 74]. The minimum 
recommended layer thickness is .1 mm (.004 in) [1]. The laser beam spot size is .20 - 
.28 mm (.008 — .011 in). The build envelope size is 250x250x250 mm (lOx lOx 10 in). 
The vertical resolution of the above platform is .0025 mm (.0001 in). Platform re- 
peatability is .013 mm (.0005 in). Scan resolution and scan repeatability are .008 mm 
(.0003 in) and .130 mm (.005 in) respectively. Part accuracies are ± .127 mm (± 
0.005 in) [2, 95]. 

The SL process produces a surface finish that is comparable to that of NC milling. 
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The process is reasonably fast and accurate [95]. A part with intricate features can be 
made with virtually no hmitations. The area of apphcation of the process is restricted 
due to expensive, smelly and toxic material. 

1.1.3 Performance Parameters of RP Processes 

The layer by layer manufacturing technologies poses a unique set of issues, irrespec- 
tive of the techniques involved in the fabrication of layers. Some of these issues are 
technical and may be as a result of software and hardware hmitations, while others 
may be due to the finishing operations. There are hmitations in CAD to RP transla- 
tions which generally take place in the CAD system itself. The best output that can 
be obtained from a machine driven by the CAD software cannot be better than the 
accuracy used in the CAD system. Several process parameters of RP which affect the 
accuracy of RP processes are as follows: 

(i) Accuracy Issues in CAD to RP Translations 
(a) Tolerfince for Tessellation 

As discussed earlier, the tessellation process is the translation of the surfaces of 
CAD model into triangular facets. In order to generate the faceted format the 
CAD system asks the user the acceptable tolerances. These tolerances can be 
the largest distance from the true surface to the triangle, or the largest distance 
from the true surface to a side of the triangle. The latter is termed as chordal 
tolerance (See Fig. 1.5). The smaller the tolerance, and more irregular the 
surface is, the larger the number of triangles required to represent the surface. 
The RP technologies has spot sizes between .127-2.5 mm (.005 - .100 in) [36]. 
Therefore, very small tolerances become meaningless. Large number of facets 
increase the .STL and slice file size and the user has to strike a balance between 
accuracy and file size. 



Figmre 1.5: Tolerances in Faceted Formatting 


(b) Convex Boundary Error 
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The tessellation results in a smaller object where the object is convex, whereas 
it results in added material where the object is concave (See Fig. 1.6). Since 
most RP systems are typically sanded to obtain good surface finish, the lo- 
cations short on material will be even more distant from the original desired 
surface, whereas the locations with added material can be sanded to accurate 
dimensions. 



(c) Flipped Normals 

Occasionally, the computed normal of a triangle obtained from translation is in 
the opposite direction, or is just plain wrong, and must be corrected. They are 
corrected by using some heuristic rules on adjoining triangles, or simply based 
on a new cross product [36]. 

(d) Mid-Line Node 

The basic rule to check the .STL file is that each edge of a triangle is shared 
by two triangles. Figure 1.7a illustrates an erroneous triangularization and 
Fig. 1.7b, a correct one. The case in Fig. 1.7a results in an incorrect model of 
the object and in a failure of the slicing algorithm. 

(e) Closure Errors 

Round-off errors during tessellation cause one point to be at multiple locations 
at the same time. Thus triangles are formed, and a thin hole is present in the 
finished model (See Fig. 1.8). In order to correct such an error, points that are 
within an infinitesimal distance, say delta, are considered, and the triangles they 
individually form are checked to see whether they share a common edge. It is to 
be noted that for a very small chordal tolerance, two legitimate points are very 
close to each other and it becomes difficult to differentiate between legitimate 
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Figure 1.7: Mid-Line Node 

points and round-off error. Round-off errors also result in overlapping facets, 
and non-manifold edges, points and faces [63]. 



Figure 1.8: Closure Error Causing Hole in Model 


(ii) Laser Beam or Nozzle Diameter 

The dimensions of the delivery systems of RP deternaine the minimum feature 
size that can be built in an RP system. For laser based systems, the minimum wall 
thickness in a part would be the width of a single cmed line. This is related to laser 
beam diameter and cure depth. The laser beam diameter used in SL is typically 
.254 mm (.01 in). For a photomask system such as Solid Ground Curing by Cubital 
Inc., the minimum feature size would be related to the resolution of the mask. A 
resolution of 200 dot per inch would result in a minimum feature size of .127 mm 
(.005 in). The minimum road width in FDM1650 for T12 tip is .305 mm (.012 in). 

The laser beam in SL systems is typically at some height above the part, and the 
part can be located at any point in the vat. In SLA 250, the laser is directed at the 
resin from a point 500 mm (20 in) above the surface of the vat and draws on an area 
of 254x254 mm (10x10 in). This causes the deviation of the beam by around 15° 
at the edges of the part. This leads to a part that is shingled and the surface finish 
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quality suffers significantly. For a typical layer thickness of .254 mm (.01 in) and an 
overcure of .1524 mm (.006 in), the variation in dimension of a part at the extreme 
corner of platform is .127 mm (.005 in) [36]. 

(iii) Staircase Effect and Adaptive Slicing 

Parts usually exhibit a visible staircase effect on slanted and curved surfaces be- 
cause of layering process. It causes the smrface finish quality to suffer significantly. 
Adaptive slicing remedies this error by allowing the part to be build at variable layer 
sizes. Orientation of the part can also be occasionally used to mitigate this error. 
Note that uniform sHcing may lead to omission of flat areas (areas parallel to slicing 
direction) and peaks from the part. Dolenc and Makela [32] have presented an adap- 
tive method handhng fiat areas and restricting the staircase effect to a user-specified 
tolerance. Kulkarni and Dutta [73] have given a procedure for adaptive slicing to 
redress the staircase effect and containment problem. Jamieson et al [61] have pre- 
sented a methodology to develop direct and adaptive slicing of CAD models for HP 
processes. Sabourin et al [105] have performed adaptive slicing using stepwise uniform 
refinement. This method is well suited for execution in parallel computer. Tata et al 
[118] have presented an adaptive shcing algorithm that can vary the layer thickness 
considering the geometry of the part. Tyberg and B0hn [122] present a new approach 
to adaptive shcing that significantly reduces fabrication time. 

Hope et al [57, 58] have presented an adaptive slicing method for improving the 
surface finish and decreasing build time. They used sloping boundary surfaces rather 
than stepped edges. 

(iv) Supports 

A majority of the RP processes require the part geometries to be supported as 
the part is being built. In contrast to conventional machining methods, RP processes 
do not require part specific fixtures and toohng, and the process is fuUy automated. 

In the SL process, base supports are built on the elevator platform and the part 
building starts on them. It ensmes that the recoater blade does not strike the platform 
upon which the part is being built. Additionally, it allows for the safe removal of the 
part from its platform and thus preserves its surface finish. 

Areas that overhang the prior layer must be supported to prevent their deforma- 
tion, known as curl. An overhang that extends beyond 1.27 mm (.05 in) will exhibit 
curl when unsupported [60]. In SL, arches or convex surfaces support themselves as 
long as the overhang between successive layers is small. “Islands” are layers of the 
part geometry that would otherwise be unconnected to any other section of the part. 
Islands must be anchored to the platform or to the part itself. 

(v) Part Orientation 
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Several factors contribute to the determination of the orientation of a build. Some 
of these factors are: surface finish, build time, shrinkage, curhng, distortion, support 
structures, roundness/fiatness, resin flow, material cost and trapped volume [93]. The 
surface finish is usually the most important issue. Higher resolution of curved surfaces 
are obtained by orienting them in such a way that their normals lie in the horizontal 
plane. The build time, which depends upon the number of layers involved, can be 
minimized by careful selection of the part orientation. 

Cheng et al [18] determine the optimal part orientation for multi-objectives such 
as part accuracy and build time. Xu et al [129] have obtained optimal the orientation 
with adaptive slicing in SL process taking into consideration build time, accuracy and 
stability of the part. Xu et al [130] have given an algorithm to determine the optimal 
part orientation for different RP processes considering accuracy, surface finish, build 
time and cost. 

There are several works on other accuracy issues, for e.g., part distortion, bench- 
mark component, part accvuacy, etc. Weidermann et al [125] investigate the influence 
of material and process parameters on part distortion to facilitate the selection of 
resins in the SL process. Lee et al [82] consider the accuracy in vertical direction 
in powder-based RP processes. Here the unsintered powder in layers, which act as 
support, undergo downward displacement due to their weight. 

Ippohto et al [59] deal with the development, manufactme and testing of a bench- 
mark component to investigate the dimensional accuracy and surface finish. Childs 
and Juster [19] have designed a benchmark component to assess not only linear ac- 
curacies but geometric tolerances, repeatability and resolution limits of several RP 
processes. Lynn-Charney and Rosen [84] have presented an empirical model for the 
accuracy of a part in the SL process as specified by a set of tolerances on the part. 
They have used a process planning method based on response surfaces that capture 
the relationship among part surfaces, tolerances, and process variables. 

(vi) Machine Accuracy 

The accuracy of a fabricator is its ability to faithfully generate a particular geo- 
metrical shape and size. Let us use the term “tool” for a curing or sintering laser 
beam, material-stream from extrusion or particle-jet nozzle in an RP system. Machine 
accuracy relates to how close to a predetermined position the RP machine can place 
its tool [14]. Accuracy is partially a statistical phenomena. A set of measurements, 
no matter how closely done, always involve some variation. This occurs because the 
factors that affect accuracy and measurement include random effects that cannot be 
ehminated and cannot be controlled. 
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Machine accuracy is different from part and process accuracy. Machine accuracy 
of an RP process involves the following matters: 

• Repeatability of the laser beam / build head position 

• Functioning of machine calibration algorithm 

• Functioning of machine levelling system 

• Functioning of the scanning system 

• Correctness of machine instructions issued by the software 

• Repeatability of the elevator platform / Z-stage platen 

There are mechanisms in RP systems to move optics / build head and elevator 
platform / Z-stage platen. Tolerances on the link lengths and clearances in the joints 
of the mechanisms cause mechanical error in their desired position. Mechanical error 
affect the accuracy and repeatability of the laser beam / build head position and 
elevator platform / Z-stage platen position. The objective of the present work is to 
find the mechanical error in RP processes. Three RP processes, LOM, FDM and SL 
are considered. Mechanical error in mechanisms has been dealt with considerably in 
the past literature [43, 28, 117, 30, 86, 102, 80]. A search of the hterature has yet to 
disclose any work on analysis and synthesis of mechanical error in RP processes. 

1.2 Analysis of Mechanical Error in Mechanisms: 
A Brief Review 

Clearances in the joints and tolerances on the fink lengths of linkages may cause 
mechanical error of appreciable magnitude. While designing a mechanism either for 
function generation or for path generation or for rigid-body guidance, it is necessary 
to take into account not only the structural error but also the mechanical error. By 
assigning stringent values of tolerances and clearances, one may get actual tolerances 
quite close to the desired one. However, the cost of manufacturing and assembling 
such a mechanism may be prohibitive. On the other hand, indiscriminate allocation 
of tolerances and clearances may lead to unacceptable deviation of the output. A 
designer would, therefore, like to assign tolerances and clearances such that the me- 
chanical error is within the specified hmits and the levels of tolerances and clearances 
axe as high as possible. 

Several attempts have been made to analyze path and function generating mech- 
anisms. There are two distinct approaches — deterministic and stochastic. The 
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deterministic approaches are based upon, worst-case analysis of individual tolerances 
[121, 53, 68, 75, 34, 5, 98, 107, 16, 39, 37, 83, 22, 120]. Deterministic approaches give 
highly conservative estimates and do not reflect the overall behavioirr of the mecha- 
nisms. Besides in most of these studies, either tolerances [53, 75, 65, 16, 39, 37, 22] 
or clearances [121, 68, 120] are considered. In contrast stochastic approaches have 
been found to be more suitable for both analysis and synthesis of mechanical error 
[43, 28, 117, 30, 86, 102, 80, 87]. There are several other statistical approaches by 
[43, 21, 10, 6]. 

The effects of both tolerances and clearances have been considered in [28, 86, 87, 
80, 83]. Rao [98] considered elastic deformation along with tolerances and clearances, 
whereas Choubey and Rao [21] and Sutherland and Roth [117] considered structural 
error along with tolerances only. 

The error analysis by Tuttle [121] is perhaps the earhest treatment of the subject. 
He gave a treatment for the estimation of average error due to clearance effects alone, 
but geometry of linkages did not enter into the analysis. Knappe [65] has given an 
approach which permits calculation of the maximum output tolerance as well as the 
relative effect of each individual tolerance on the output. Hartenberg and Denavit 
[53] estimated the mechanical error on function-generating four-bar linkage on the 
basis of maximum allowable tolerance on each link separately; and the maximum 
error was shown to be the sum of errors due to individual links. Garrett and Hall [43] 
have applied a statistical method for mechanical error analysis of function-generating 
mechanisms. The effects of manufacturing tolerances and clearances are presented in 
the form of mobility bands for the hnkages. The effects of clearances on mechanical 
error was taken into account by Kolhatkar and Yajnik [68]. They treated joint clear- 
ances as “equivalent clearance hnks” and determined the maximum output error in 
four and six bar function-generating mechanisms. They drew a number of simple con- 
clusions for four-bar hnkage. Lakshminarayana and Narayanamurthi [75] have shown 
that for complex mechanisms it is convenient to base the analysis on individual loop 
closure equations rather than on complex input-output relationship for the complete 
linkage. The total output deviation is obtained by adding the absolute values of the 
deviations due to individual hnk lengths. 

Dhande and Chakraborty [28, 15] used a stochastic model for mechanical error 
analysis of function-generators. Dubowsky et al [34] utihzes optimization technique 
to identify the dimension which contributes to the mechanical error. The probabilistic 
analysis by Baumgarten and van der Werff [10] utilizes a finite-element description of 
the kinematic chain. A band of error which would contain 99% of the coupler path of a 
population of randomly assembled mechanisms is plotted. Chatterjee and Mallik [16] 
derived the equations for the position error of a coupler curve as functions of crank 
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angle, the link lengths and the error in link lengths. Malhk and Dhande [86] used 
a stochastic approach for mechanical error analysis of path-generating mechanisms. 
Fenton et al [39, 22] have presented a method for error analysis of multi-loop planar 
mechanisms. Lin and Chen [83] presented a generalized error analysis methodology. 
They presented an error matrix to account for the effects of error motions due to 
joint clearances. Ting and Zhu [120] presented a kinematic model and applied N-bar 
rotatabihty laws to identify the worst position and orientation errors due. to the joint 
clearances of linkages and manipulators. 

Several authors have considered sensitivity of mechanisms to the tolerances and 
clearances. Coit and Riley [23] have analyzed the sensitivity of the output of an 
inverted slider-crank straight line generator to the link length tolerances. A link 
length ratio at which tolerances have the most desirable effects on the linkage can 
be selected. Faik and Erdman [37] have proposed a method to synthesize a four- 
bar Hnkage with a prescribed sensitivity value using synthesis solution space. They 
compute the maximum allowable tolerances on link lengths to achieve the desired 
accuracy on link angles. Lee and Gilmore [80] presented the “ effective link length” 
model which considered tolerances on link length, radial clearance and pin center 
location as random variables. They presented a method to determine the sensitivities 
of the output variables with respect to the random variables. 

There are several other approaches to consider the effect of tolerances and clear- 
ances on mechanical error in mechanisms. Sukhija and Rao [115] have apphed in- 
formation theory for optimum synthesis of path-generating mechanisms. The theory 
utihzes the concept of entropy, a measure of uncertainty. Optimum allocation of tol- 
erances on link length is also carried out. Ting and Long [119] used Taguchi’s concept 
to determine the performance quality and tolerance sensitivity of mechanism in the 
early stages of design synthesis. They introduced a Jacobian matrix to relate the per- 
formance variation and dimensional tolerances. Rao and Hati [99] have considered 
the design of a spherical four-revolute function-generating mechanism, with the ob- 
jective of minimizing structural error and manufacturing cost using game theory. Lee, 
Gilmore and Ogot [81] have considered dynamically driven planar kinematic chains, 
either open-loop or closed-loop, with uncertainties in tolerances and clearances. They 
presented probabilistic models and method to allocate tolerances on link lengths and 
radial clearances. Shi [108] has used the reliability concept for synthesis of mechanical 
error in spatial hnkages. 

There are several works on the error analysis of cam, journal bearings, spatial 
mechanisms and other mechanisms [27, 29, 124, 116, 49, 51, 30, 106, 4]. Weihang and 
Qixian [126] have established the stochastic models of spatial joints by introducing 
the concept of “clearance characteristic element” and “clearance space”. Probability 
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distribution and stochastic characteristic of accumulated kinematic error are analyzed 
using Monte Carlo simulation. Yufeng and Jiangqin [132] used matrices to obtain 
tolerances on structural and kinematic parameters of robot manipulators. Zhu and 
Ting [134] presents probabilistic analysis to examine the performance uncertainty due 
to joint clearances in robots. Hati and Rao [54] and Wu and Dalai [128] have given 
their approach to find optimum machining conditions. 

Many researchers have considered optimal tolerance allocation in assembly [111, 
94, 88, 89, 7, 78, 79, 47]. Balling et al [6] have done worst-case analysis of tolerances 
in design optimization. Graves [46] has reviewed four basic tolerancing formulae that 
have appeared in the literature and shown that each solves a slightly different problem. 
He has focused on the more general problem of determining the assembly tolerance 
and component tolerance simultaneously to minimize the total cost of assembly and 
use of a product. 

Several attempts have been made for optimum design of mechanisms [20, 42, 91, 
92, 41, 50, 96, 17, 9, 70, 76, 77]. Root and RagsdeU [104] have given a survey of 
optimization techniques applied to the design of mechanisms. Barker and Baumann 
[8] developed the concept of four-bar solution space. Optimization methods for engi- 
neering design is discussed in [12, 35]. 

1.3 Synthesis of Mechanical Error in Mechanisms: 
A Brief Review 

A significant development in the treatment of tolerances and clearances is their opti- 
mal allocation so as to restrict the mechanical error within specified limits. Dhande 
and Chakraborty [28, 15] used a stochastic model and an equivalent linkage model to 
allocate tolerances and clearances in four-bar function-generators. Sutherland [117] 
has given a method for synthesis of mechanism taking into account structural and 
mechanical error due to tolerances. The dimensions and tolerances of a mechanism 
can be obtained for a given maximum allowable function generating error while min- 
imizing the manufacturing cost. Bakthavachalam and KimbreU [5] have considered 
synthesis of four-bar path-generating mechanisms as an optimization problem under 
inequahty constraints. Equality constraints are modified by introducing tolerances 
and clearances and thereby, the difficulty in satisfying the equality constraints ex- 
actly are eliminated. By this modification, the objective function is also changed. 
The penalty . function approach is used. 

Rao [98] has suggested an iterative method for the synthesis of mechanism taking 
into account the effect of link deformations, tolerances and joint clearances. However, 
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tolerances and joint clearances were specified prior to the synthesis of mechanism. 
Choubey and Rao [21] have suggested a method for minimizing the structural er- 
ror together with the mechanical error due to manufacturing tolerances on the link 
dimensions. Nominal link length are obtained prior to tolerance allocation. The 
mechanical error is treated as a deviation of structural error. Tolerances are then al- 
located with reference to the position of maximum error to hmit the mechanical error 
below a specified value. Sharfi and Smith [107] proposed a method for tolerance and 
clearance allocation in multi-loop planar mechanisms based on the output sensitivity 
with respect to the link lengths. Tolerances and clearances are allocated as a result 
of the synthesis of mechanical error. 

Mallik and Dhande [86] developed a stochastic model for the synthesis of mechan- 
ical error in four-bar path-generating linkages. They analyzed the mechanical error in 
the path of a coupler point for the three-sigma band of confidence level. A synthesis 
procedme to allocate tolerances and clearances so as to restrict the output error in 
the path of coupler point within specified limits is developed. They found that the 
mechanical error of the coupler-point path is dependent on whether one considers 
the original mechanism or its cognate [52, 109] mechanisms. Rhyu and Kwak [102] 
have presented a procedure for optimal stochastic design of mechanisms considering 
tolerances and clearances. A weighted sum of the mechanical error and the manu- 
facturing cost is minimized for the optimal allocation of tolerances and clearances. 
Fenton et al [39, 22] have presented a method for error analysis and tolerance synthesis 
of multi-loop planar mechanisms. 

1.4 Objective and Scope of the Present Work 

In the present work mechanical error in RP processes is studied using the stochastic 
approach. The coordinates of a point on the work surface traced by the laser beam or 
the tip of the extruder head is expressed as a function of random variables involved 
in the process. If there are n random variables Vi, V 2 , . . . , Vn involved in the process 
then the coordinates of the point are given by x{Vi, V 2 , , Vn), V 2 , . . . , Vn) 

and ziy-i, V2, ... , Vn). The methodology to associate any given range of a dependent 
variable x, y or z with its corresponding probability is given. For a given probability, 
the range of a dependent variable is evaluated from the variance of the dependent 
variable. The variances of the dependent variables x, y and z are expressed in terms 
of the influence coefficients (^) ' 

Chapter 1 of the thesis gives an overview of RP processes. It describes three 
RP processes LOM, FDM and SL in brief. The accuracy issues in RP processes are 
discussed. A review of the literature pertaining to analysis and synthesis of mechanical 
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error in mechanisms is presented. 

Chapter 2 of the thesis presents a unified approach for studying the effects of 
manufacturing errors in RP processes. It gives random errors associated with the link 
lengths, joints in hnkages, coordinates of a point and orientation of a link. Stochastic 
model of an RP process is discussed. 

Chapter 3 performs analysis of mechanical error in LOM. The geometric model 
for description of the LOM process is presented. This gives expression of a point on 
the contour of a slice on the LOM sheet. The work surface and the laser beam is 
expressed parametrically and the intersection of the laser beam and the work surface 
is found. Newton-Raphson method is used to solve the transcendental equations. 
Expressions for the influence coefficients and the partial derivatives needed for them 
are derived. A numerical example is given for a set of input values. The variances 
D[x], D[y] and D[z], and their sum are fisted for a grid of points on the work surface. 
The three-sigma bands of error in tracing several fines on the work surface is plotted. 
This is the band with the probability 0.9973. 

Chapter 4 performs analysis of mechanical error in FDM. It begins with the geo- 
metric modeling for description of the FDM process. The expression for a point on 
the contour of a slice traced by the nozzle tip is found. The influence coefficients and 
the partial derivatives needed for them are derived. A numerical example is presented 
for given tolerances and clearances, and other input values. The variances and their 
sum are listed for a grid of a points traced by the nozzle tip. The three-sigma bands 
of error in tracing few example curves is plotted. 

Analysis of mechanical error in SL is presented in Chapter 5. Parametric repre- 
sentation of the incident ray is found. For a desired reflection from the galvanometer- 
driven mirror, the orientation of the mirror and the intersection of the incident ray 
and the mirror is found. Expression for a point on the contour of a slice drawn by the 
laser on the resin surface, is derived. A numerical example for a set of input values 
gives the fist of variances. The three-sigma bands of error for few example curves are 
also plotted. 

Chapters 2, 3, 4 and 5 deal with the analysis of mechanical error in RP processes. 
The synthesis part is the inverse of the above problem. In synthesis, the designer 
has to decide the levels of tolerances on the random variables Vi for certain allowable 
tolerance limits on the dependents variables x, y and z. Chapter 6 deals with the 
optimal allocation of tolerances and clearances in RP processes. Since the constraints 
are linear functions of design variables, therefore, the optimization is done using 
Lagrange multiplier technique. Optimal allocation of tolerances and clearances is 
done in the LOM, FDM and SL processes. Optimization is also done using Genetic 
Algorithm, GA in short, for FDM and SL. The optimum values of variances of random 
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variables obtained using Lagrange method and GA are compared. 

The conclusions of this work are presented in Chapter 7. The scope of further 
work has also been discussed in this chapter. 



Chapter 2 

STOCHASTIC MODELING OF 
MECHANICAL ERROR 


In any electromechanical system the issue of accuracy is important from the 
point of view of successful operation of the system. If the electro-mechanical system 
involves mechanical elements such as cams, linkages, gears, screw and nut assemblies, 
sliding elements, ball bearings, etc., the final performance gets affected due to the 
manufacturing errors in terms of tolerances and clearances of these machine elements. 
In order to evaluate the performance of an electromechanical system, it is necessary 
to study the effect of such manufacturing errors on the final performance of the 
system. The parameters like tolerances and clearances are stochastic in nature. The 
modeling of the process using a set of variables which are stochastic in nature is 
now a well established aspect of design methodology in general and for the design of 
electro-mechanical systems in particular. In the present chapter, a unified approach 
for studying the effects of manufacturing errors in RP processes has been presented. 
Using that unified approach, specific studies have been carried out for three major 
RP processes, namely. Laminated Object Manufacturing (LOM, in short). Fused 
Deposition Modeling (FDM, in short), and Stererolithography (SL, in short) process. 
The details of these studies are given in Chapters 3, 4 and 5. 

2.1 Introduction 

In the LOM process, the sheet material is supphed from a continuous roll on one side of 
the machine and is taken up on the opposite side (Fig. 1.2). The sheet material moves 
between two idle rollers. An XY positioning device contains the optical attachment. 
The laser beam emanates from the optical attachment. The laser beam follows the 
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contour of a slice and cuts the sheet. In the FDM process, the extruder head is 
controlled in the XY-plane (Fig. 1.3). In the SL process, the galvanometer-driven 
mirrors deflect a laser beam emanating from a laser source (Fig. 1.4). The laser beam 
moves across the surface of a photocurable liquid acrylate resin in a vat and traces 
the contour of a slice. 

Geometric model of an RP process consists of several elements involved in the 
process. Geometric model of the LOM process consists of rollers, work sturface, loca- 
tion of the laser somce, etc. Model of the FDM process consists of links, hinge joints, 
location and orientation of Z-stage platen, etc. Model of the SL process consists of 
location of the laser source, orientation of mirrors, location and orientation of the 
elevator platform, etc. 

Tolerances and clearances in these elements are random in nature. These random 
errors are discussed in Section 2.2. 


2.2 Random Errors in the Elements of the Geo- 
metric Model of an RP Process 

The tool of an RP system traces the contour of the shce of a part on a platform. 
There are several mechanisms, comprising of hnks and hinges, responsible for the 
motion of tool and platform. A schematic diagram is shown in Fig. 2.1. The error at 
a point on the contour depends upon the error in the elements of these mechanisms. 
The random errors associated with the link lengths, joints in linkages, coordinates of 
a point and orientation of a link are as follows. 



Figure 2.1: An RP Process 
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2.2.1 Random Error in a Revolute Joint and a Revolute 
Dyad 

Figure 2.2 shows a revolute dyad having links of nominal lengths li and I 2 and pin 
joints 12, 23 and kl. Let €i be the tolerance coefficients defined as the tolerance per 
unit nominal length of the hnk. The actual lengths Ri are given by 

Ri = li €ili i = 1,2 (2-1) 

where Ri are now random variables. 



Figure 2.2: A Revolute Dyad 


There is some clearance always present between the pin of a link i and the race of 
its adjacent link j as shown in Fig. 2.3. Let us choose a two dimensional, Cartesian, 
local coordinate system so that the a;L-axis is coincident with the center line of 
the hnk ending in the race and the origin is at the center of the race. The axis of the 
pin always lies within a circle of radius Vij and center at the origin of the coordinate 
system where rij is the radial clearance of the hinge between and links. 
When the mechanism is in motion the pin axis does not touch the circumference of 
the circle due to hydrodynamic lubrication. Further, if certain foreign particles are 
present in the race or if the whole mechanism vibrates, the pin axis may randomly 
lie anywhere inside the circle. Due to these operating conditions it may be safely 
assumed that the probability that the pin axis lies at a point inside the circle is the 
same for all points [28]. If (xL, y'^j) are the coordinates of the pin axis then the 
probabihty density function is given by 




1 

TrrJ- 

0 


if </+»:/ <4 

it 4^ + Vt/ > 4 


( 2 . 2 ) 


The effects of tolerances on the link lengths and clearances in the hinges of a 
revolute dyad are shown in Fig. 2.4. The equivalent dyad which takes into account 



Chapter 2. Stochastic Modeling of Mechanical Error 



Figure 2.3: Stochastic Model of a Revolute Joint 


the stochastic effect is shown by the dashed hnes. The link lengths of the equivalent 
dyad are given by 

^ii = i^+ 

Or, 

Rii = Ri + 4 i = 1,2 (2.3) 

j = l,2 

So there are five random variables involved in the model of a revolute dyad. If 
these are denoted by I^i, V^, . . . , V5 then 


Vi = Ri i = l,2 

Vz = X12 
V4 = X'2Z 
Vs = X',, 


(2.4) 


As has been observed in practice the link lengths are assumed to be normally 
distributed. Therefore the mean values m[Vi] and the variances D[Vi] of the random 
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Figure 2.4: Stochastic Model of a Revolute Dyad 


variables are given as follows [28]: 


m[Vi] = 

k 

mi = 

( ^i^i\ 

\Y J 

mlVi] = 

mm 

D[Vi\ = 

^12 

4 

niui = 

^23 

4 

Dim = 

ill 

A 


i = 1,2 
i = l,2 
mlVs] = 0 

(2.5) 


2.2.2 Random Error in a Spatial Revolute Joint 

Figure 2.5a shows a revolute joint in a spatial mechanism. It is assumed that there is 
no wobbling of the pin in the socket. The socket and the pin are connected to hnks 1 
and 2 respectively. Oc is the point of intersection of the axis of link 1 and that of the 
socket. Consider now a plane r perpendicular to the socket and passing through Oc- 
This plane intersects the pin axis at Op. The distance OcOp is the clearance error in 
the revolute joint. 7 denotes the angle made by link 1 with the plane r. 

A coordinate system <S'(i,j,k) is set up with its origin at Oc (Fig. 2.5b). The 
unit vector k is directed along the axis of the socket. The vector OcOp is uniquely 
determined by two parameters r and a. These parameters are random in nature. 
Their ranges are given by 0 < r < c and 0 < a < 2%; where c is the radial clearance. 
The hmiting area for the random vector OcOp is a circle with its center at Oc and 



Chapter 2. Stochastic Modeling of Mechanical Error 


25 



Figure 2.5: Clearance Model for a Revolute Joint in a Spatial Mechanism 


radius c. Dhande and Chakraborty [30] have found the mean and the variance of the 
random error in the spatial joint as discussed in the following. 

As has been observed in practice, it is assumed that the probability that at any 
instant the axis of the pin lies inside or on the circumference of the circle is uniform. 
If /(r, a) is the joint probability density function and if dA is an elemental area at the 
tip of the error vector then the probabihty f(r, a) dr da is proportional to dA. Hence, 


/(r, a)drda = k^dA 


( 2 , 6 ) 


where is a proportionality constant. 

The constant k 2 and the joint probability density function f(r, a) is given by 

1 


k2 = 

/(r,Qi) = 


TTC^ 


TTC^ 


The effective clearance error e along the link axis is. given by 

rcosa 


€ = 


COS7 


The mean and the variance of the variable e are given as 

m[e] = 0 

■ ■ 

£>[£| = 


4008^7 


(2.7) 


( 2 . 8 ) 


(2.9) 

( 2 . 10 ) 
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For a planar mechanism the angle 7 = 0°. Therefore 

e = r cos a 

m[e] = 0 (2.11) 

CW = T 

So expressions (2.8), (2.9) and (2.10) reduce to those for the random variable and 

its mean and variance, respectively, for a revolute pair in a planar mechanism as in 
(2.5). 

2.2.3 Random Errors in the Coordinates of a Point 

Consider a point P(xi,yi, Zi) in the model of an RP process, the coordinates of which 
are in error and tolerances are specified on its coordinates. Tolerances are specified 
in absolute manner in contrast with the tolerance coefficients for link lengths. Then, 
these coordinates are random variables in the stochastic model of the RP process. Let 
Ci, Cj+i and 61+2 be fbe absolute tolerances on coordinates Xi, pi and Zi respectively. 
The corresponding random variables Vi, Vi-i-i and Vi+2 are then given by 

Vi = Xi + ei 

Vi^\ = yi + Ci+i (2-12) 

Vi +2 = -Zi + 61+2 

The coordinates of the point are in error due to locating the point in space while 
fabricating and assembling the RP machine. These are observed to be normally 
distributed in practice. Assuming normal distribution the means and the variances 
of the coordinates are given by 


m[VSl 

= Xi 


= Vi 

>n|V;+2] 

= Zi 

on 


0 (1^+, 1 


BW+jI 

= (¥)' 


2.2.4 Random Errors in the Orientation of a Link or an Axis 

Consider a link AB lying in 3D space (Fig. 2.6). Let XYZ be the base coordinate 
frame. 



Chapter 2. Stochastic Modeling of Mechanical Error 


27 



Let Yi be an axis along link AB with origin at A. Let ¥3 be an axis parallel to Yj, 
but passing through the origin 0 of the base coordinate frame. The orientation of a 
vector OY3 can be expressed by two successive rotations about two axes in space. So 
we can express the orientation of Y 3 by two successive rotations a and about base 
coordinate axes X and Y respectively. Since Yi is parallel to Y3, therefore, the two 
angles a and /? describe the orientation of Yi also. 

If the orientation of link AB is in error and ej and €^^.1 be the absolute tolerances 
on angles a and /3 respectively then the actual angles are given by 

Vi — Q! + Cj 

+ Cj+i 

Obviously Vi and are random variables. The angles giving orientation of a link 
are observed to be normally distributed in practice. Assuming normal distribution of 
the random variables Vi and Vi+i, their means and variances are given by 

m[Vi] =. a 

m[Vi+i] = /3 

Dm = (I)' (2-15) 

2.2.5 Stochastic Model of an RP Process 

The coordinates of a point on the work surface traced by the tip of the extruder head 
or the laser beam can be expressed as a function of random variables involved in the 


(2.14) 
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corresponding process. If there are n random variables Vi, V2, involved in the 

RP process under consideration then the coordinates of the point are given by 

X = x(Vi,V2,...,Vn) 

y = (2.16) 

^ = ^(^1,^2,. ...K) 

Since the dependent variables x, y and z are functions of random variables, there- 
fore, any given range of a dependent variable may be associated with the correspond- 
ing probability if the probability densities or at least certain numerical characteristics, 
such as means m[T^] and variances D[Vi] of random variables Vi are known. Means 
and variances of several random variables which may be involved in an RP process 
have been discussed earlier. 

To evaluate the numerical characteristics of x we proceed as follows. As the 
theorems of random variables are well estabhshed for linear functions, the right hand 
side of Eq. (2.16) is expanded in Taylor series about the means of the random variables 
and the terms having order two and above are neglected. As a result we get 

" ( dx\ 

X c:i x{m[Vi], i = l,2,...,n)-t-^^|^— j (P* - m[K]) (2.17) 

The term is a constant for a particular Vi- Let us say 

Therefore, x can be expressed as 

xc:^x(m[Vi], i = l,2,...,n)-t-X^Ai(K--w[l^]) (2.19) 

t=l 

In the present case Vi, i = l,2,...,n; are independent variables. Independent 
variables have zero covariance and, therefore, they are uncorrelated [69]. Hence the 
mean value of x is given by 

m[x] —x{m [Vi] , i = 1,2, . . ., n) (2.20) 

The expression (2.19) on substituting (2.20) into it reduces to 

x = m[x] + Y^Ai(yi-rn[Vi]) 

1=1 

X — m\x\ = Ai(Vi — m[Vi]) 

i=l 


( 2 . 21 ) 
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The expression for the variance D[x] of the dependent variable x can be derived 
as follows [56]. Variance D[x] is defined as [71] 

D[x] = E(x — m[x\)'^ 

where E(w) is the expected value of the dependent variable w. Erom Eq. (2.21) we 
get, 


D[x] = E 


Y,MVi-min) 


1 2 


Li=l 

== E\ «-mK|)(Vi-miH)) 

= EE 

i=i fc=i 


where 


and 

Cjk = Cov(V^-, Vfc) 1 < (j, k) < n 


If the n random variables Vj are uncorrelated, 


Cov(F^,Vfc) = 


Var(^-), 

0 , 


j=k 

jVk 


i.e., the covariance matrix [C] is diagonal. 

In the present case, all the n random variables are independent and hence uncor- 
related. Therefore, 

D[x\ = E{x — m[x])^ 

= Var(V5) 

i=l 

= EA^ mi 

i=l 


So expressions for numerical characteristics m[x] and D[x] of the dependent vari- 
able x are given by (2.20) and (2.22) respectively. 
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Similarly, numerical characteristics of the dependent variables y and z can be 
obtained. Taylor series expansions of y and z neglecting terms of order two and 
above are given by 


y ~ y{m[Vi], 
z cz. z{m[Vi\, 


” / B'>i\ 

i = 1.2 + (V;-m|l/,l) 

i=l 


(2.23) 


The numerical characteristics of y and z are given by 


Let 


QjI 


m[y\ 

m[z\ 

D[y] 

D{z] 


= 1 = 1,2, ...,n) 

= z{m[Vi\, i=l,2, ...,n) 






(2.24) 


( dxY _ ( dyV 



2 


m 


The terms ai, bi and Ci are called influence coefficients of x, y and z respectively. 

To associate an exact probability for any dependent variable x, y or z to lie within 
a certain range, it is necessary to determine distribution functions F(x), F(y) or F(z). 
Even if the probability density functions of all the random variables were known 
exactly, it is a difficult task to evaluate the distribution functions F(x), F(y) or F(z) 
[28]. 

The evaluation of F(x), F(y) and F(z) can be avoided if we use the central limit 
theorem regarding functions of random variables. This theorem states that iiVi, 

. . . , Vi are mutually independent random variables with finite mean and variance then 
the sum 


Si = J^Vi 

i 


tends to a normal variable if no single variable contributes to the sum as i tends to 
infinity. It is known that if the function is of the form given in (2.20) or (2.23) and if 
i > 5, then the dependent variable x, y and z may as well be taken as normal [123]. 

To evaluate the variances D[x], D[y] and D[z], the partial derivatives of x, y and z, 
respectively, with respect to the variables Vi are to be evaluated. Once the variance of 
a dependent variable has been found for the RP process to trace a particular point on 
the work surface, the range of dependent variable is evaluated for three-sigma band 
of confidence level (with probability 0.9973). 
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ANALYSIS OF MECHANICAL 
ERROR IN LOM 


3.1 Geometric Modeling for Description of the LOM 
Process 

The Laminated Object Manufacturing (LOM) process comprises of a sheet material 
moving between two idle rollers and a laser beam striking the sheet material as shown 
in Fig. 1.2. The sheet material is supplied from a continuous roll on one side of 
the machine and is taken up on the opposite side. An optical device moves in the 
horizontal plane with the help of an XY positioning device. The optical attachment 
emits the laser beam in vertically downward direction and the laser beam traces the 
contour of a slice and cuts the sheet [ 2 , 63, 31, 45]. The LOM process has been 
discussed in detail in Section 1 . 1 . 2 . 1 . 

The geometric model of the LOM process would allow to express a point on the 
contour of a slice. The geometric model is discussed below (Fig. 3.1). The rollers are 
modeled as two cylinders Ci and C 2 with centers of the front ends at Oi(ai, bi,ci) 
and 02 (^ 2 , 621 ^ 2 ), respectively, with respect to a base frame XYZ. The cylinders Ci 
and C 2 have radii ri and r 2 , respectively, and lengths li and I 2 respectively. The axes 
of the two cylinders are not parallel due to the errors in the location of centers Oi 
and O 2 and the errors in the orientation of the two axes. Ideally, the laser beam 
should be oriented in the vertically downward direction. However, due to mechanical 
errors in the location and orientation of the optical attachment, the actual orientation 
of the laser beam is oblique in space. The actual inclination of the laser beam can 
be described by means of two angles 6 and <j) as shown in Figs. 3.1 and 3.8a. This 
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general orientation of the laser beam would allow to find the partial derivatives of the 
coordinates of a point on the work surface with respect to 9 and <p. In computation, 
the nominal values of 9 and 0 will be kept such that the mean position of laser beam 
is vertical. 



Figure 3 . 1 : Geometric Model of the LOM Process 

The source of the laser beam EF is at E(a3, 63,03) and it strikes the work surface 
ABDC at F{x, y, z). Ideally, the axes of the two cylinders should be along Y direction, 
but due to mechanical errors they are inclined from the Y direction. A coordinate 
frame XiY\Z\ is attached at the center of the front end of cylinder Ci with Xi and Yi 
axes on the plane of the front end. The frame XiYiZi is obtained by two successive 
rotations and / 5 i about lines passing through Oi and parallel to Aand Z directions 
respectively. Another coordinate frame X2Y2Z2 is attached at the center of the front 
end of the cylinder C2 with X2 and Y2 axes on the plane of its front end. The frame 
X2Y2Z2 is obtained by two successive rotations 02 and (32 about fines passing through 
O2 and parallel to AT and Z directions respectively. 

Let pi(ti,u) and q2 (r,s) be the parametric representation of the cylindrical sur- 
faces Cl and C2 with respect to the coordinates frames XiYiZi and X2Y2Z2 respec- 
tively (Fig. 3 . 2 ). The parameters v and s are along the axes of Ci and C2 respectively. 
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Figure 3.2: The Cylindrical Surfaces pi('u,u) and q 2 (r, s) 


The surface pi(w, v) in the frame XiYiZi is given by 

pi(u, u) = ricosw risinu —liv ij 0 < u < 27r, 0<u<l (3.1) 

The surface q 2 (r, s) in the frame X 2 Y 2 Z 2 is given by 

q 2 (r, s) = r 2 Cosr r 2 sinr —I 2 S l] 0 < r < 27r, 0<s<l (3.2) 

The parametric representation qi(r, s) of the cylinder C 2 can be expressed in the 
frame XiYyZi as follows. Let [Tq] be the homogenous transformation matrix of the 
frame XiY^Zi with respect to the frame XYZ. Let [Tq] be the transformation of the 
frame X 2 Y 2 Z 2 with respect to the frame XYZ. Therefore, the frame X 2 Y 2 Z 2 with 
respect to XiYiZi, as shown in Fig. 3.3, is given by 

[r^] = [t^] [i?] 

[rf] = [r^] (3.3) 

Hence the surface qi(r, s) in the frame XiYiZi is given by 

qiCr.s) = q2(>',s) [t?] 


(3.4) 
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Figure 3.3: Transformations [Tq], [Tq^] and [Tf] 

3.1.1 The Transformation Matrices [Tq^], [Tq] and [Tf] 

The frame XiYiZi is obtained from the frame XYZ by a rotation about X by 90°, 
a rotation about X by Q!i, a rotation about Z by pi and a translation 
successively, as shown in Fig. 3.4. The frame XjYjZi from XY Z is obtained by a 90° 
rotation about X. The frame XjiYuZu is obtained from XjYjZj after two successive 
rotations, ai and Pi about X and Z directions respectively. The transformation [Tq] 
is given by 

[^ToJ = Rot{X, 90°) Rot{X, ai) Rot{Z, pi) ^ Trans{Base,ai,bi, Cl) 

This yields [103] 

[To^j = Rot(X, 90°) Rot{X, cki) Rot{Z, Pi) Trans{Base, ai, 6 i, ci) 

Or, 

cos Pi sin Pi 0 0 

sincKisin/?! — sin ai cos /?i coso;i 0 
cos Oil sin /?! — coscKiCOS/di — sinQ;i 0 

fli bi Cl 1 

Similarly the matrix [Tq] can be obtained as 

[Tq^] = Rot(X, 90°) Rot{X, 0 : 2 ) Rot{Z, P 2 ) Trans{Base, 02 , 62 , C 2 ) 
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Figure 3.4: The Coordinate Frames XYZ and XiYiZi 


Or, 



cos jd 2 sin ^2 0 0 

sinQ!2sin/32 — sina2Cos/32 cos 0:2 0 
cos a 2 sin P 2 — cos 0:2 cos (32 — sin 0:2 0 
0,2 ^2 C 2 1 


Substituting [Tq] and [Tq] in Eq. (3.3) it follows that 


(3.6) 


(I?) 


cos /?2 sin/52 0 0 

sin a 2 sin P 2 — sin 0:2 cos cos ot 2 0 
cosa2sin/32 —cos 0 : 2 cos /^ —sin 0:2 0 
02 h ^2 1 _ 

cos /5i sin a!i sin /5i cos ai sin /3i 0 

sin pi — sin ai cos /5i — cos 0:1 cos /3i 0 

0 cos Oil — sin Oil 0 

— tti cos Pi — bi sin Pi {-Oi sin ai sin Pi {-Oi cos ai sin Pi 1 

+bi sin Q!i cos Pi +bi cos ai cos Pi 

—Cl cos Oil) CisincKi) 


(3.7) 
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3.1.2 Mathematical Representation of the Work Surface 

The LOM sheet from which the contour of the slice is cut is the work surface. This 
surface is to be modeled mathematically so as to be able to find a point on the contour 
of a slice. A point on the contour is obtained by the intersection of the laser beam with 
the work surface. Consider the following three alternatives to approximate the work 
surface. In the first alternative the line AC, tangent to the two cylinders, is found 
such that the end points A and C lie on the front ends of Ci and C 2 respectively. 
This gives two points A and C on the work surface. Then two straight lines AB 
and CD are drawn passing through A and C, respectively, and parallel to the axes of 
the corresponding cylinders. The work surface can be modeled by a bifinear surface 
ABDC with the four corner points A, B, C and D as shown in Fig. 3.5. 


Bilinear Surface 



Figure 3.5: Mathematical Model of the Work Surface 

The second alternative is to take the two tangents at both front and rear ends 
of the cylinders. But then AB and CD will, be curved lines and wrapping has to 
be considered to find them. For an approximation, the bilinear surface in the first 
alternative is good enough as the work smrface is not highly curved, and it is a plane 
when there is no error. The third alternative is to take the highest points on both 
ends of the two cylinders giving four corner points. A bifinear surface can be obtained 
with these four points. Obviously, the first alternative is a better representation of a 
si^face tangent to the cylinders as the approximate surface in that case consists of at 



Chapter 3. Analysis of Mechanical Error in LOM 


37 


least one tangent line. 

So the work surface is modeled as in the first alternative. The line AC, tangent 
to the two cylinders is found such that its end points lie on the front ends of the 
cylinders. To approximate the surface, two lines AB and CD are drawn parallel to 
the axes of the corresponding cylinders and passing through the points A and C 
respectively. Let the points A, B, C and D correspond to four corner points of a unit 
square in parametric space, p(0,0), p(0,l), p(l,0) and p(l,l) respectively. A bilinear 
surface is constructed from these four points as corner points (Fig. 3.6). 



Figure 3.6: The Bilinear Surface ABDC in XYZ Object Space 


The coordinates of these four corner points are found in XYZ object space by the 
desired transformation. The bilinear surface ABDC in the frame XYZ is given by 




p(0.0) p(0,l) 
P(liO) p(l.l) 



(3.8) 


where the four comer points p(0,0), p(0,l)) p(liO) and p(l,l) are expressed in the 
object space. 

The surface is, in fact, a hyperbolic paraboloid, and therefore, it is a doubly ruled 
surface. 


3.1.3 The End Points of the Edge AC of the Work Surface 

The line AC in Fig. 3.5 is tangent to the two cylinders such that its end points A 
and C he on the front ends of C\ and C 2 respectively. Therefore, the normals to the 
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surfaces pi(w,u) and q 2 (r,s) are orthogonal to the line AC at A and C respectively. 
The other two conditions are that the end points A and C of the line should lie on the 
front ends of the cylinders Ci and C 2 respectively. Let A(a;i,yi,0) and C{x 2 ^y 2 ^Z 2 ) 
be the two points in the frame XiYiZi, as shown in Fig. 3.7. 


Figure 3.7: Tangent Points A and C in the frame X-CY^Zx 

The parametric representation of the line AC in the frame XiYxZi is given by [38, 
pp. 55-56] 

l(t) = [ Pi 0 1 ] -I- 1 [ 0:2 - xi y^ — Pi Z2 0 0 < t < 1 (3.9) 

The tangent to the line l(t) is obtained by differentiating Eq. (3.9) with respect to t: 

l*{t) =^X2-Xi y2-y\ Z2 0 (3.10) 

The surface pi(n, n) in the frame XiYiZi is given by Eq. (3.1), i.e., 

Pi(w, u) = ri cosw risiuu —Ixv 1 

The surface qi (r, s) in the same frame XiYxZx is given by 

qi (r,s) = r 2 COsr r 2 sinr —I 2 S 1 T? (3-11) 

There are four conditions (a), (b), (c) and (d) to find the tangent points A and C, 
as discussed below. The values of t on the line l(t) at the points A and C are 0 and 
1 respectively. The values of v and s at the points A and C, respectively, are both 0. 
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(a) The normal to pi(w, n) is orthogonal to the hne l(t) at A(t=0, v=0), i.e., 


(Pixp^).l‘ 


t=0,'u=0 


0 


where the parametric derivatives p“ and p| are the tangent vectors to pi {u, v) [38, 
90]. This yields 

( 2/2 - yi)sin« + (xa - Xi) cosw = 0 (3-12) 

(b) The normal to qi(r, s) is orthogonal to the line l(t) at C(t=l, s=0), i.e., 

where the parametric derivatives and qf are the tangent vectors to qi (r,s). 
This yields 


— ail sinr + a 2 i cosr — ai 2 sinr + a 22 Cosr — ai 3 sinr + a 23 cosr 
—CL^ih ~o^2h ~ciszh 

X 2 -X 1 y2-y\ Z2 


= 0 
(3.13) 


where Oy is the element in the row and column of the matrix [Tfj. 

(c) The line l(t) touches pi(u,u) at A(t=0, v=0), i.e., 

Pi(«,0) -1(0) =0 

This yields 


Xx = ri cos u 
yx = f'x sin u 

(d) The line l(t) touches qi (r,s) at C(t=l, s=0), i.e., 

qi(r,0)-l(l)=0 

This yields 


(3.14) 

(3.15) 


X2 

= ’’2(ctii COS r + a2i sin r) + a 4 i 

(3.16) 

yt 

= f 2 (ai2 cos r + a22 sin r) + 042 

(3.17) 

Z2 

= r 2 (ai 3 cosr + otz sin r) + a 43 

(3.18) 


Substituting a?!, t/i, 0:2, ^2 and ^2 from Eqs. (3.14), (3.15), (3.16), (3.17) and 
(3.18) respectively, into Eqs. (3.12) and (3.13) gives two transcendental equations in 
two unknowns u and r. 
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3.1.4 Solving for u and r by Newton-Raphson Method 

The two equations, Eqs. (3.12) and (3.13) are to be solved for the two variables u 
and r. The two equations, being transcendental, have to be solved by an iterative 
technique. Newton-Raphson method [97] is adopted to solve the two equations as 
follows: 

Let the left hand sides of Eqs. (3.12) and (3.13) be expressed as Fi and F 2 re- 
spectively. Therefore, the two transcendental equations in two unknowns u and r are 
given by 


Fi = (s /2 - Vi) sin w + ( 0:2 - xi) cos u 
F 2 = 


where Eq. (3.13) is written as 


/ 

9 

h 



p 

9 

r 


=z 

u 

V 

w 



as 







f 

9 

h 



P 

9 

r 



u 

V 

w 


(3.19) 

(3.20) 


= 0 


The terms xi, pi, X 2 , y 2 and Z 2 in the above are given by Eqs. 

El 


The system of equations is given by [F] 
of [F] is given by 

[J] = 


F 2 


(3.14) to (3.18). 

" 0 ' 


0 


where 


dFi 

du 


dFx 

dr 

dF2 

du 


+ 


(s/ 2 -S/i)cosu + 

dx2 dxi 

du du 

\ > 

dy2 _ 

dr dr J 

f 9 h 

p q r 
u V w 


■ dFi 

dFi 

du 

dr 

dFi 

dFi 

du 

dr 

dpi _ 

dyi'^ 

du 

du ) 


. The Jacobian 


(3.21) 


— — smu — ( 0:2 — xij sinu 


cosu 


sinw + 


'dx2 _ dxi 
dr dr 


cosu 


(3.22) 

(3.23) 


du 


The partial derivatives of the determinant F 2 can be expressed as the sum of the 
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three determinants [71, pp. 243] as follows 

m 

du 


f 

o' 

h' 


/ 

9 

h 


f 

9 

h 

p 

q 

r 

+ 

p' 

4 

r' 

+ 

P 

q 

r 

u 

V 

w 


u 

V 

w 


u’ 

v' 



(3.24) 


where prime denotes partial derivatives with respect to u. Since the first and second 
rows of F 2 do not depend upon u, therefore, 

f 9 h 


m 

du 


p q r 
u' v' w' 

The second row of does not depend upon r too. Therefore, 

m 

dr 


(3.25) 


r 

9' 



/ 

9 

h 

p 

q 

r 

+ 

p 

q 

T 

u 

V 

w 


u' 

v' 



(3.26) 


where the prime denotes partial derivatives with respect to r. 

So the Jacobian of [F] can be expressed analytically. 

There are four tangents to the cylinders at the front end. If our initial guess for 
both u and r are close to ^ 72 , then the desired result, i.e., the one for which both 
Hi and y 2 are positive, is obtained (Fig! 3.7). Newton-Raphson method, on solving 
the two transcendental equations, gives the values of the parameters u and r for the 
tangent AC. Let these values of u and r at the points A and C, respectively, be given 
by 

U = Ua 
r = rc 


(3.27) 


3.1.5 The Four Corner Points of the Work Surface in the 
Base Frame XYZ 


The four corner points of the work surface ABDC are described in Section 3.1.2. 
Once the tangent line AC is found, the corner points p(0,0), p(0,l), p(l,0) and 
p(l,l) correspond to {ua,v — 0), {ua,v = 1), {rc,s = 0) and (rc,s = 1) respectively. 
Hence the corner points p(0,0) and p(0,l) in the frame XYZ are given by 


p(0,0) 


tiCosua risiuuA 0 1 



. p(0) 1) . 


TxOOSUa risinw^ -l-y 1 


0.0 


The other two corner points p(l,0) and p(l,l) in the frame XYZ are given by 


(3.28) 


P(l)O) 


r 2 cos rc ra sin rc 0 1 

tI 



ra cos rc rasinrc —h 1 



And the bilinear surface is given by Eq. (3.8). 


(3.29) 
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3.1.6 Parametric Representation of the Laser Beam 

The optical attachment emits the laser beam to draw a contour on the work surface 
and cut the sheet. To find the intersection of the beam and the work surface, it is 
necessary to represent the incident beam mathematically. 

Let EF be the incident laser beam emitted from the optical attachment at E(a3, 63 , C3) 
and intersecting the XY plane at G as shown in Fig. 3.8. 



Figure 3.8: The Incident Laser Beam 


The parametric equation of the line EG is given by 


QW 


ag 63 C3 1 


+ t 


C3 tan 6 cos (j) C3 tan ^ sin (j) — C3 0 


The above equation holds for a vertical laser beam too, by substituting 9 
any value of 0, say 0 = 0°. 


(3.30) 
0° and 


3.1.7 Intersection of the Laser Beam and the Work Surface 

The laser beam EG in the base frame XYZ is given by 


Q(t) = 1 ^ as 63 C3 1 +t C3tan(9cos0 C3tan6'sin0 -C3 0 


The surface ABDC in the frame XYZ is given by 


q(u,u;) 


1 — w u 


p(0,0) p(0,l) 

P(l7 0) P(l)l) 



1 — w 


w 


as in Eq. (3.8). The defining points p(0,0), p(0,l), p(l,0) and p(l,l) are given by 
Eqs. (3.28) and (3.29). 

Hence the point of intersection of Q(t) and q(u,w) is given by [38, 90, 133] 


q(w,w) = Q(t) 


0 <t <00 


(3.31) 
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A frame X*Y*Z* is obtained by a transformation with respect to the frame XYZ 
such that the ^*-axis coincides with GE. See Fig. 3.9. The transformation of the 
base frame XYZ with respect to the X*Y*Z* is given by 

= Rot(Y\e) ^ Rot{Z2,-4>) -^Trans{-OG) 



Figure 3.9: Intersection of the Laser Beam and the Work Surface 
Let the above transformation matrix be simply represented by [T], 

Hence 

[T] = rrans(-OG) -4>) Rot{Y*, B) 

The translation OG is given by 


Let 

where 


OG = as + C3 tan 9 cos 0 63 + C3 tan 0 sin ^ 0 0 I 

OG = [ 5'a: Qy 0 0 
Pa; = 03 + C 3 tan 6 cos (j> 


(3.32) 


and 


Pj/ = &3 + C 3 tan 6 sin <() 
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Therefore, 

Or, 


[T] = Trans{-g^, -Qy, 0) Rot{Z 2 , -<p) Rot{Y*, 9) 


(3.33) 


IT] 


cos 6 cos 4) 

— sincf) 

— sin 6 cos (p 

0 ' 

cos 6 sin (f) 

Q,OS(p 

— sin 6 sin (p 

0 

sin 6 

0 

oos6 

0 

— cos 0 (a 3 cos (f> + hz sin (p 
+cz tan d) 

sin (p ’-bs cos cp 

sin 0 (as cos (p + bz sin (p 
+cz tan 6) 

1 


(3.34) 


The transformed surface is given by 


Let 


where 


cC{u,w) = q(u,-u;)[T] 

X y z I [T] 

q*(M,u;)=[ x*(u,w) y*(u,w) z*{u,w) 1 


(3.35) 


x*{u,w) y*{u,w) z*(u,w) 1 ] = [ x(u,w) y(u,w) z(u,w) 1 [T] (3.36) 


and 


p(0,0) p(0,l) 
p(l,0) p(l,l) 



1 — w 


w 

3 


.37) 


x(u,w) y(u,w) z(u,w) 1 =[l — u uj 

The transformed line is given by 

Q*(<) = Q(*) m 

The line EG is coincident with the Z*-axis. Hence the x and y coordinates of the 
line EG in the frame X*Y*Z* are both zeroes. Therefore, Q*(^) is given by 


Q*(t) = [0 0 z{t)\Q. 1 


(3.38) 


Equating right hand sides of Eqs. (3.35) and (3.38) gives the point of intersection 


as 


x*{u,w) = 0 
y*{u,w) = 0 


(3.39) 

(3.40) 


The defining points p(0,0), p(0,l), p(l,0) and p(l,l), and the transformation 
matrix [T] do not depend upon the parameters u and w. It can be seen that the 
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expressions for x* and y* are non-linear in u and w. Eqs. (3.42) and (3.43) can be 
solved for u and w by Newton-Raphson method, an iterative technique [97]. 

For the intersection 


Fi — x*(u,w) ~ 0 
F 2 = y*{u,w) = 0 


The Jacobian [J] of [F] = 


Fi 

F 2 


is given by 


r dF, 


{J] = 


du 

m 

du 


dFi - 
dw 

m 

dw . 


(3.41) 

(3.42) 


(3.43) 


The Jacobian [J] can be expressed analytically as the expressions for Fi and F 2 involve 
only multiphcation terms of u and w. Since the elements of [T] do not depend upon 
u and w, therefore, the partial derivatives of F\ and F 2 with respect to u are given by 


■ dFi 

dF2 

dz* ^ ' 

_a_ 



. du 

du 

du ^ J 

du 

X* y* z* 



du 


X y z 1 I [Tj 

[-1 1] 


p(0,0) p(0,l) 

p(l,0) p(l, 1) 



1 — w 

J 

w 

L J 


[T] 


and the partial derivatives of Fi and F 2 with respect to w are given by 


dFi dF 2 dz* 


dw dw dw 


dw 

dw 


[ X* y* z* 1 
[ X y z 1 ] [T] 


= ’ 1 


u u 


p(O.O) p(O.l) 
p(l>0) p(l,l) 



’ -1 ' 


1 

• 


m 


(3.44) 


(3.45) 


Hence the elements of Jacobian [J] are given by Eqs. (3.44) and (3.45). Solving 
Eqs. (3.41) and (3.42) by Newton-Raphson method gives u and w. Substituting u 
and w in Eq. (3.37) gives x, y and z. Hence the point of intersection P is given by 

q = [ X y 2 ; 1 

where q is the position vector of a point on the contour of a slice cut by the laser 
beam. 
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3.2 Stochastic Error Analysis of the LOM Process 


Chapter 2 discusses stochastic modeling of an RP process. Following the approach 
given there, the stochastic error analysis of the LOM process can be carried out as 
follows. 

The dimensions in error in the LOM process are shown in Fig. 3 . 10 . The dimen- 
sions in error are the coordinates ai, cj, a2, C2, as, bs and cg; the angles ai, 

Pi} ^2} ^2} ^ and <j), and the link lengths L, r2 and Is- The error in aj is Aai, the 
error in 6i is A6i, and so on. There are 19 dimensions in error. Let Vi, i=l to 19 , be 
the random variables involved in the process. Using the approach given in Chapter 2, 
the random variables involved in the LOM process are as follows: 


Vx = Ui + €i , 

Vs —Pi +^5} 
Vs = <22 T £$) 
Vi 2 = /?2 + ei2, 

Vl5 = as + 6X5, 

Vig = <P + ei9. 


V2 = ^>1 + 62) 

Ve = ri(l + ee)) 

Vg = 62 + 69 ) 

Vi 3 = rsil + 613), 

Vl6 — bs + 6x6, 


V3 = Cl + 6 s, 

V^ = lx(l H- 67), 

Vio = C2 + 6x0, 

Vx 4 = 12(1 + 614), 

Vl7 = C3 -I- 617, 


V4 — <3!l -F 64, 

l"u = a2 + £n. 

V18 = 6 + 618 , 


where ei, 62, 63, 63, 69, eio, 615, 6x6 and 617 are absolute tolerances in meters (m) on 
coordinates; 64, 65, en, 612, ei8 and 619 are absolute tolerances in radians on angles; 
and €6, 67, 613 and 614 are tolerances per unit nominal lengths on Vx, lx, f'2 and I2 
respectively. The random variables kj are independent of each other. 

The mean values of the random variables m[K]) i=l to 19 , are ax, bx, Cx, ax, Px, 
O) k, 0^2) b2, C2, a2, P2, T2, I2, as, bs, C3, 0 and (j) respectively. 

The variances D[I^] are given as follows: 

D[VA = (f i = 1 to 5 , 8 to 12 , 15 to 19 

Dm = 

Dm = (3.47) 

D[Vn] = 

Divu] = 

The coordinates of a point on the work smrface traced by the laser beam are given 
by 


X = xiVx,V 2 ,...,Vxg) 
y = y{Vx,V 2 ,...,Vxg) 
z = ziVx,V 2 ,...,Vxg) 


( 3 . 48 ) 
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Figure 3.10: Geometric Model of the LOM Process Showing Errors in Dimensions 

The means and the variances of the dependent variables x, y and z are given as 
follows: 

m[x] = x{m [Vi] , i = 1,2, . . ., 19) 

m{y\ = y(m[l^], i = 1,2, . . ., 19) 

m[z\ = 2: (m [Vi] , i = 1,2, . . ., 19) 

D[A om (3.49) 

i=l \ ^/ m 

™ - S (S)’ 

™ - £ ( 5 )’ *' 

The partial derivatives are evaluated at the mean values of the variables Vi as 
denoted by the subscript m. The mean values in[Vi] of the variables Vi are as given 
earlier in this subsection. The necessary expressions for the influence coefficients 
{MT’ iwY derived as in the following. 

The coordinates of a point on the work surface is given by Eq. (3.8): 


- 

r . 1 

1 

p(0. 1) 


1 — w 

y z 1 

1 

r-H 

II 

p(l,0) 

p(l,l) 


w 
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Therefore, 

X = g[l] = pii{l — u){l — w) - u)w + pziu{l — w) piiuw (3.50) 

y = <?[2] = Pn{'^ - u){l - w) + p 22 {l - u)w Jr pz 2 u{l - w) + p^ 2 '<J-'^ (3.51) 

2 : = q[3] = fii 3 (l - tt)(l - w) + p 2 z{l - u)w + P 33 «(l — ui)+ pizuw (3.52) 

where pij is given by 


[P] = 


(3.53) 

The expressions for p(0,0), p(0,l), p(l,0) and p(l,l) in the above equation are 
given in Eqs. (3.28) and (3.29). The coordinates Pij in the above equation are with 
respect to the base frame XYZ. Therefore, the points on the surface q are also given 
with respect to the frame XYZ. 

Equations (3.5) and (3.6) give 


pll 

P 12 

PlZ 

1 ' 


'p(o,o) ■ 



ri cos u 

Ti sin u 

0 

1 ’ 

|r„‘] 

P 21 

P 22 

P2Z 

1 


p(0,l) 



ri cos u 

Ti sin u 


1 

Pzi 

PZ2 

Pzz 

1 


P(l.O) 



T2 cos r 

r2 sin r 

0 

1 ’ 

[T}] 

_ P41 

P42 

P4Z 

1 


p(l.l) 


1 

. i 

r2 cos r 

r2 sin r 

-h 

1 




= To^(ai,6i,Ci,Q!i,A) 

= Tq {a2,b2,C2,a2,p2) 


Therefore, 


Pii and p 2 i fn(ai, 61 , Cl, Oil, /3i, ri, k, u) i = 1 to 4 

= 1 to 4 


(3.54) 

(3.55) 

(3.56) 

(3.57) 


Psi and p 4 i -)■ fn(a 2 , 62 , C 2 , Oi 2 , 02 , ’' 2 , ^ 2 , r) 

The parameters u and r in the expressions for pij are the solution of Eqs. (3.12) 
and (3.13). Rewriting Eqs. (3.12) and (3.13) 

( 2/2 — Pi) sinu + (^2 — ®i) cos rx = 0 


—an sinr + a 2 i cosr 
—cLzih 

X2 — Xi 


—ai 2 sinr + 022 cos r 

—Cl32h 

P2 - Pi 


— ai 3 sin r + a 23 cos r 

— 033^2 

Z2 


= 0 


where 


X\ = Ti cos u 

= Tisinu 

X 2 — r 2 (an cos r + a 2 i sin r) + O 41 

y 2 = r 2 (ai 2 cos r + a 22 sin r) + a 42 

Z 2 = r 2 (ai 3 cos r + an sin r) + a 43 
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Therefore, 

X = g[l] = pii{l — u){l — w) + p 2 \{l — u)w + pziu{l — w) + piiuw (3.50) 

y = g[2] = pi 2 (l -«)(!- to) + P22(1 - u)w + P 32 w(l — + Pa^uw (3.51) 

^ = 9[3] = Pi3(l — t«)(l - w) + P23(1 - u)w + P33w( 1 — ti^) + Pazww (3.52) 

where pij is given by 


b] 


(3.53) 

The expressions for p(0,0), p(0,l), p(l,0) and p(l,l) in the above equation are 
given in Eqs. (3.28) and (3.29). The coordinates pij in the above equation are with 
respect to the base frame XYZ. Therefore, the points on the surface q are also given 
with respect to the frame XYZ. 

Equations (3.5) and (3.6) give 


Pn 

P 12 

Pl3 

1 ' 


■ p(0,0) ■ 



Ti COS U 

Ti sin u 

0 

1 ' 

Pol 

P21 

P 22 

P 23 

1 


p(0, 1) 



Ti COS U 

ri sin u 

■“^1 

1 

Psi 

P32 

Pzz 

1 


p(l,0) 



r 2 COS r 

r 2 sinr 

0 

1 ’ 

[I?] 

_ P41 

P42 

P4Z 

1 





r 2 COS r 

r2 sin r 

—k 

1 


■‘■0 

^2 

0 


ri 


— T(j (Qi\ , 5l , Cj , OC-y , (3y ) 
= Tq (02, 62 ) C2, 0:2) Pz) 


(3.54) 

(3.55) 


Therefore, 


1 to 4 
1 to 4 


(3.56) 

(3.57) 


piiandp2i fii(oi,6i,ci,Q!i,/?i,ri,ii,w) 

pzi and pAi -> fn(a 2 , 62 , C 2 , 0 : 2 , P 2 , ^ 2 ) k,r) 

The parameters u and r in the expressions for Pij are the solution of Eqs. (3.12) 
and (3.13). Rewriting Eqs. (3.12) and (3.13) 

( 2/2 “ 2 / 1 ) sinu + (x 2 — a:i) cos w = 0 


where 


= 0 


—ail sinr + 

021 cos r 

—012 sinr + 022 

cos 

r —013 sin 

— 031/2 


—0,32k 


—033/2 

a;2 — xy 


2/2 -2/1 


Z2 


aji 

= ri cos u 




yi 

= risinw 




X2 

= ^2(^11 cos r + 

021 

sin r) + 041 


2/2 

= r2(oi2cosr + 

022 

sin r) + 042 


Z 2 

= r2(oi3Cosr + 

^23 

sinr) + 043 


48 
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In the above expressions {xi^yi) and {x 2 ^y 2 iZ^ are coordinates with respect to 
the frame XxYxZi. 

Let Eq. (3.13) be represented by 

/ll /l2 /l3 

/21 /22 /23 = 0 (3.58) 

fzx hz fzz 

Differentiating Eqs. (3.12) and (3.58) with respect to the random variables Vi gives 

5 ( 2/2 — 2 / 1 ) . , / ^ du d(x 2 -xx) 

4 smw + {y 2 — yx) cosw-;;:^^^ H irrTT : — - cosu 


dVi 


dVi 


dVi 


— {x 2 — Xx) sintt 


du 

d\^ 


0 


fix f'x 2 fxz 


/ll /l 2 fxz 


/ll fx 2 fxz 

/21 /22 fzz 

+ 

fix fiz fiz 

+ 

fzx fzz fzz 

fzx fz2 fzz 


fzx fz 2 fzz 


fzx fzz fzz 


(3.59) 


(3.60) 


where prime denotes derivative with respect to V 

Equations (3.59) and (3.60) involve derivatives of the elements of the matrix j, 
given by Eq. (3.7). Let in Eq. (3.7) be represented by 

®11 ^12 ^13 0 

Ci2X <222 <223 0 

OiZl 0'Z2 <233 0 
<241 ®42 <243 1 

Equation (3.7) gives \Tl] in terms of the mean values of the variables V^. In terms of 
the variables Vi, Uij is given by 


[r’] = 


(3,61) 


a 




= aij{Vx,V2,...,V^,V^,V,...,Vx2) 


(3.62) 


In addition to Fi, V 2 , . . . , I 5 and 1^, Lg, . . . , V 12 the coordinates of Xx,yx,X 2 , y 2 and 
Z 2 are functions of Vg = ri and V 13 = r 2 (and, of course, u and r). Hence Eqs. (3.12) 
and (3.13) which are two equations in two variables u and r, involve Vx to Vi 3 except 
Vj, i.e., 

u = u{Vx,V2,...,V^,Vi,V,,...,Vxz) (3.63) 

r = r{Vx,V2,...,V6,Vs,Vg,...,Vx3) (3.64) 

Also 

^i,yi,X2,y2,Z2 i'xx{Vx,V2,....,V^,V^,VQ,...,Vxz) (3.65) 

Partial derivatives of the elements of [Tf] and ®i, j/i, ^ 2 , yz and Z 2 with respect to 
the variables Vi are derived in the following subsections. These derivatives are needed 
in solving Eqs. (3.59) and (3.60) for ^ and ^ for all Ws. The derivatives || and 
^ are needed for the influence coefEcients ^ and 
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3.2.1 Partial Derivatives of the Elements of [Tf] 

Equation (3.7) gives [T^^] in terms of the mean values of the random variables. The 
elements of [Tf] are given by Eq. (3.61). The expressions for each a^- is written 
comparing Eqs. (3.7) and (3.61). The partial derivatives of a^- is derived at the mean 
values of random variables in the following. 

an = cos(l^ - V 12 ) 

Hence, 

= - sin (A - /? 2 ) 

= sin(,0i-/?2) (3.66) 

= 0 i = 1 to 4, 6 to 11 , 13 to 19 


^<^11 

dV^ 

daxx 

dVn 

dan 

dVi 


ai 2 = sin V 4 sin(V 5 — Vn) 


Hence, 


Hence, 


dai2 


dV4 

dai2 


m 


dVs 

dai2 

dVn 

dai2 

dVi 


\m 


m 


m 


= cos O'! sin (jdi — P 2 ) 

= sin Q!i cos (^1 — P 2 ) 

= — sin Oil cos {Pi — P 2 ) 

= 0 i = 1 to 3 , 6 to 11, 13 to 19 


«13 = cos V 4 sin(l^ — V 12 ) 


(3.67) 


daiz 


dV4 

daiz 

dVs 

daiz 

dVi2 

daiz 


m 


m 


m 


dVi 


im 


— sin Oil sin {Pi — P2) 
cos Q!i cos {Pi — P2) 

— cos ai cos {Pi — P 2 ) 

0 i = 1 to 3 , 6 to 11, 13 to 19 


(3.68) 



a2i = sin Vii sin(Vi2 - V5) 


Hence, 


da2i 


dVs 

dci2i 


m 


dVi^ 

da2i 

'm2 

da2i 


m 


m 


dVi 


m 


= - sin 0:2 cos (/?2 - Pi ) 

= cos 0:2 sin (/?2 —A) 

= sin Q!2 cos (p2 — Pi) 

= 0 i = 1 to 4 , 6 to 10 , 13 to 19 


( 3 . 69 ) 


a22 = sin V4 sin V12 cos(V^ — V12) + cos V4 cos Vii 


Hence, 


da22 

'm 

da,22 

'm 

8022 

8022 


m 


Im 


m 


8V12 

80,22 


8Vi 


— cos ai sin 0:2 cos (A — P2) ~ sin ai cos 0:2 
= — sin Oil sin 012 sin (A — A) 

= sin Q!i cos Q!2 cos (A — A) — cos Oil sin 0:2 
= sin CKi sin q; 2 sin (A — A) 

= 0 i = 1 to 3 , 6 to 10 , 13 to 19 


( 3 . 70 ) 


023 — sin Vii cos V4 cos(V5 — V12) — sin V4 cos Vn 


Hence, 


8023 

'm 

8023 

'm 

8023 

Wi 

8023 

8023 


Im 


im 


m 


Im 


8Vi 


m 


— sin Q!i sin 0:2 cos (A — A) “ cos oii cos Q!2 

— cos ai sin a2 sin (A — A) 

cos ai cos 0:2 cos {Pi — A) + sin ai sin Q!2 
cos ai sin 0:2 sin (A “ A) 

0 i =; 1 to 3 , 6 to 10 , 13 , to 19 


sfTOra* 


»rr7cfhr 
^rfc<r spo A 



( 3 . 71 ) 
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Hence, 


Hence, 


Hence, 
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asi = cos Vii sin(yi2 - V 5 ) 


dasi 


dVs 

dasi 

dVn 

dazi 

dViz 

dazx 


m 


m 


\m 


dVi 


m 


= - cos Q!2 cos {/32 - A ) 

= — sin 0:2 sin (A — A ) 

= cos q' 2 cos (A — A) 

= 0 i = 1 to 4, 6 to 10, 13 to 19 


(3.72) 


agj = sin V 4 cos Vn cos(V5 — V 12 ) — cos V 4 cos Vu 


daz 2 

daz 2 

dVz 

daz 2 


m 


m 


dVxx 

daz2 

W 2 

daz 2 

dVi 


\m 


m 


\m 


cos a\ cos Q!2 cos (/3i — A) + sin ax sin az 

— sin ax cos az sin (A — A) 

— sin ax sin 0:2 cos (A ~ A) “ cos ax cos 0:2 
sin ax cos 0:2 sin (A ~ A) 

0 i = 1 to 3, 6 to 10, 13 to 19 


(3.73) 


azz = cos V4 cos Vxx cos(p5 — H12) + sin V 4 sin 


dazz 

dV 4 

dazz 


m 


dVz 

dazz 

dVxx 

dazz 

% 

dozz 


\m 


\m 


m 


dVi 


\m 


— sin ax cos 0!2 cos (A ~ A) + cos ax sin 0:2 

— cos ax cos Q'2 sin (A “ A) 

— cos Q!i sin 0:2 cos (A ~ A) + cki cos q ;2 
cos ax cos Q '2 sin ( A “ A) 

0 i = 1 to 3, 6 to 10, 13 to 19 


(3.74) 


52 
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a4i = {Vs - V,) cos H + (^9 - Vi) sin 


Hence, 


da4i 


dVi 

da^x 

ddix 

dCL4,x 

W 

ddix 


\m 


\m 


\m . 


Im 


dV^ 

dan 

~dV 


Im 


Im 


— cos Px 


— sin Px 

- (02 - ax) sinPx + {h - h) cosPx 
cosPx 

sinPx 

0 i = 3 , 4 , 6, 7 , 10 to 19 


(3.75) 


an = {Vs — Vx) sin V4 sin V5 — (Vg — V2) sin V4 cos V5 + (Vio — Vz) cos V4 


Hence, 


5042 


dVx 

5042 


dVi 

5042 


dVs 

5o42 


dV4 

5o42 


dVs 

5o42 


5^8 

5o42 


dVg 

5o42 


51^10 

5o42 


dVi 


= — sin ax sin Px 
= sin CKi cos Px 
= — cos Q!i 

= (02 — Oi) cos ax sinPx — {^2 ~ ^1) cos ax cosPx — (c2 — Ci) sin ax 
= (02 — Oi) sin ax cos Px + (62 — ^*1) sin ai sin px 
= sin Q!i sin /?i 
= — sinQ;i cos/ 3 i 

= cosai 


m 


0 


i = 6 , 7 , 11 to 19 


( 3 . 76 ) 
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043 = {Vs - Vl) cos V 4 sin Vs - (Vg - V 2 ) cos 1^4 cos ^5 - (i/10 - Vs) sin V 4 


Hence, 


da^s 

dVi 

da4s 

dVs 

da4s 

dVs 

da43 

dV4 

da43 

dVs 

da4s 

dVs 

^043 

dVg 

(9043- 

9^10 

da43 


m 


m 


m 


m 


m 


m 


m 


m 


dVi 


\m 


— cos (Xi sin Pi 
cos O'! cos Pi 
sin O'! 

— {as — Oi) sinai sin/3i + (62 — 5i) sinai cos/3i — (cg — Ci) coscti 
{as - ai) cos 0:1 cos Pi + (62 - h) cos ai sin pi 

cos Oil sin Pi 

— cos CKi cos Pi 

— sin ai 

0 i = 6, 7, llto l9 

(3.77) 


aj 4 = constant j = 1 to 4 


Hence, 


daj4 



j = 1 to 4, i = 1 to 19. 


(3.78) 


This completes the derivation of partial derivatives of the elements of matrix [Tf] 
needed for solving Eqs. (3.59) and (3.60). The partial derivatives of Xi, yi, xs, ys 
and zsi which are also needed to solve the two equations, are derived in the following 
subsection. 


3 . 2.2 Partial Derivatives of ri, 2/1, 2:25 2/2 and ^2 

The coordinates xi, yi, xs, ys and zs are given by Eqs. (3.14), (3.15), (3.16), (3.17) 
and (3.18) in terms of the mean values of the random variables Vi. These coordinates 
involve the parameters u and r. Prom Eqs. (3.63) and (3.64), 

u = u{Vi,Vsi. ■ ■ iVsjVs^Vdy.. ■ ,Vis) 

= r{Vi, Vs,. , Ve, Vg, Vg, • • - 1 ^la) 


r 
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Also 

dij = aij (Vi jVi, . . . ,V 5 ,Vs,Vg, . . . , V12) 

Keeping these in mind the partial derivatives of the coordinates at the mean values 
of the random variables are derived in the following: 

(0 

Xi =V$ cos u 


Hence, 


dxi 


dVe 

dxi 


dVi 

dxi 


dVi 


cosu — ri Sinn 


du 

dVe 


= —ri sin n 


= 0 


du 

dVi 


i = 1 to 5, 8 to 13 
i = 7, 14 to 19 


(3.79) 


(ii) 


Vi = Vgsinn 


Hence, 


dyi 


dVe 

dyi 

dVi 

dyi 


dVi 


= sinn + ricosn 


du 

Wb 


Ti cosn 


du 

dvi 


= 0 


i = 1 to 5, 8 to 13 
i = 7, 14 to 19 


(3.80) 


(iii) 


x-i. = ti3(aii cosr + 021 sinr) + 041 

Since aij is not a function of V13 = V2, but r is a function of V13 = r2, therefore, 


it follows that 
dX2 


dVe 

dx2 


dVi 


dx2 


= r 2 (— Oil sin r + 021 cos r) 
= r 2 (— Oil sin r + 021 cos r) 

da4i 


dr 

m 

dr 

Wi 


+M 


da 


11 


dVi 


cosr + 


da2i 


+ 


dVi 


dx2 


dVi 


ail cos r + 021 sin r + r 2 (— On sin r + O 21 cos r ) 


dVi 

i = 1 to 5, 8 to 12 
dr 


sinr) 


dVis 


0 


i = 7, 14 to 19 


(3.81) 
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(iv) 


y2 — Viz{ai2 COST + 022 sinr) + a42 


Hence, 

dy2 


dVe 

dy2 

dVi 


dy2 


^2 (-012 sin r + a22 cos r) 
r2 (— Oi2 sin r + 022 cos r) 


+ 


<9o42 


dr 

dr 

m 


dVi- 


+ r2{ 


dan 


dVi 


cosr + 


da22 


dVi 


sinr) 


i = 1 to 5 , 8 to 12 


dVn 

dy2 

dVi 


an cos r + a22 sin r + r2 (— O12 sin r + 022 cos r) 

0 


dr 


i = 7, 14 to 19 


( 3 . 82 ) 


(v) 


^2 = ^13 (oi3 COS r + 023 sin r) + O43 


Hence, 

dz2 


dV^ 

dz2 


dVi 


dZ2 


= ’'2 (—013 sin r + 023 cos r) 
= 'f'2 (— Oi3 sin r + a23 cos r) 


dr 

dr 

Wi 


+ 


5043 


+ r 2 ( 


dan 


dVi 


cosr + 


5 o 23 


dVi 


dV^2 

dZ2 

Wi 


ai3 cos r + 023 sin r + ?'2 (— 013 sin r + O23 cos r) 

0 


dVi 

i = 1 to 5 , 8 to 12 
dr 


sinr) 


5^13 
i = 7 , 14 to 19 


( 3 . 83 ) 


3.2.3 Partial Derivatives of the Elements fij 

Consider Eq. ( 3 . 58 ) which is a representation of Eq. ( 3 . 13 ). The derivatives of the 
elements fij in Eq. ( 3 . 58 ) are needed in solving Eqs. ( 3 . 59 ) and ( 3 . 60 ). Expressions 
for fij are written and differentiated in the following: 

(i) The elements /ij for j=l, 2, 3 are given by 


f^j = —aij sin r + 02j cos r 


j = 1 , 2,3 
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Hence, 

dfij 

dVi 


{ • \ / do/T'i dciO'i 

(-ay cosr - a 2 j smr)~ + ^ smr + cost) 

oVi oVi dVi 


j = 1,2,3. i = 1 to 5, 8 to 12 

df\j / dv 

— = (-aycosr-aajsinr)— j = 1,2,3. i = 6,13 

dhi 


dVi 


0 


j = 1,2,3. i = 7, 14 to 19 


(ii) The elements f 2 j for j=l,2,3 are given by 

f2j = O-Zj 


j = 1,2,3 


Hence, 


df2i da- 


'3i 


dVi dVi 


j = 1,2,3. i = 1 to 19 


(ui) 


fai =12-11 


Hence, 


dhi ^ dx2 dxi 
dVi ~ dVi dVi 


i = 1 to 19 


where ^ and are given by Eqs. (3.79) and (3.81) respectively. 


(iv) 


fz2 =y2-y\ 


Hence, 


dfz2 dy2 dyx 


i = 1 to 19 


dVi dVi dVi 

where and are given by Eqs. (3.80) and (3.82) respectively. 


(v) Finally, 
Hence, 


fzz = 22 


dfzz dz 2 


dVi dVi 

where ||| is given by Eqs. (3.83). 


i = 1 to 19 


(3.84) 


(3.85) 


(3.86) 


(3.87) 


(3.88) 



Chapter 3. Analysis of Mechanical Error in LOM 


58 


3.2.4 Partial Derivatives of the Parameters u and r 


The partial derivatives of the parameters u and r with respect to a variable Vi are the 
solution of Eqs. (3.59) and (3.60). The two equations are to be solved one by one for 
each variable Vi to obtain ■^r and i=l to 19. Rewriting the equations gives 


d{y 2 -yi) . ^ f ^ du ^ d(x 2 -xi) 

av, - - ») ““"av; av, ““ 


(a;2 -a;i)sinM^ 


fii fi 2 fn 


fn fi2 fi3 


fn fi 2 fiz 

/21 /22 /23 

+ 

fn fn fn 

+ 

hi hz hz 

/31 /32 fz3 


hi hz hz 


CO 


0 

0 


i = 1 to 19 


where prime denotes derivative with respect to Vi, i=l to 19. 

The terms in the above two equations contain the partial derivatives ^ and 
There is no multiplication of these two derivatives. Hence, the above two equations 
have linear terms of and ^ for a particular Vi. The two equations are not Hnear 
in u and r because the trigonometric functions of u and r are also present in the 
equations. But the two equations are solved for particular values of the parameters 
u and r. The parameters u and r in these equations correspond to the tangent line 
AC, as shown in Fig. (3.5). Their values are given hy u = ua and r = re- Therefore, 
the two equations are linear in ^ and ^ for a particular Vi. 

Even though Eqs. (3.59) and (3.60) are linear in ^ and it is a tedious process 
to express them explicitly in terms of these partial derivatives. This is because the 
derivatives of coordinates Xi, yi, X 2 , y2 and Z 2 involve and there are third 
order determinants involved, and moreover the process has to be repeated for the 
19 variables as the derivatives of the coordinates xi, yi, etc, are different for these 
variables. Therefore, an iterative technique is needed to solve the two equations. The 
two equations are solved by Newton-Raphson method [97]. Since it is difficult to 
compute the Jacobian matrix [J] analytically, the Jacobian is found numerically by 
the finite difference method. 

Solving the two equations by Newton-Raphson method give ^ and ^ at u = ua 
and r = tq. The two equations are solved one by one for each random variable Vi. 
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3.2.5 Partial Derivatives of the Elements of [Tq^] 

The dependent variables x, y and z are given by Eqs. (3.50), (3.51) and (3.52), and 
involve the elements pij. The matrix [p] having elements pij is given by Eq. (3.53). 
Rewriting it gives 

ricosu risinn 0 1 

ricostt risinw —li 1 

r2Cosr rasinr 0 1 

r2Cosr r2sinr —I2 1 

The derivatives of pij need the derivatives of matrices [Tq] and [Tq ] . The derivatives 
of the elements of [Tq], [Tq] and [p] are found in the present and following subsections. 

The matrix [Tq], in terms of the mean values of the random variables, is given by 
Eq. (3.5). Writing [Tq^] in terms of the variables Vi yields 

cos V5 sin 0 0 

sinV^sinVs -sinl^cosVs COSV4 0 
cosV^sinVs — COSV4COSV5 — sinV^ 0 
Vi V2 ^3 1 



(3.89) 




It follows that 

[r„‘] (3.90) 
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Differentiating [T^] with respect to the variables Vi yields 


d [Tb^] 


dVi 


d [T^] 


dV2 


d [Tb^] 


dVs 


d m] 


dV. 


d [To^] 


dVs 


d [T^] 


dVi 


0 0 0 0 
0 0 0 0 
0 0 0 0 
10 0 0 

0 0 0 0 
0 0 0 0 
0 0 0 0 
0 10 0 

0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 10 


0 


cos Oil sin (3i 
- sin Oil sin Pi 
0 

— sin Pi 


0 

— cos ai cos Pi 
sin Oi cos Pi 
0 


0 0 

— sin ai 0 

— cos ai 0 

0 0 


cos Pi 0 0 
sino!iCos/3i sinaisin/^i 0 0 
cos ai cos Pi cos ai sin /3i 0 0 
0 0 0 0 


= 0 i = 6 to 19 


m 


(3.91) 
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3.2.6 Partial Derivatives of the Elements of [T^] 

The matrix [Tq] in terms of the mean values of the random variables, is given by 
Eq. (3.6). Writing [Tq] in terms of the variables Vi yields 


qn2 

^0 


2 _ 


cos Vi 2 

siu V12 

0 

0 

sin Vii sin V 12 

— sin Vii cos V12 

cos Vii 

0 

cos Vji sin V 12 

— cos Vn cos Vi2 

— sin Vii 

0 


Vg 

V^io 

1 


It follows that 


rp2 

*^0 


= 7i(Vs,V,,...,Vn) 

Differentiating [Tq ] with respect to the variables Vi yields 


d [T 2 ] 


dVs 


d m 


dVg 


d [Ti] 


dVo 


^V^, 


d 


dVi2 


d [Til 


dVi 


0 0 0 0 
0 0 0 0 
0 0 0 0 
10 0 0 

0 0 0 0 
0 0 0 0 
0 0 0 0 
0 10 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 10 

0 0 0 0 
cos a 2 sin ^ — cos 02 cos — sin 02 0 
- sin 02 sin /32 sin 02 cos ^ -cos 02 0 


0 


0 0 
— sin/32 cos/32 0 0 

sinQ!2Cos/^ sinQ!2sin^2 0 0 
cos 02 cos /^ cos 02 sin ^2 0 0 
0 0 0 0 

i = 1 to 7, 13 to 19 


0 0 


(3.92) 


(3.93) 


(3.94) 
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3.2.7 Partial Derivatives of the Elements of [p] 


The matrix [p] in terms of the mean values of the random variables is given by 
Eq. (3.53). Differentiating it with respect to the variables Vi yields, 


d\p] ^ 

dVi dVi 


Vgcosu Vesinu 0 1 

Eg cos u Vgsiiiu —14 1 



Vis cos r Vs sin r 0 1 

Vs COST Vs sin r — V 4 1 



where as mentioned earlier, 


u 

r 



= u(V,V,...,V,V,V,...,V3) 

= Vi2) 


Keeping these in mind, it follows that 

du 


d\p\ 


dVi 


d[p] 


dVe 


d\p\ 


dVr 


— ri sinn- 


du 

dv.*- “ 

• n n 

— ri sinu-;;:77+ ricoswTrTT 0 0 


+ 


■dVi' dVi 

■ n n 

-r;jsmr-7^ r2Cosr^7r U U 

dVi oVi 

dr n n 

-r2smr— r2Cosr— U U 


in] 


ricosu risinu 0 1 

ricosw Tisintt —h 1 

0 0 0 0 


[n 

din] 


dVi 


0000 

du 

-ri sinu-; 


du 


i = 1 to 5 


+ cos It ri cosu f + sinu 0 0 

dVe dVs 

• n n 

-risiau-^ + COSU riCosu— + smu 0 0 
dVs . dVe 

. n n 

-r2smr-^ 0 U 

aV dV^ 

. n n 

-r^smr— rjccsr^ 0 0 

0 0 


in] 


0 0 

cos Oil sin /3i cos oci cos Pi sin ai 0 

0 0 0 0 

0 0 0 0 


(3.9.51 
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d\p] 

dVi 


d\p] 


dV^z 


d\p] 

dVi4 

d\p] 

dVi 

Hence, 


+ 


. du r. 

-nsmu— ricosu— 0 0 

u\/% 

. du du ^ ^ 

-r^Bxnu— r.cosu— 0 0 

.dr 9 t 

-r2smr— rzcosr— 0 0 
dVi dVi 

dr n n 

-r2smr— r2COsr— 0 0 
dVi _ dVi 

0 0 0 o' 

0 0 0 0 


[^0^] 


r2Cosr r2sinr 0 1 

rzcosr r2sinr —I 2 1 
du du 


d[Ti] 


-ri sin u 


-ri sm u- 


dV,z 

au 


Ti cos u 


n cosu 


dV^z 

du 

dWz 


dV^z 
dr 

—r 2 sin r + cos r rzcosr 


dVi 

0 0 

0 0 
dr 


i = 8 to 12 


wi 

+ sinr 0 0 


0 


dVi3 ' "" ■ aji3 

dv ot 

r 2 sin rTr— — h cos r rz cos r-^ — h sin r 0 0 
dViz dViz 

0 0 0 0 

0 0 0 0 

0 0 0 0 

— cos 012 sin ^2 cos az cos /?2 sin 0:2 0 

i = 15 to 19 


(3.96) 


[p] = P (^i> • • • » (3.97) 

The partial derivatives of the elements of [p] have- been derived as in Eqs. (3.95) 
and (3.96) above. 


3.2.8 Partial Derivatives of u and w 

The parameters u and w on the surface q(u,w) are given by Eq. (3.39) and (3.40) 

x*(u,iy) = 0 
y*{u,w) = 0 
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where, from Eqs. (3.36) and (3.37) 


[ x{u,w) y{ 

u,w) z{u,w) 1 j 

[T] 

\ 1 — u u 

■p(0,0) p(0,l)' 

1 — w 

L J 

_p(l,0) p(l,l)_ 

w 


[T] 


(3.98) 


Hence u and w are given by 


1 I 

o' 

'ft 

p(o, 1) 

1 —w 

(T| = 

r _ . . 1 

1 — u u J 

_ p(l,0) 

P(l.l) _ 

w 

1 — 1 

0 

0 


From Eq. (3.97), [p] — p {Vx, V 2 , . . . , Vu). Equation (3.34) gives [T] in terms of the 
mean values of variables Vi. 


[T] —T{a3,bz,cz,6,(p) 

[T], in terms of the variables Vi, is given by 

[T]=T{Vxs,Vie...,Vx,) (3.100) 

Hence, from Eq. (3.99) it follows that 


u = u{Vx,V2,...,Vx^) (3.101) 

w = wiVi,V2,...,Vxg) (3.102) 


where u and w are parameters of the surface q(u,w). 
Differentiating Eqs. (3.39) and (3.40) yields 

dx* 


dVi 

dy* 

dVi 


= 0 


= 0 


where, using Eq. (3.100), 


+[ X 


y z 


1 ; 


m 

dVi 


dr* 

dp* 

dz* 

0 


dr 

dp 

S 0] 

dVi 

dVi 

dV, J 


. dVi 

dVi 

dVi J 

dr* 

dp* 

0 

dVi J 


dr 

dp 

0 

dVi J 

. Wi 

dVi 


. m 

dVi 


[T] 

[T] 


(3.103) 

(3.104) 


i = 1 to 14 


(3.105) 


i = 15 to 19 


Partial derivatives of the dependent variables x, y, and z, and those of the elements 
of [T] are derived in the next two subsections. For each Vi, i=l to 19, Eqs. (3.103) 
and (3.104) are solved for ^ and |^. This gives ^ and ^ for all i, i=l to 19. 

Although Eqs. (3.103) and (3.104) are linear in ^ and || it is quite cumbersome 
to solve them by collecting the coefficients of ^ and So an iterative technique is 
adopted. They are solved by Newton-Raphson method with the Jacobian computed 
numerically [97]. 
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3.2.9 Partial Derivatives of the Elements of [T] 

The transformation matrix [T] in terms of the mean values of the random variables 
is given by Eq. ( 3 . 34 ). Writing [T] in terms of the variables Vi yields 


[T] = 


COS Vis cos Vi 9 

-sinlig 

— sin Vis cos ViQ 

0 ■ 

cos Vis sin Vig 

COS ^19 

— sin Vis sin ^19 

0 

sin Vis 

0 

cos ids 

0 

— cos Fis (Ei 5 cos Ei9 + Vi 6 sin Vig 
+Vi7tanyi8) 

Vis sin Vis — Vie cos Vis 

sin Vi 8 (Vi 5 cos Vis + Vie sin Vig 
+Vi7tanVi8) 

1 


( 3 . 106 ) 


Therefore, as in Eq. ( 3 . 100 ) 


[T]=r(Fa5,Vi6,...,Ei9) 


Hence, 


d[T] 

dVi 


\m 


m 

avis 


m 


m 

91^16 


m 


m 

9^7 


m 


am 


ai^is 


Im 


am 

dVis 


m 


= 0 


i = 1 to 14 


0 0 0 0 

0 0 0 0 

0 0 0 0 

— cos 6 cos <p sin (f> sin 6 cos 4 > 0 



0 

0 


0 

0 


0 

0 


0 

0 


0 

0 


0 

0 

8 

1 

3 in</> - 

- cos^ 

sin 

0sin0 

0 

0 

0 

0 

0 ' 



0 

0 

0 

0 



0 

0 

0 

0 



— sin^ 

0 sin 

9 t& 3 x 6 

0 




— sin 6 cos <j> 

0 

— cos 9 cos <p 

0 ‘ 

— sin 6 sin (j) 

0 

— cos 0 sin ^ 

0 

cos 9 

0 

— sin^ 

0 

sin 9 {as cos (p + hs sin (p) 
— C3 cos 9 

0 

cos 9 {as cos 0 + 63 sin (p) 
C3 sin 0(1 + sec^0) 

0 

— cos 9 sin Ip 

— cos 4 > 

sin 0 sin 0 

0 

cos 9 cos <p 

— sm<p 

— sin 0 cos <p 

0 

0 

0 

0 

0 

— cos 9 {-az sin ^ 
+63 cos (^) 

as cos <p-Vhz sin 4 > 

sin0(— a3 sm0 
+63003^) 

0 


( 3 . 107 ) 
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3.2.10 The Influence Coefficients (^)^and 

The dependent variables [a: y z 1 ] = [9] are given by Eqs. (3.50), (3.51) and 
(3.52). It follows that 


q[j] =PijO- ~ ■“)(! ~ +P 2 i(l - u)w + P3jU{l ~ w) +P4juw 


where, as has been mentioned earlier 


j = 1,2,3 
(3.108) 


Pij = Pij{Vi toVu) 


Hence, differentiating 
dq[j] dpij 


u = w(Vi,F2 ,.-.,1 i9 ) 

W = w{Vi,V2, . . . ,Vi9) 
in Eq. (3.108) with respect to Vi gives 


dVi 


dVi 


(1 - w)(l -w) + ~ 'w)w + ^^■«(1 - li^) + 


dVi 


dVi 


+ [-Pli(l - w) - P2jW +P3i(l - w) +P4jW[ 


du 


+ [-Pli(l - u) +P2jil-u)-pyU + P4jU] 


9q[j. 

dVi 


dVi 

dVi 

du 


dVi 


i = 1 to 14 


-P2jVJ +J>3i(l - W) +P4jW] — 

, - . , . , dw 

+ [-Pli(l - u) +P2/(1 - u) -P3ju+P4ju] — 


i = 15 to 19 

(3.109) 

where and |^, i=l to 19, have been derived earlier. Substituting, we get the 

partial derivatives ^ and i=l to 19, from Eqs. (3.109) above. The influence 
coefficients (^)^, i=l to 19, are then obtained by squaring these 

partial derivatives. 


3.3 Numerical Results 

Numerical results are obtained with the following input values (See Fig. 3.1): 


tti = 0.1 m 

61 = 0.1 m 

Cl = 0.1 m 

ai = 0.5^ 


A = 1-0° 

ri = 0.03 m 

li = 0.50 m 

02 = 1.04 m 

62 = 0.1 m 

C2 = 0.1 m 

0 

P 

rH 

1 

II 


A = -0.5° 

r2 = 0.03 m 

^2 = 0.05 m 

C3 = 0.6 m 

11 

p 

0 

0 

11 

0 


Absolute tolerances 

on x-coordinates: 

€1 = eg = ei5 

= 0.0001 m 
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Absolute tolerances on y-coordinates: 62 = eg = eig = 0.0001 m 

Absolute tolerances on z-coordinates: 63 = do = ei? = 0.0001 m 

Absolute tolerances on angles: d = £5 = en = ei2 = Cis = ejg = 0.008° 

Tolerance per unit nominal length on radii and lengths of the two rollers: 

€7 = eg = ei3 = ei4 = 0.005 m/m 

The nominal values of the angles 0 and (f) are both zero, so that the laser beam 
is vertical. The standard deviation cr of error at a point on the work surface is given 
by [110] 

cr = VDM + D[y] + D[z| 

= VSSUM, say. 

where 

DSUM = D[x] + D[y] + D[z] 

Hence, 

DSUM = 0-2 

i.e., DSUM is the square of the standard deviation of error at a point. 

The variances D[x], D[y] and D[z] are obtained at several points on q(u,w), the 
bilinear surface ABDC. These points are obtained by varying the parameters u and 
w in the ranges, 0 < « < 1.0 and 0 < tw < 1.0 with increments of 0.1. Table 3.1 
tabulates the variances D[x], D[y] and D[z], and their sum DSUM for these values of 
u and w. 

Prom the table it is observed that the variances D[x], D[y] and D[z] are of the 
same order. The variance D[z] varies significantly over the work surface compared to 
D[x] and D[y]. The variance D[y] doesn’t vary over the work surface. The influence 
coefl^cient = 1.0 throughout the work surface as can be easily seen; and 

(^) > other i’s, are negligible compared to it. Hence, the variance D[y] is equal 
to the variance pig throughout the work surface and, therefore, it doesn’t vary. 

The variance D[x] varies slightly over the work surface. The influence coefficient 
~ work surface, as can be easily seen. Here Vis is the x-coordinate 

of the source of the laser beam. For mfVig] = 0° and m[Vig] = 0°, increasing Vig 
increases the x-coordinate, whereas increasing Vig doesn’t change the roordinates of 
a point. That is why is significantly large but (^)^, (^)m 

are neghgibly small. The influence coefficients varies sfightly over the work 

surface causing small variation of the variance D[x] over it. Table 3.2 lists the influence 

coefficients at (u=1.0, w=1.0) for the input values mentioned above. 

Prom the error analysis above, the maximum value of the sum of variances DSUM 
is found at (u=1.0, w=l.G). 
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QX 

m 

m 

mgn 

D[y] 

DH 

DSUM 

10-® m 2 

Qk 

m 

Qy 

m 

D [ x ] 

Dly ] 

Dfz ] 

10-® m 2 

DSUM 

10-^mP 

o.ob 

' o.od 

1.5896 

1.1111 

2.50028 


6.50 

6.50 

1.5918 

1.1111 

1.26460 

3.9676 

0.00 

0.10 

1.5887 

1.1111 

2.50050 

5.2003 

0.50 

0.60 

1.5923 

1.1111 

1.27104 

3.9744 

0.00 

0.20 

1.5878 

1.1111 

2.50074 

5.1997 

0.50 

0.70 

1.5927 

1.1111 

1.27862 

3.9825 

0.00 

0.30 

1.5870 

1.1111 

2.50098 

5.1991 

0.50 

0.80 

1.5932 

1.1111 

1.28733 

3.9916 

0.00 

0.40 

1.5861 

1.1111 

2.50124 

5.1984 

0.50 

0.90 

1.5936 

1.1111 

1.29716 

4.0019 

0.00 

0.50 

1.5852 

1.1111 

2.50150 

5.1978 

0.50 

1.00 

1.5941 

1.1111 

1.30811 

4.0133 

0.00 

0.60 

1.5843 

1.1111 

2.50178 

5.1972 

0.60 

0.00 

1.5896 

1.1111 

1.29991 

4.0007 

0.00 

0.70 

1.5834 

1.1111 

2.50206 

5.1966 

0.60 

0.10 

1.5903 

1.1111 

1.30042 

4.0019 

0.00 

0.80 

1.5825 

1.1111 

2.50235 

5.1960 

0.60 

0.20 

1.5910 

1.1111 

1.30202 

4.0042 

0.00 

0.90 

1.5817 

1.1111 

2.50265 

5.1954 

0.60 

0.30 

1.5918 

1.1111 

1.30471 

4.0076 

0.00 

1.00 

1.5808 

1.1111 

2.50296 

5.1949 

0.60 

0.40 

1.5925 

1.1111 

1.30848 

4.0121 

0.10 

0.00 

1.5896 

1.1111 

2.05010 

4.7508 

0.60 

0.50 

1.5932 

1.1111 

1.31332 

4.0176 

0.10 

0.10 

1.5890 

1.1111 

2.05036 

4.7505 

0.60 

0.60 

1.5939 

1.1111 

1.31921 

4.0242 

0.10 

0.20 

1.5884 

1.1111 

2.05078 

4.7503 

0.60 

0.80 

1.5953 

1.1111 

1.33415 

4.0406 

0.10 

0.30 

1.5878 

1.1111 

2.05136 

4.7502 

0.60 

0.90 

1.5960 

1.1111 

1.34317 

4.0503 

0.10 

0.40 

1.5871 

1.1111 

2.05210 

4.7503 

0.60 

1.00 

1.5968 

1.1111 

1.35322 

4.0611 

0.10 

0.50 

1.5865 

1.1111 

2.05300 

4.7506 

0,70 

0.00 

1.5896 

1.1111 

1.45008 

4.1508 

0.10 

0.60 

1.5859 

1.1111 

2.05405 

4.7511 

0.70 

0.10 

1.5906 

1.1111 

1.45041 

4.1521 

0.10 

0.70 

1.5853 

1.1111 

2.05526 

4.7517 

0.70 

0.20 

1.5916 

1.1111 

1.45159 

4.1543 

0.10 

0.80 

1.5847 

1.1111 

2.05663 

4.7524 

0.70 

0.30 

1.5926 

1.1111 

1.45360 

4.1573 

0.10 

0.90 

1.5840 

1.1111 

2.05815 

4.7533 

0.70 

0.40 

1.5935 

1.1111 

1.45643 

4.1611 

0.10 

1.00 

1.5834 

1.1111 

2.05981 

4.7543 

0.70 

0.50 

1.5945 

1.1111 

1.46009 

4.1657 

0.20 

0.00 

1.5896 

1.1111 

1.69996 

4.4007 

0.70 

0.60 

1.5955 

1.1111 

1.46456 

4.1712 

0.20 

0.10 

1.5893 

1.1111 

1.70034 

4.4007 

0.70 

0.70 

1.5965 

1.1111 

1.46984 

4.1774 

0.20 

0.20 

1.5889 

1.1111 

1.70121 

4.4012 

0.70 

0.80 

1.5975 

1.1111 

1.47592 

4.1845 

0.20 

0.30 

1.5886 

1.1111 

1.70257 

4.4022 

0.70 

0.90 

1.5985 

1.1111 

1.48279 

4.1924 

0.20 

0.40 

1.5882 

1.1111 

1.70441 

4.4037 

0.70 

1.00 

1.5994 

1.1111 

1.49045 

4.2010 

0.20 

0.50 

1.5878 

1.1111 

1.70674 

4.4057 

0.80 

0.00 

1.5896 

1.1111 

1.70032 

4.4011 

0.20 

0.60 

1.5875 

1.1111 

1.70954 

4.4081 

0.80 

0.10 

1.5909 

1.1111 

1.70044 

4.4024 

0.20 

0.70 

1.5871 

1.1111 

1.71281 

4.4111 

0.80 

0.20 

1.5921 

1.1111 

1.70104 

4.4043 

0.20 

0.80 

1.5868 

1.1111 

1.71655 

4.4144 

0.80 

0.30 

1.5934 

1.1111 

1.70214 

4.4066 

0.20 

0.90 

1.5864 

1.1111 

1.72075 

4.4183 

0.80 

0.40 

1.5946 

1.1111 

1.70372 

4.4094 

0.20 

1.00 

1.5861 

1.1111 

1.72541 

4,4226 

0.80 

0.50 

1.5959 

1.1111 

1.70578 

4.4128 

0.30 

0.00 

1.5896 

1.1111 

1.44986 

4.1506 

0.80 

0.60 

1.5971 

1.1111 

1.70831 

4.4165 

0.30 

O.lO 

1.5895 

1.1111 

1.45037 

4.1510 

0.80 

0.70 

1.5984 

1.1111 

1.71132 

4.4208 

0.30 

0.20 

1.5894 

1.1111 

1.45172 

4.1523 

0.80 

0.80 

1.5996 

1.1111 

1.71480 

4.4255 

0.30 

0.30 

1.5894 

1.1111 

1.45391 

4.1544 

0.80 

0.90 

1.6009 

1.1111 

1.71875 

4.4307 

0.30 

0.40 

1.5893 

1.1111 

1.45692 

4.1573 

0.80 

1.00 

1.6021 

1.1111 

1.72315 

4.4364 

0.30 

0.50 

1.5892 

1.1111 

1.46076 

4.1610 

0.90 

0.00 

1.5896 

1.1111 

2.05066 

4.7514 

0.30 

0.60 

1.5891 

1.1111 

1.46540 

4.1656 

0.90 

0.10 

1.5911 

1.1111 

2.(»056 

4.7528 

0.30 

0.70 

1.5890 

1.1111 

1.47086 

4.1710 

0.90 

0.20 

1.5927 

1.1111 

2.05063 

4.7544 

0.30 

0.80 

1.5889 

1.1111 

1.47711 

4.1771 

0.90 

0.30 

1.5942 

1.1111 

2.05086 

4.7561 

0.30 

0.90 

1.5888 

1.1111 

1.48416 

4.1841 

0.90 

0.40 

1.5957 

1.1111 

2.05125 

4.7580 

0.30 

1.00 

1.5887 

1.1111 

1.49199 

4.1918 

0.90 

0.50 

1,5972 

1.1111 

2.05179 

4.7601 

0.40 

0.00 

1.5896 

1.1111 

1.29982 

4.0005 

0.90 

0.60 

1.5987 

1.1111 

2.( e 250 

4.7623 

0.40 

0.10 

1.5898 

1.1111 

1.30041 

4.0013 

0.90 

0.70 

1.6003 

1.1111 

2.( K 336 

4.7647 

0.40 

0.20 

1.5900 

1.1111 

1.30210 

4.0032 

0.90 

0.80 

1.6018 

1.1111 

2.05438 

4.7673 

0.40 

0.30 

1.5902 

1.1111 

1.30488 

4.0061 

0.90 

0.90 

1.6033 

1.1111 

2.05555 

4.7700 

0.40 

0.40 

1.5903 

1.1111 

1.30873 

4.0102 

0.90 

1.00 

1.6048 

1,1111 

2.05687 

4.7728 

0.40 

0.50 

1.5905 

1.1111 

1.31366 

4.0153 

1.00 

0.00 

1.5896 

1.1111 

2.50110 

5.2018 

0.40 

0.60 

1.5907 

1.1111 

1.31964 

4.0214 

1.00 

0.10 

1.5914 

1.1111 

2.50088 

5.2034 

0.40 

0.70 

1.5909 

1.1111 

1.32668 

4.0287 

1.00 

0.20 

1.5932 

1,1111 

2.50067 

5.2050 

0.40 

0.80 

1.5910 

1.1111 

1.33475 

4.0369 

1.00 

0.30 

1.5950 

1.1111 

2.50048 

5.2066 

0.40 

0.90 

1.5912 

1.1111 

1.34386 

4.0462 

1.00 

0.40 

1.5968 

1.1111 

2.50029 

5.2082 

0.40 

1.00 

1.5914 

1.1111 

1.35400 

4.0565 

1.00 

0.50 

1.5986 

1.1111 

2.50012 

5.2098 

0.50 

0.00 

1.5896 

1.1111 

1.24983 

3.9506 

1.00 

0.60 

1.6003 

1.1111 

2.49995 

5.2114 

0.50 

0.10 

1.5901 

1.1111 

1.25043 

3.9516 

1.00 

0.70 

1.6021 

1.1111 

2.49980 

5.2131 

0.50 

0.20 

1.5905 

1.1111 

1.25221 

3.9538 

1.00 

0.80 

1.6039 

1.1111 

2.49966 

5.2147 

0.50 

0.30 

1.5910 

1.1111 

1.25518 

3.9572 

1.00 

0.90 

1.6058 

1.1111 

2.49953 

5.2164 

0.50 

0.40 

1.5914 

1.1111 

1.25931 

3.9618 

1.00 

1.00 

1.6076 

1.1111 

2.49941 

5.2181 


Table 3.1: Variances and Their Sum at Several Points in uw Parametric Space on the 


Bilinear Surface ABDC in LOM 
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The Influence Coefficients 

i 



i '1 

\9Vi)m 



1 

-7.553e-14 

1.981e-13 

-4.784e-09 

2 

4.797e-13 

1.982e-13 

7.651e-06 

3 

-7.335e-14 

-3.546e-13 

4.383e-04 

4 

-7.548e-14 

1.982e-13 

-3.443e-07 

5 

-7.548e-14 

1.982e-13 

-4.869e-ll 

6 

-7.327e-14 

-3.546e-13 

4.382e-04 

7 

-7.546e-14 

1.982e-13 

-2.400e-15 

8 

-6.110e-12 

-2.942e-13 

1.372e-02 

9 

-7.553e-14 

1.981e-13 

-7.651e-06 

10 

-7.769e-14 

1.964e-13 

-4.383e-04 

11 

4.797e-13 

-3.570e-13 

3.940e-06 

12 

-7.552e-14 

1.982e-13 

7.192e-06 

13 

-8.3l5e-14 

-3.811e-13 

9.997e-01 

14 

-3.896e-13 

3.760e-12 

-1.490e-09 

15 

l.OOOe-fOO 

1.369e-13 

-1.372e-02 

16 

2.392e-13 

l.OOOe-00 

-1.734e-02 

17 

-7.546e-14 

1.982e-13 

-2.400e-15 

18 

4.787e-01 

1.689e-13 

-6.570e-03 

19 

-7.546e-14 

1.982e-13 

-2.400e-15 


Table 3.2: The Influence Coefficients at (u = 1.0, w = 1.0) in uw Parametric Space 
on the Bihnear Surface ABDC in LOM 

The three-sigma bands of error in tracing several lines on the bilinear surface 
q(u,w) is obtained. The three-sigma band is the band of 3cr and -3cr values of error at 
a point. Figure 3. 11 shows the three-sigma band for a horizontal line in uw parametric 
space from (.1, .1) to (.9, .1). Figure 3.12 shows the three-sigma band for a vertical 
line from (.9, .1) to (.9, .9). Figure 3.13 shows the three-sigma band for an inclined 
line from (.1, .1) to (.9, .9). The three-sigma bands for the horizontal and inclined 
lines decreases and then increases. The three-sigma band for the vertical hne is 
almost constant. This is because the sum of variances DSUM doesn’t vary much in 
the vertical direction in uw parametric space as can be seen from Table 3.1. From 
Table 3.1 and the three-sigma band it is seen that the overall error doesn’t vary much 
over the work surface and it is immaterial where on the work smface, the contour is 
drawn. 




Error (in mm) 
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Parameter s along the Line 


Figure 3.11: The 3a band for a Horizontal Line in uw Parametric Space from (.l,.l) 



OA 0.6 

Parameter s along the Line 


Figure 3.12: The 3a band for a Vertical Line in uw Parametric Space from (•9,T) 









Error (In mm) 
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Figure 3.13: The 3cr band for an Inclined Line in uw Parametric Space from (.l,.l) 
to (.9, .9) 



Chapter 4 


ANALYSIS OF MECHANICAL 
ERROR IN EDM 


4.1 Geometric Modeling for Description of the FDM 
Process 

The Fused Deposition Modeling (FDM) process comprises of an extruder head moving 
in the horizontal plane (Fig. 1.3). The extruder head slides along a link in X direction. 
Figure 4.1 shows the motion of extruder head on the head rails. This hnk itself shdes 
along another link in Y direction. When the part is built progressively, the Z-stage 
platen moves downwards. The tip of the extruder head deposits material on the 
Z-stage platen. It draws the contour of a slice, and fills the contour with roads 
[63, 95, 13, 45]. The FDM process has been discussed in detail in Section 1. 1.2.2. 

The geometric model of the FDM process is as follows. There is a kinematic 
chain consisting of links and other members causing the movement of the extruder 
head in FDM (Fig. 4.2). The extruder head shdes on member 3. The sub-assembly 
consisting of member 3 and the extruder head shdes on hnk 1. The members DE 
and EG constitute a rigid hnk. The members DE and EG are called the sliding 
member 2 and member 3 respectively. The shding member 2 is taken to be paraUel 
(lO Zq direction. There is a shding pair between hnk 1 and the sliding member 2. It 
is called pair 12. 

The present model is being done for a particular layer of the part being built. 
Consider the Z-stage platen for that layer. A base frame XoYoZo is attached at the 
front end of the Z-stage platen such that the surface of the foam foundation is the 
^oFb plane (Fig. 4.2). The frame XoVo^o is fixed in space, it does not move with the 



Chapter 4. Analysis of Mechanical Error in PDM 


73 



Figure 4.1: Motion of Extruder Head in FDM 


platen. 

Ideally, link 1 and Yq direction should be parallel, and member 3 and Xq direction 
should be parallel. However, due to mechanical errors the members 1 and 3 are 
oblique in space. A frame XiYiZi is attached at the front end A of link 1 such that 
Yi axis is along link l(Fig. 4.2). The frame XjYiZi is obtained by a rotation ai 
about Ao-axis, a rotation /?i about ko-axis and a translation (oi,fci,ci) in XqYqZq 
successively. 

Another frame X 2 Y 2 Z 2 is attached at the left end E of member 3 such that the 
A 2 -axis is along member 3 (Fig. 4.2). The frame A 2 k 2‘^2 is obtained from XiYiZi as 
follows. The frame X 2 Y 2 Z 2 is obtained by a translation AD along link 1, a translation 
DE along link 2, a rotation 0:2 about a direction parallel to Xi-axis and passing 
through E and rotation ^ about a direction parallel to ki-axis and passing through 
E successively. 

The extruder head is modeled by two members 4 and 5 such that they are parallel 
to k 2 and Z 2 directions respectively. As shown in Fig. 4.2, the line segments HI and 
IQ are members 4 and 5 respectively. The nozzle tip is at the point Q. The members 
4 and 5 are rigidly connected to each other and, therefore, make a single link. The 
extruder head being a rigid body can always be modeled by two sudi members parallel 
to ^2 and Z 2 directions. There is a sliding pair between Unk 3 and member 4. This 
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Figure 4.2: Geometric Model of the FDM Process 


pair is called pair 34. 

The kinematic chain is an open link mechanism. The left and right head-rails are 
fixed links and, therefore only the head-rail on the left is considered. The nominal 
lengths of various links and members are as follows. The sHding member 2 is con- 
strained to move between two points B and C on link 1 such that the lengths of AB 
and ~BD are hfixed and hmax respectively. At any position D of the sliding member 2 
the lengths of link 1 is 1 % fixed + h 

The nominal length of the sliding member 2 is I 2 . The sliding member 4 is con- 
strained to move between two points F and G on member 3 such that the lengths 
of 'EF and TG are k/ixed and kmax respectively. At any position H of the sliding 
member 2 the length of member 3 is hfixed + h- The nominal length of members 
4 and 5 are I 4 and h respectively. The dimensions are taken in S.I. units, i.e., the 
lengths in meters and the angles in radians. 
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4.2 Stochastic IVIod.©! of the FDJVE Process 

Chapter 2 discusses stochastic modeling of the RP process. Following that approach, 
the stochastic modeHng of FDM is done as follows. The sliding members 2 and 4 are 
driven by belt and pulley drive mechanism. There is an error in positioning 2 and 
4 by the corresponding drive mechanism. Let ci be the absolute error, in meters, in 
positioning the sliding member 2 along link 1 by the drive mechanism. Let €2, €4 and 
65 be the tolerances per unit nominal length in m/m on members 2, 4, 5 respectively. 
Let €3 be the absolute error, in meters, in positioning the sliding member 4 along 
member 3 by the drive mechanism. The actual lengths Ri of links are given by 


Ri 

= hfixed + + Cl 

( 4 . 1 ) 

R2 

= ^2(1+ ^2) 

(4.2) 

Rz 

= hfixed + h + ^3 

( 4 . 3 ) 

R4 

= 14(1 -(- 64 ) 

( 4 . 4 ) 

jRs 

= 15(1 +65) 

( 4 . 5 ) 


where Ri are now random variables. 

The effective clearance error e along a link axis in a spatial revolute joint is given 
by Eq. (2.8). It is as follows. 


e = 


rcosa 
cos 7 


where r and a are parameters, random in nature as illustrated in Fig. 2.5. 7 is the 
angle made by the link with a plane perpendicular to the socket axis. 

Although the pairs 12 and 34 are sliding pairs, they comprise of the cylindrical 
elements socket and pin. Therefore the clearance in them is given by Eq. (2.8), as 
above. 

Let eg be the effective clearance error in pair 12 along member 3, i.e., along X 2 
direction (Fig. 4.3). Therefore, 1 

( 4 . 6 ) 

COS76 

where rj and ctb are the parameters for eg. They are random in nature with their 
ranges given by 0 < rj, < C12 and 0 ^ 27r; where C12 is the radial clearance in pair 

12. It is assumed that the radial clearances in both the left and right head-rails are 
equal, and they are C12. 75 is the angle made by member 3 with a plane perpendicular 
to link 1. 

Let €7 be the effective clearance error in pair 12 along the sliding member 2, i.e., 
along Zo direction (Fig. 4.3). Therefore, 

rg cosag (4.7) 


€7 = 


cos 7a 
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Figure 4.3: Equivalent Kinematic Chain 


where and aa are the parameters for e^. They are random in nature with their 
ranges given by 0 < To < Ci 2 and 0 < oia < 27r. 7 a is the angle made by the sliding 
member 2 with a plane perpendicular to hnk 1. 

Let €8 be the effective clearance error in pair 34 along the sliding member 4 which 
lies along ¥2 direction (Fig. 4.3). Therefore, 

Vc cos Qc 

cs = 

cos 7c 

where rc and ckc are the parameters for eg. They are random in nature with their 
ranges given by 0 < rc < C34 and 0 < qjc < 27r; where C34 is the radial clearance in 
pair 34. % is the angle made by the sliding member 4 with a plane perpendicular 
to member 3. Since the shding member 4, being parallel to T^-axis, Hes in a plane 
perpendicular to member 3, therefore, 

7o = 0° 


Hence, 


es = Te cos ac 


(4.8) 
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direction (Fig. 4 . 3 ). Therefore, 


rj cos ad 
cos 7 d 


where tj, and cxd are the parameters for eg. They are random in nature with their 
ranges given by 0 < < C34 and 0 < < 27 r. 7^ is the angle made by member 5 

with a plane perpendicular to member 3 . Since member 5 , being parallel to Z2-axis, 
lies in a plane perpendicular to member 3 , therefore, 


Hence, 


7 <f = 0° 


£9 = Td cos ad 


The effective clearance errors q, i= 6 to 9 , are random variables. 

The extruder head is fixed by screws to the slider which slides on member 3 . This 
head attachment causes an error in Z2 direction. Let ejo be the absolute error, in 
meters, in Z2 direction at the head attachment. This modifies the length of member 


The nozzle tip is screwed to the extruder head. This nozzle tip attachment too 
causes an error in Z2 direction. Let en be the absolute error, in meters, in Z2 direction 
at the nozzle tip attachment. The error ey too modifies the length of member 5 . The 
errors €10 and en are also random variables. 

Let £12, ei3, 614 and 615 be the absolute tolerances, in radians, on angles ai, /?i, a2 
and 132 respectively. The actual angle is then given by the sum of the nominal angle 
and its corresponding tolerance. The actual angles !are now random variables. 

So, in all, there are fifteen random variables involved in the present model. If 
these are denoted by Pi, V2, . . . , Vis then 

Vi = Ri i = 1 to 5 

Vi = €i i = 6 to 11 

Vi2 = ai + ei2 (4-10) 

Via = /?! + cia ' 

Vi4 = 0!2 + ei4 

Vis = yda + ^15 

These random variables are independent of each o^her. 

The equivalent lengths of the members in terms of the random variables are given 

as follows: 

i?ii = Ri = Vi 

R21 = Hz + e? = '^2 + ^4 

R21 = R3 + eg = V3 + ^6 

i ?41 = i ?4 + €8 = ^4 + ^8 

iJsi = i ?5 + €9 + 610 + fill = Vs+Vg + Vio + Vii 


( 4 . 11 ) 
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where Rn the equivalent length of the member. The kinematic chain with 
equivalent lengths of members is shown in Fig. 4.3. The above is a good approximation 
for Ri\ and Rzi- When members 1, 2 and 3 are perpendicular to each other, as in 
the ideal case, then the expressions for i? 2 i and R 31 holds exactly. 

When the first layer of the part is being built, the nozzle tip is buried .01 in into 
the foam foundation on the Z-stage platen. The foam foundation deforms as the first 
few layers are made, and the layer comes out as traced by the nozzle tip. So the error 
in the Z-stage platen does not give rise to any error in the part. 

4.2.1 Mean Values of the Random Variables 

Prom Section 2.2.1, the mean values of actual member lengths are given by 

^[1^] “ Wfixed “I" ^1 

m[V2] = h 

m[V 3 ] = hfixed + h 

m\y^ — I 4 
mlVs] = k 

Section 2.2.2. which discusses effective clearance error in a spatial revolute joint, 
yields 

m[Vi\ =0 i = 6, 7, 8, 9 

€10 and Cii are absolute error on member lengths. They are normally distributed 
with zero nominal value. Hence 


m[Vi\ =0 i = 10,11 

Section 2.2.4 yields 


m[Vi2] 

= ai 

m[Vi3] 

= A 

m[Vi4] 

= 0:2 

mjVis] 

= A 


4.2.2 Variances of the Random Variables 

Variances pi of the random variables Vi are discussed below. Let e< be the tolerance 
per unit nominal length of a hnk. From Section 2.2.1, the variance of the actual 

length of a link is given by 

Pi = Wif 

where k is the nominal length of the link. 
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Therefore, if is the absolute tolerance on the nominal length of a link then the 
variance is given by 

Pi = (et/3)^ (4.12) 

Hence for the random variables in the present model, it follows that 

Pi = i^/Sy i = 1,3,10,11 

Pi = i = 2,4,5 

Section 2.2.2 yields 


Pe = 

^12 


4 cos 7a 

P7 = 

42 


4 cos 76 

P% == 

^34 

4 


ci 

P9 = 

4 


From Section 2.2.4, 

Pi = (ei/3)2 i = 12, 13, 14, 15 

The angles 7a and jb, which appear in the expressions for pe and p 7 , respectively, are 
derived in the following subsections. 


4.2.3 The Angle ja 


7a is the angle made by the sliding member 2 with a plane perpendicular to link 1. 
A frame X 5 Y 5 Z 5 parallel to XiYiZi is attached at D as shown in Fig. 4.4. 

So 7o is the angle made by link 2, i.e., the line segment DE with the X^Zs-plane. 
Figure 4.5 shows a frame XYZ with 0 as its origin. A line OP from a point P(x,y,z) 
makes an angle 7 with the XZ-plane. The angle 7 is given by 


7 = tan ^ 


■\/z^ + 


(4.13) 


Using this, the angle 7a can be derived as follows. The point E in X 2 Y 2 Z 2 frame, 
being its origin, is given by 

[e 2]=[0 001] 

! 

Therefore, the point B in X5F5.Z5 frame is given by 


[es] = [ 0 0 0 1 j 
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Z-stage platen 


Figure 4 . 4 : The Angles 7^ and 75 


The transfornaation [Ts ] involves a translation i?2i along Zio direction. So the trans- 
formation [Tq] of frame X5Y5ZS with respect to XoYoZq is found and using it the 
transformation [Tg ] is derived. The transformation [Tq] is given by 



Or, 



Hence, 


[t/] = [If] [If] 


-1 


( 4 . 14 ) 


The transformation [T^] and [T^] axe derived in the following subsections. They 
are evaluated at the nominal values of link lengths and angles for calculating [es]. 
Both [To^] and [Tq®] involve Ei, the actual length of link 1 . The mean value oi Ri is 
hfixed + k. However, the angle 7a does not depend upon the length of link 1 and it 
can be taken to be any arbitrary positive value. So the mean value of Ri is taken to 
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Figure 4 . 5 : Angle 7 Made by a line OP With the XZ-plane 

be hfixed: be., k is taken to be zero. Hence [65] is given by 

W = [ 0 0 0 1 ] [T|]^ ( 4 . 15 ) 

where [Tf ] is evaluated at = 0 . The subscript m signifies that [Tf] is being evaluated 
at the mean values of link lengths and angles. 

Hence, from Eq. ( 4 . 13 ), the angle 7 a is given by 

7.=tan-‘ (4.16) 

yeL +^lx 

where 65^, Ssy and 65^ are the x, y and z coordinates of [65] respectively. 

4.2.4 The Angle 7 {, 

76 is the angle made by member 3 with a plane perpendicular to link 1 . A line DS 
is drawn parallel to member 3 and passing through D as shown in Fig. 4 . 4 . So 75 is 

the angle made by DS with the A6.Z5-plane. 

A frame XyYrZj parallel to X2Y2Z2 is attached at point D. The angle 76 does not 
depend upon the length of member 3 and it can be taken any arbitrary positive value. 
So the mean value of R3 is taken to be h fixed) b®.) ^3 is taken to be zero. Hence the 
point S m. X7Y7Z7 frame is given by 

[ st ] = [ hfixed 001 ] 

The point S in X5Y5Z5 frame is given by 

[55] = N [ 2 f] 


( 4 . 17 ) 
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where 

= Rot(X5,Q;2) — >■ Rot(l5,;02) 


Or, 

[r/] = Rot(X5,a2)Rot(y5,^2) 

C0SP2 0 — sin/^ 0 

sinQ!2siii^2 cos 0:2 sin 02 cos 0 ( 4 - 18 ) 

cosQ!2sin/32 —sin 0:2 cos 02 cos 0 
0 0 0 1 
The transformation [Tj ] is evaluated at the mean values of angles 0:2 and ^2- Hence 



[^ 5 ] ^Zfixed 0 0 1 J [-^]^ 

where m signifies that [Tg’] is being evaluated at the mean values of angles 0:2 and p2- 
The angle 7^ is then given by 

76 = tan-^ - (4.20) 

Y 

where 55*, s^y and s^z are the x, y and z coordinates of [55] respectively. 


4.3 The Transformation Matrix [Tq^] 

The frame XiYiZi is obtained by a rotation 0:1 about Xo axis, a rotation about 
axis and a translation (ai, 61, Ci) with reject to the base frame XoVoZo (Fig. 4.6). 
The transformation [T(J] is then given by [103] 

= Rot(J!iro, CKi) — ^ Rot(lo,/3i) Trans(Rflse, fli, 61, Cl) 

Or, 

= Rot(Xo, Q:i)Rot(ro, / 5 i)Trans(Rase, oi , 61, ci) 

Hence 

cos A 0 -sin/3i 0 

sinQ!isin/3i cos 0:1 sinaicos^i 0 

coscKisinA — sinai cosaicosj^i 0 

ai 61 1 

4.4 The Transformation Matrix [T^] 

The frame X^Y^Z^ in Fig. 4.4 have axes parallel to those of J^iF!.^!. Another frame 
XeFeZe is drawn parallel to X^Y^Z^ with origin at E as shown in Fig. 4.7. 
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Figure 4.6: The Transformation [Tq 


The frame X 5 Y 5 Z 5 with respect to the base frame XqYoZq is given by 

[2j] = [r„‘] -4 'ivaiis(yi,;!u) 

Or, 

[To®] = TransiYuRti) [tJ' 

The frame X^YeZe with respect to the base frame XqYqZo is given by 

j^T^j = [Tq] -> Trans(Base,0,0,-i?2i) 

Or, 

[Tq] = [Tq] Trans(5ase,0,0,-i?2i) 

Substituting [T^] from Eq. (4.22) in the above equation yields 

r^] = Trans(yi, i^ii) [?o] Trans{Base,0A -R 2 i) 


(4.22) 


(4.23) 


The frame ^ 212^2 is obtained by two successive rotations 0:2 and about Xe 
and Ye axes respectively. Hence the frame X 2 T 2^2 with respect to XeYeZe is given 
by 


= Rot(X6,Q:2) Rot(ie,/32) 
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Z-stage platen 


Figure 4.7: The Transformation [Tq] 


Or, 


Rot(X6)Qi2) Rot(F6)/?2) 

cos /?2 0 — sin/^ 0 

sin 012 sin /32 cosq !2 sina 2 Cos/^ 0 
cos 0 ( 2 sin ^ —sin 012 cosq! 2 Cos/% 0 
0 0 0 1 
Now the frame X2Y2Z2 with respect to the base frame XqYoZq is given by 

[ro‘] = [2j] ^ [3|] 

Or, 

[r„^J = [Ti] [< 

Substituting [T^] from Eq. (4.23) into [T^] above yields 
[r^] == Trans(yi,i?ii) [ro^] 

Let [Tyi] and [Ts] be the two transformation matrices such that 

[T^] = fr|] Trans(Fi,i?ii) 


(4.24) 
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Therefore, 


[Tb] = [T(JJ Trans{Base,0,0,-R2i) 

cos /?2 0 — sin,^ 0 

siiiQ !2 sm ^2 cos 0:2 sina2Cos/32 0 
coso! 2 sin ^2 —sin 0:2 cos 0 : 2 cos 0 

0 Rii 0 1 


(4.25) 


cos/3i 0 — sin/?i 0 

j j _ ' sin q;i sin Pi cos ai sin ai cos pi 0 

cosaisin/?! -sinai cos 0:1 cos /?i 0 

^^1 61 Cl — i?2l 1 

The transformation [Tq ] is now given by 

[Jj] = [TxlIJil 

where [T^] and [Tb] are given by Eqs. (4.25) and (4.26) respectively. 


(4.26) 


(4.27) 


4.5 The Point Q on the Contour of a Slice Traced 
by the Nozzle Tip 

The nozzle tip traces the contoxrr of a slice. It deposits the material over the previous 
layer of a part as the part is progressively built. When the first layer is being built 
the nozzle tip is buried .01 inch into the foam foundation [113]. So the nozzle tip is 
always in touch with either a previously built layer or the foam foimdation. Therefore 
the coordinates of the nozzle tip are the coordinates of a point on the part being built. 
The coordinates of the nozzle tip Q in the frame is given by 

[<ftl = [iJ3l -fti 1 ] (*-28) 

Therefore, the coordinates of the nozzle tip Q in Xoio^o fi^ame is given by 

[?) = W [lo] (<-29) 

Substituting [Tq] from Eq. (4.27) yields 

The above equation gives the coordinates of a point on the contour of a shoe traced 
by the nozzle tip. Let 

I 9 I = [ I V « 1 ] 


(4.31) 
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4.6 The Means and the Variances of the Depen- 
dent Variables x, y and z 


From Eqs. (4.30) and (4.31) the coordinates of a point on the contour of a slice traced 
by the nozzle tip are given by 


X — X 

y = i/(Vi,V 2 ,...,yi 5 ) (4.32) 


Using the approach given in Chapter 2, the means and the variances of the de- 
pendent variables x, y and z are as follows: 


m[x] = x{m\y^^ i = 1,2, ..., 15) 

m[y\ = y(m[V-], i = 1,2, . . ., 15) 

m[z] = 2 (m [Vi] , i = 1,2, . . ., 15) 



(4.33) 


The partial derivatives are evaluated at the mean values of the A^riables Vi as de- 
noted by the subscript m. The mean values m|Vi] of the variables ^Vi are as given in 
Section 4.2.1. The influence coefficients , (^) (^) ^6®^ the partial 

derivatives of matrices [T^] and [Tb] and are derived in the following subsections. 


4.7 Partial Derivatives of the Elements of [Ta] 


The transformation matrix [T^j is given by Eq. (4.25). Expressing it in terms of the 
random variables V^, it follows that 


[Ta] 


cos Vis 0 —sin Vis 0 

sin Vi 4 sin Vis cosVi 4 sin Vu cos Vis 0 

cosVi 4 sinFis -sinVi4 cos Vi 4 cos Vis 0 

0 Vi 0 1 


Hence the elements of [T^] are functions of Vi, Vi4 and Vis only, i.e.. 


(4.34) 


{Ta]=-Ta{Vi,Vu,Vis) 


(4.35) 
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vari” 


d[TA] 

dVi 




d[TA] 

dVis 

d[TA] 

dVi 


0 0 0 0 
0 0 0 0 
0 0 0 0 
0 10 0 


0 0 0 0 
• sin Vu cos Vu cos V15 0 
cos 1^14 - sin Vi4 cos Vu 0 


cos Vi4 sin Vu 
— sin Vu sin V15 

0 0 

— sinVis 0 — cosVis 0 

sin Vi4 cos Vu 0 — sin V14 sin 1^5 0 
cos Vu cos Vu 0 — cos Vu sin V15 0 
0 0 0 0 

0 i = 2 to 13 


0 0 


( 4 . 36 ) 


4.8 Pcirtial Derivatives of the Elements of [Tb] 

The transformation matrix [Tb] is given by Eq. ( 4 . 26 ). Expressing it in terms of the 
random variables Vi, it follows that 


[Tb] = 


cos Via 
sin Via sin Via 
cos Via sin Via 

ai 


0 — sin Via 0 

cos V\<2 sin Via cos Via 0 
- sin Via cos Via cos Via 0 

bi Cl — Va — V 7 1 


( 4 . 37 ) 


Hence the elements of [Tb] are functions of Va, V, Va and Via only, i.e., 

[Tb]=Tb{V2,V7,Vi2,Vi3) (4-38) 

Therefore, the partial derivatives of the elements of [Tb] with respect to the vari- 
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ables Vi are given by 
d[TB] _ ^ 


dV2 


dlTs] 

dVi2 


d[TB] 

dVn 

dm 

dVi 


dVr 

0 0 0 
0 0 0 
0 0 0 
0 0 - 


0 

0 

0 

1 0 


0 0 
cos Vi 2 sin Via — sin Via 
- sin Vi 2 sin Via - cos Via 
0 

-sinV^ia 0 
sin Via cos Via 0 
cos Via cos V^a 0 
0 0 


0 0 

cos Via cos Via 0 
- sin Via cos Via 0 
0 0 0 
-cos Via 0 
sin Vi 2 sin Fia 0 
cos Via sin Fia 0 
0 0 


(4.39) 


= 0 


i = 1, 3 to 6, 8 to 11, 14 to 15 


4.9 The Influence Coefficients (|^)'and (^)' 

The point Q on the contour of a slice is given by Eq. (4.30), 

[q] = fe] [2a] [2^b] 

where [ga] is given by Eq. (4.28), that is 

[ga] = [ i?3i “-^41 ~2?5 i 1 ] 

Taking partial derivatives of [q] with respect to a variable Vi yields 




dVi 
i = 1 to 19 


(4.40) 


The derivatives of the elements of the matrices [Ta] and [Tb] with respect to the 
variables Vi have been found in the previous subsections. The position vector [ga] can 
be expressed in terms of the variables V as follows: 


tel 


Hence 


V 3 + Ve -V 4 - Vs -Vs -Vg- Vio -Vii 1 ] 
[ga] = g2(F3, Vi,.v. , •.,^ 11 ) 
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Taking partial derivatives of [92] with respect to the variables Vi yields 


dVs 

dV4 

d[q2] 

dVs 

dVe 

d[q2] 

dVs 

d[q2] 

dVg 

d[q2] 

afe] 

dVi 


=[ 1000 ] 
=[0-100] 
= [0 0 - 10 ] 
=[1000] 

= . [ 0 -1 0 0 ] 
=[00-10] 
=[00-10] 
=[00-10] 
= [0 0 0 0 ] 


i =1, 2, 7, 12 to 15 


(4.41) 


The X, y and z coordinates of Eq. (4.40) gives and respectively. Squar- 
ing them gives the influence coefficients (^) ' influence 

coefficients are obtained by substituting [ 92 ]) [Ta]) [Tb] and their partial derivatives 
into Eq. (4.40). 


4.10 Numerical Results 

Numerical results are obtained with the following input values (See Figs. 4.2 and 4.3): 


ai = —0.15 m 

fixed = 0.09 m 

fixed = 0.12 m 

Oi = 1.0° 


= -0.10 m 

hmax = 0.25 m 

hmax = 0.25 m 

A = 0 . 0 ° 


Cl = 0.22 m 

I 2 = 0.07 m 

I 4 = 0.01 m k = 

0(2 = 0 . 0 ° = 


0.15 m 
- 1 . 0 ° 


Absolute Error in positioning member 2 along link 1: ei - 0.0001 m 

Tolerance per unit length on links 2 , 4 and 5; 62 — 64 — £5 — 0.001 m/m 

Absolute error in positioning member 4 along member 3. £3 — O.CXXll m 

t> j. , 1 • • 10. Ci2 = 0.0001 m 

Radial clearance in pair 12 . 

D 1 - 1 T • • QA. C 34 = 0.0001 m 

Radial clearance in pair 34. 

Absolute tolerance in -^2 direction at the head attachment: ejo = 0.0002 m 

Absolute tolerance in Z 2 direction at the nozzle tip attachment, en 
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Absolute tolerances on angles: 

The standard deviation cr of error at a 


ei2 = ei3 = ei4 = ejg = 0.03° 
point on the work surface is given by [110] 


(T 


\/D[x] + D[y] + D[z] 

VDSim, sa; 


where 

DSUM=:D[x] + D[y] + D[z] 

Hence, 

DSUM = 0-2 

i.e., DSUM is the square of the standard deviation of error at a point. 

The variance D[x], D[y] and D[z] are obtained at the nozzle tip Q by varying the 
link lengths k and ^3 in the ranges 0 < < .25 and 0<k< .25 with increments of 

.025 each. Table 4.1 lists the variances D[x], D[y] and D[z], and their sum DSUM for 
these values of li and Z3. 

From the table it is seen that the variances D[x] and D[y] are of the same order 
and they vary slightly over the work surface. The variance D[z] is an order higher. 
This is because the influence coefficients have larger magnitudes compared to 

and ■ Table 4.2 lists the influence coefficients for k = .25 m and k = 

.25 m for the input values mentioned above. The variance D[z] varies significantly 
over the work surface causing DSUM to vary over it. The maximum value of DSUM 
is found at li = .25 m and Z3 = .25 m. 

The three-sigma bands of error in tracing several curves by the nozzle tip is ob- 
tained. The three-sigma band is the band of 3a and — 3cr values of error at a point. 
Figure 4.8 shows the three-sigma band for tracing a curve by the nozzle tip whose 
projection on the JAolo-plane is a line from (.05, .05) to (.20, .20). The three-sigma 
band increases with parameter's along the line. Figure 4.9 shows the three-sigma 
band for tracing a cmrve whose projection on the Xoio^pl^’ne is a circle with center at 
(.15, .1) and radius .05 m. The three-sigma band decreases with angle psi along the 
circle and then increases. From Table 4.1 it is observed that the overall error varies 
appreciably over the work surface. 
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Variances 




and Their Sum 

-"dR- 


Variances and Their Sum 




DSUM 


10 




0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.025 

0.026 

0.025 

0.025 

0.025 

0.025 

0.025 

0.025 

0.025 

0.025 

0.025 

0.050 

0.050 

0.050 

0.050 

0.050 

0.050 

0.050 

0.050 

0.050 

0.050 

0.050 

0.075 

0.075 

0.075 

0.075 

0.075 

0.075 

0.075 

0.075 

0.075 

0.075 

0.075 

0.100 

0.100 

0.100 

0.100 

0.100 

0.100 

0.100 

0.100 

0.100 

0.100 

0.100 

0.125 

0.125 

0.125 

0.125 

0.125 


0.000 

0.025 

0.050 

0.075 

0.100 

0.125 

0.150 

0.175 

0.200 

0.225 

0.250 

0.000 

0.025 

0.050 

0.075 

0.100 

0.125 

0.150 

0.175 

0.200 

0.225 

0.250 

0.000 

0.025 

0.050 

0.075 

0.100 

0.125 

0.150 

0.175 

0.200 

0.225 

0.250 

0.000 

0.025 

0.050 

0.075 

0.100 

0.125 

0.150 

0.175 

0.200 

0.225 

0.250 

0.000 

0.025 

0.050 

0.075 

0.100 

0.125 

0.150 

0.175 

0.200 

0.225 

0.250 

0.000 

0.025 

0.050 

0.075 

0.100 


4.9343 

4.9264 
4.9187 
4.9109 
4.9032 

4.8954 

4.8877 
4.8801 

4.8724 
4.8648 
4.8572 
4.9304 
4.9226 

4.9148 

4.9070 

4.8993 

4.8916 

4.8839 
4.8763 
4.8686 

4.8610 

4.8534 

4.9265 

4.9187 

4.9109 

4.9032 

4.8955 

4.8878 

4.8801 

4.8725 

4.8648 

4.8572 
4.8497 
4.9226 

4.9149 

4.9071 

4.8994 

4.8917 

4.8840 
4.8763 
4.8687 

4.8611 

4.8535 
4.8459 

4.9188 

4.9110 

4.9033 

4.8956 

4.8879 

4.8802 
4.8725 

4.8649 

4.8573 
4.8497 
4.8422 
4.9149 

4.9072 

4.8995 
4.8917 

4.8841 


4.9657 

4.9619 

4.9581 

4.9543 

4.9506 

4.9468 

4.9431 

4.9394 

4.9358 

4.9321 

4.9285 

4.9618 

4.9580 

4.9542 

4.9505 

4.9467 

4.9430 

4.9393 

4.9356 

4.9320 

4.9284 

4.9248 

4.9580 

4.9542 

4.9504 

4.9466 

4.9429 

4.9392 

4.9355 

4.9318 

4.9282 

4.9246 

4.9210 

4.9541 

4.9503 

4.9465 

4.9428 

4.9391 

4.9354 

4.9317 

4.9281 

4.9244 

4.9208 

4.9172 

4.9502 

4.9465 

4.9427 

4.9390 

4.9353 

4.9316 

4.9279 

4.9243 

4.9207 

4.9171 

4.9135 

4.9464 

4.9426 

4.9389 

4.9352 

4.9315 


1.8052 


1.8463 

1.8951 

1.9514 

2.0154 

2.0870 

2.1662 

2.2530 

2.3474 

2.4495 

2.5591 

■ 1.8197 

1.8608 

1.9096 

1.9659 

2.0299 

2.1015 

2.1807 

2.2675 

2.3619 

2.4639 

2,5736 

1.8380 

1.8791 

1.9278 

1.9842 

2.0482 

2.1198 

2.1989 

2.2858 

2.3802 

2.4822 

2.5918 

1.8601 

1.9012 

1.9499 

2.0063 

2.0703 

2.1418 

2.2210 

2,3078 

2.4023 

2.5043 

2.6139 

1.8859 

1.9271 

1.9758 

2.0322 

2.0961 

2.1677 

2.2469 

2.3337 

2.4281 

2.5302 

2.6398 

1.9157 

1.9568 

2.0055 

2.0619 

2.1258 


2.7952 

2.8352 

2.8828 

2.9380 

3.0008 

3.0712 

3.1493 

3.2350 

3.3282 

3.4292 

3.5377 

2.8089 

2.8489 

2.8965 

2.9517 

3.0145 

3.0849 

3.1630 

3.2487 

3.3420 

3.4429 

3.5514 

2.8264 

2.8664 

2,9140 

2.9692 

3.0320 

3.1024 

3.1805 

3.2662 

3.3595 

3,4604 

3.5689 

2.8477 

2.8877 

2.9353 

2.9905 

3.0533 

3.1238 

3,2018 

3.2875 

3.3808 

3.4817 

3.5902 

2.8728 

2.9128 

2.9604 

3.0156 

3.0785 

3.1489 

3.2270 

3.3126 

3.4059 

3.5068 

3.6154 

2.9018 

2.9418 

2.9894 

3.0446 

3.1074 




m 

0.125 

0.125 

0.125 

0.125 

0.125 

0.125 

0.150 

0.150 

0.150 

0.150 

0.150 

0.150 

0.150 

0.150 

0,150 

0.150 

0.175 

0.175 

0.175 

0.175 

0.175 

0.175 

0.175 

0.175 

0.175 

0.175 

0.175 

0.200 

0.200 

0.200 

0.200 

0.200 

0.200 

0.200 

0,200 

0.200 

0.200 

0.200 

0.225 

0.225 

0.225 

0.225 

0.225 

0.225 

0.225 

0.225 

0.225 

0.225 

0.225 

0.250 

0.250 

0.250 

0.250 

0.250 

0.250 

0.250 

0.250 

0.250 

0.250 

0.250 


10“®m- 
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1 — 8 , 




DSUM 


10 - 


0.125 

0.150 

0.175 

0.200 

0.225 

0.250 

0.000 

0.025 

0.050 

0.100 

0.125 

0.150 

0.175 

0,200 

0.225 

0.250 

0.000 

0.025 

0.050 

0.075 

0.100 

0.125 

0.150 

0.175 

0.200 

0.225 

0.250 

0.000 

0.025 

0.050 

0.075 

0.100 

0.125 

0.150 

0.175 

0.200 

0.225 

0.250 

0.000 

0.025 

0.050 

0.075 

0.100 

0.125 

0.150 

0.175 

0.200 

0.225 

0.250 

0.000 

0.025 

0.050 

0.075 

0.100 

0.125 

0.150 

0.175 

0.200 

0.225 

0.250 


4.8764 

4.8688 
4.8611 

4.8536 

4.8460 

4.8384 
4.9111 

4.9034 
4.8956 

4.8803 

4.8726 

4.8650 

4.8574 

4.8498 

4.8423 

4.8347 
4.9073 

4.8996 

4.8919 

4.8842 

4.8765 

4.8689 

4.8613 

4.8537 

4.8461 

4.8385 
4.8310 

4.9035 

4.8958 

4.8881 

4.8804 

4.8727 

4.8651 

4.8575 

4.8499 

4.8424 

4.8348 
4.8273 

4.8997 

4.8920 

4.8843 

4.8766 

4.8690 

4.8614 

4.8538 
4 . 8462 . 
4.8387 
4.8312 
4.8237 

4.8959 

4.8882 

4.8805 
4.8729 
4.8653 
4.8577 
4.8501 

4.8425 
4.8350 
4.8275 
4.8200 


4.9278 

4.9242 

4,9205 

4.9169 

4.9133 

4.9098 

4.9426 

4.9388 

4.9351 

4.9277 

4.9240 

4.9204 

4.9168 

4.9132 

4.9096 

4.9061 

4.9388 

4.9350 

4.9313 

4.9276 

4.9239 

4.9203 

4.9166 

4.9130 

4.9095 

4.9059 

4.9024 

4.9349 

4.9312 

4.9275 

4.9238 

4.9202 

4.9165 

4.9129 

4.9093 

4.9057 

4.9022 

4.8987 

4.9312 

4.9274 

4.9237 

4.9201 

4.9164 

4.9128 

4.9092 

4.9056 

4.9020 

4.8985 

4.8950 

4.9274 

4.9237 

4.9200 

4.9163 

4.9127 

4.9091 

4.9055 

4.9019 

4.8984 

4.8948 

4.8913 


2.1974 

2.2766 

2.3634 

2.4578 

2.5599 

2.6695 

1.9492 

1.9903 

2.0390 

2.1594 

2.2309 

2.3101 

2.3969 

2.4913 

2.5934 

2.7030 

1.9865 

2.0276 

2.0763 

2.1327 

2.1967 

2.2682 

2.3474 

2.4342 

2.5286 

2.6307 

2.7403 

2.0276 

2.0687 

2.1175 

2.1738 

2.2378 

2.3094 

2.3885 

2.4753 

2.5698 

2.6718 

2.7814 

2,0725 

2.1137 

2.1624 

2.2187 

2.2827 

2.3543 

2.4335 

2.5203 

2.6147 

2.7167 

2.8263 

2.1213 

2.1624 

2.2111 

2.2675 

2.3314 

2.4030 

2.4822 

2.5690 

2.^4 

2.7654 

2.8751 


3.1778 

3.2559 

3.3416 

3.4349 

3.5358 

3.6443 

2.9345 

2.9745 

3.0221 

3.1401 

3.2106 

3.2887 

3.3743 

3.4676 

3.5685 

3.6771 

2.9711 

3.0111 

3.0587 

3.1139 

3.1767 

3.2472 

3.3252 

3.4109 

3.5042 

3.6051 

3.7136 

3.0114 

3.0514 

3.0990 

3.1542 

3.2171 

3.2875 

3.3656 

3.4513 

3.5446 

3.6455 

3.7540 

3.0556 

3.0956 

3.1432 

3.1984 

3.2612 

3.3317 

3.4098 

3.4954 

3.5887 

3.6897 

3.7982 

3.1036 

3.1436 

3.1912 

3.2464 

3.3092 

3.3797 

3.4578 

3.5434 

3.6367 

3.7377 

3.8462 


Table 4.1: 
h and k 


Variances and Their Sum at the Nozzle Tip Q in FDM Obtained by Varying 
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The Influence Coefficients 

i 

(fiS)m 



1 

O.OOOe+00 

9.998e-01 

1.745e-02 

2 

O.OOOe+00 

O.OOOe+00 

-l.OOOe+00 

3 

9.998e-01 

-3.046e-04 

1.745e-02 

4 

O.OOOe+00 

-9.998e-01 

-1.745e-02 

5 

1.745e-02 

1.745e-02 

-9.997e-01 

6 

9.998e-01 

-3.046e-04 

1.745e-02 

7 

O.OOOe+00 

O.OOOe+00 

-l.OOOe+00 

8 

.O.OOOe-l-00 

-9.998e-01 

-1.745e-02 

9 

1.745e-02 

1.745e-02 

-9.997e-01 

10 

1.745e-02 

1.745e-02 

-9.997e-01 

11 

1.745e-02 

1.745e-02 

-9.997e-01 

12 

O.OOOe+00 

1.377e-01 

3.325e-01 

13 

-1.377e-01 

O.OOOe+00 

-3.726e-01 

14 

1.745e-04 

1.502e-01 

-7.379e-03 

15 

-1.435e-01 

6.502e-03 

-3.725e-01 


Table 4.2: The Influence Coefficients for li = .25 and k = .25 in FDM 



Figure 4.8: The 3o- band for Tracing a Curve by the Nozzle Tip whose Projection 
ZoFo-Plane is an Inclined Line from (.05, .05) to (.20, .20) 
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Figure 4.9: The 3<r Band for Tracing a Curve by the Nozzle Tip whose Projection in 
XoFo-Plane is a Circle With Center at (.15, .10) and Radius .05 m 




Chapter 5 


ANALYSIS OF MECHANICAL 
ERROR IN SL 


5.1 Geometric Modeling for Description of the SL 
Process 

The StereoLithography (SL) process comprises of a Isiser which emits a laser beam 
(Fig. 1.4). An elevator platform rests below the surface of a photocurable liquid 
acrylate resin in a vat. Initially, the elevator platform is located at a distance from 
the surface equal to the thickness of the first bottom-most layer of the part to be built. 
The laser beam, deflected by the galvanometer-driven mirrors, scribes the contour of 
the slice of a part. It then hatches the interior of the contoin:. The liquid when 
exposed to the beam solidifies. The elevator platform is then lowered by a distance 
equal to the thickness of the next slice. The subsequent layers are produced similarly 
[31, 60, 63, 45]. The SL process is discussed in detail in Section 1 . 1 . 2 . 3 . 

The geometric model of the SL process is as follows. The laser beam emitted 
by the laser strikes the mirror GHIJ and gets deflected towards the resin surface 
(Fig. 5.1). It is assumed that there is a reflection from a single mirror. The ray EP 
is the centerline of the beam incident on the mirror. The ray FQ is the centerline 
of the beam reflected from the mirror (Fig. 5.1). A freime XoYqZq with AoTo-plane 
being horizontal is chosen as the base frame. 

The elevator platform moves vertically up and down. Therefore the bottom of the 
current layer being built is horizontal as is the liquid surface. Let as and 63 be the x 
and y coordinates of the front left end of the elevator platform. Let the surface of the 
liquid is given by Zq — C3. A coordinate frame XiYiZi is attached at Oi (as, 63 , 03 ) 
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with X,K,-plan,. horizontal. Hence the -plane is the free surface of the liquid 
The cortf^sponciinK axes of XoKoZo and XiViZ^ are parallel. 



figure 5 . 1 : Geometric Model of the SL Process 


The ray EP emanating from E(oi,6i,Ci) strikes the mirror at F and is reflected 
towards the point Q on the XiKrplane. Ideally the ray EF is along Xi direction. 
However, due to mechanical erroiB it is oblique in space and its direction is given by 
two angles & and ^ aa shown in Fi^. 5.1 and 5 . 3 . Nominal values of 0 and (f) axe taken 
such that the ray EP is along Xq direction. 

The corner point G(rt2,i!>2,C2) of the mirror is taken to be fixed in space. The 
mirror rotates in spatx^ such that the point G is fixed. A firame X2Y2Z2 is chosen 
with origin at G and the XgFa-plane being the mirror plane. The frame X2Y2Z2 is 
obtained front the frame XqYqZo by a rotation a about Xo-axis, a rotation (3 about 
ko-axia and a translation (02,621^2) in the base frame XqYoZq. The angles a and /3 
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ar6 obt^uiKHl stK.h tJifit t»lK' riiy EF is r6fl0ct6d towards the desired point Q on the 
work surfaces XiYi* 

5.2 Stochastic Model of the SL Process 

Chapter 2 ciiscuswis stochastic modeling of an RP process. Following the approach 
given there, the stochiistic modeling of the SL process is done in the following. 

Due to mechanical errors, the locations of the source E(ai,6i,ci) and the corner 
point G{a 2 , ^ 2 , ^‘ 2 ) ^^^e in error {Fig. 5.2). The ray EF is not along Xo direction due to 
mechanical tirrors and t.hfire are tolerances on the angles 9 and 4>. When oriented to 
achieve the (ie.sired refhietion, the mirror is not exactly rotated by the angles a and 
/5 but the angh^s et and are achieved with some tolerances. 



Y 
* • 


X| ' 

FigiiWi 5.2: Geometric Model of the SL Process Showing Error in Dimensions 
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T\w plan.' is tlH‘ resin kivcl and is given by Zq = cg in the base frame. The 
liquid rt'siu shrinks upon polymerization. Therefore the resin level decreases after 
laser drawinj'; of f'adi layer. A sensors checks the resin level upon completion of laser 
drawiiif; of eaeh layei'. In the (!vent f.liat the sensor detects a resin level that is not 
within lh<' foleranee limits, resin is displaced in the vat using a plunger. When the 
resin kn'el lies wit hin t he tol(;ranc(; limits, the levelling operation is completed [60, pp. 
13|. Hence rg, the /^o-coordinate of the point G, is in error and there is a tolerance 
in it. Th<^ t.olerane(' in eg t.akCvS into account the error in positioning the platform in 
l?odirec.t ion an<l the error in the t.hickness of the current layer swept by the recoater 
blade also. 

So, in all, t-here are <iimenfiiona in error. Therefore there are eleven random 

varialdfus K, i I !<» H, involved in the process corresponding to these eleven dimen- 
sions. I'liese random variahkis are independent of each other. The random variables 
are giv(ni as follows: 

Vl fp 1 « 1 * ^2 1 ^2i Kl “ Cl + €3, V4 = ^ + €4, Vs = + 65, 

V(j (12 i V7 " 1 C?) 14 = Cg + €3, = O' + 69, Vxo = P + ClO) 

^ 1 Cg f » n I 

(5.1) 

where („ I 1 ,2,.'l,d,7,M and 1 1 we absolute tolerances in meters on respective coor- 
dinates. i 4,5,11 and 10, are absolute tolerances in radians on angles 9, a and 
P respectively. 

'flu' rnearj vahws of the wiriables V*, i = 1 to 11, are given by 

mjKil *'*»! = h, mfVs) - cj, m[V4] = 9, mfVs] = 

m|V«j aa, mfVyl = mjVs) = ca, mfVg] = a, m[yio] = P, (5-2) 

rnjl/al .. eg, 

The variance /^(Idj of the variable V< is given by 

“ (1)^ i = 1 to 11 (5-3) 


5.3 The Incident Ray EF 

Gemsider the ray BP in Fig. 5. 1 which emanates from the source E and strikes the mir- 
ror GIH.I at, F. f/rt \tf\ he a unit vector along EF in the base frame. The orientation 
of EF i.M given by two angles 0 and (p shown in Fig. 5.3. Due to mechanical errors 
the ray EF is oblique in space and it intersects the Xo>^o-plane at say, K. Therefore 
the unit vector |f/| is given by 
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lf\ { sin sin 0 ~cos0 0 


(5.4) 


I I i# 


0 


0 

- ^ / 


/ 



\'^ ' /I 

\ /EK sinQsinct) 

n/ 

f K ',iii0(O', <{> 


Figiirf^ 5.3: The vector EF 
The ray EF is then {^iven hy (38j 

(r(01 [ej ) f,(g/) 0<t<oo 


where 

(c| = [ ai bi Cl 1 


(5.5) 

(5.6) 


5.4 The Transformation Matrices [Tq] and [T 2 ] 


The friiine in obtained by a rotation a about Xo axis, a rotation /? about Yq 

axis an<i a trauHlnticm ( 03 , 63 , 03 ) in the ba* frame XqTo^o (Fig. 5.4). The transfor- 
mation j'/’lj is Kiwn iiy [103) 

[7;?] s Eot{Xo, a) -4 Rot(ro, p) 'B:ans(Bose, 02, 62, ca) 


Or, 


This yields 


rj] ^ Rot(Xo,«)Rot(Fo,i5)Tran8(Base, 02 , 62 , 02 ) 


€080 
sin Of sin ^ 
mn a sin 0 


0 —8X0.0 0 

caosof sinacos/? 0 
-sina co8aco80 0 
62 C2 1 


(5.7) 
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Figure 5 . 4 : The Transformations [T^] and [T^ 


The transformation (Jf), the inverse of [Tq] is given by 


cos/? 

sin a sin P 

cos a sin /? 

0 ' 

0 

COS a 

~ sin a 

0 

— sinj^ 

sin p cos 

cos a cos /? 

0 

-a2<»8^ 
-f cg sin /? 

“■02 sin 0! sin /? — 62 cos a 
—C2 sin a cos/? 

-02 cos a sin /? + 62 sin a 
—C2 cos a cos (5 

1 


5.5 Solving for the Mirror Plane and the Intersec- 
tion Point F 

This section finds the nornaal to the mirror plane such that the ray EF from the 
given source E is reflected towards the point Q on the resin surface. The mirror plane 
passes through the fixed point G. The normal gives the orientation of the plane of 
mirror in space. This section also finds the point of intersection F of the incident ray 
and the mirror plane. 

The point E as in Eq. ( 5 . 6 ) is given by 

[e] = [ ai h Cl 1 

All the position vectors including [e] above and other derivations in this section 
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0 7r/2.0, <^ = 0.0 . (5.11) 

At th<‘S(‘ noiiiinal valiH’s flu- pknc of rcfloction BFQ, as shown in Fig. 5.5, has 
a norinnl vector with zero .Xo-ftwrciinatc;. That is, the plane of reflection is normal 
to the plnne. 'rhis i.s n dcgom-rato case and should be taken care of separately, 
rntlu'r than soiviiif; for t he mirror plane with arbitrary values, 0 and <p. This facihtates 
ol>taininf, t he solnf ion, Iz‘t he the Xo-coordinate of point F. Since ray EF is along 
X() direct ion, therefore, the point F is given by 

I/I = [ A l>i c. 1 ] (5.12) 

U'i the tmit vector FNi, normal to the mirror plane GHIJ at F, is given by 

|/nl| [ zi /? C 0 ] (5.13) 

The mirror plain* pamrs through G(a2, 62,03) and FNj is a unit normal vector to 
it. A point l*(x,y,z.) lies on the mirror plane GHIJ if and only if GP is perpendicular 
to PN|; that is, if ami only if 

PNi,GP = 0 (5.14) 

This gives ilie mirror plane as 

.4(* - tta) f - 62) 4- C(a: - C2) = 0 (5.15) 

Since the point F lies on the mirror plane ^ven by Bq. (5.15), therefore, 

™ **2) .^{^1 “■ ^) "h C{ci — C2) = 0 (5.16) 

Moreover, j/ni] is a unit normal vector. Therefore, 

^34.534.(72 = 1 (5.17) 

Let the plane of rtdieetion, which passts through B, F and Q, is given by 


Ls + My + Nz ~ I 


(5.18) 
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Figure 5.5: The Mirror Plane and the Plane of Reflection 

Since the points E, Q and F lie on the plane given by Eq. (5.18), therefore, it 
follows that 

Lax + Mbx + Ncx = 1 (5.19) 

Lq, + Mqy + Nq. = l (5-20) 

and 

Lfa: + Mbx + JVci = 1 (5-21) 

Equations (5.19), (5.20), (5.21) constitute a set of three equations in four un- 
knowns L, M, N and /a,. Subtracting Eq. (5.21) from Eq. (5.19) yields 

L = 0 (5-22) 

So the above set reduces to the following two equations in two unknowns M and 
N: 


Mbx + Ncx 


1 


(5.23) 
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Mqy + Nq^ = 1 

Solving Eqs. (5.23) and (5.24) yields 

hqz - ciqy 

and 

AT = ~ 

hqz - Ciqy 

The plane of reflection EFQ is now given by 


(5.24) 


(5.25) 


(5.26) 


My + Nz = l (5.27) 

The coefficients A, B and C of the mirror plane and the coordinate of the point 
F are found in the following. 

First Law of Reflection: 

Prom the first law of reflection, the normal vector FNi to the mirror plane, and 
hence the point Ni, lies on the plane of reflection EFQ given by Eq. (5.27). The point 
Ni is given by 

N=[/.+^ h+B Ci-iC l] (5.28) 

Since Ni lies on the plane (5.27), therefore 


M{bi + B) + N{ci + C) = 1 

Subtracting Eq. (5.19) from Bq. (5.29) and substituting L = 0 yields 

MB +NC=0 


where 

Qz 

Substituting Eq. (5.30) into Eq. (5.17) yields 

Or, 


(5.29) 


(5.30) 

(5.31) 


(5.32) 




(5.33) 
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Second Law of Reflection: 

prom the second, law of reflection, the angle of incidence is equal to the angle 
of reflection 72 (Fig. 5.5). That is 

7i =72 

Hence, 

COS71 = cos 72 

FE.FNi _ FNj.FQ 
|FE||FNi| “ |FNil|FQ| 

This yields 

-A = ~ ~ ^1) + ~ 

\/ (<Ix - /*)^ + (qv - h)^ + {q^ - ci)2 

Squaring both sides of the above equation gives 

42 — ~ fx) + ^{Qy ~ ^1) + C{qz — ci)}^ 

{9® - fx^ + {qy - h)^ + {qz - ci)2 

Consider the following equation, 

A = ~ ~ ~ 

y(g* - /x)^ + (9y - hf + («z - ci)2 

The solutions (yli, Si,^!) and (yl2,52,C'2) satisfying Eqs. (5.35) and (5.37), re- 
spectively, satisfy Eq. (5.36) also. Moreover, if (Ai,Bi,Ci) and {A2, 82,02) satisfy 
Eqs. (5,35) and (5.37), respectively, then (-/li, -Bi, -Ci) and {-A2, —B2, -C2) also 
satisfy Eqs. (5.35) and (5.37) respectively. Equation (5.35) as well as Eq. (5.16) are 
homogeneous in A, B and C. It follows that there are four values of the unit nor- 
mal vector [/nj] = ( A B C 0 ] possible on solving Eqs. (5.16), (5.30), (5.32) and 
(5.36). 

Squaring Eq. (5.35) to obtain Eq. (5.36) is equivalent to finding the unit normal 
vector FNi such that 

r(±FE).FNi]'' _ r FNi.(±FQ) l^ 

_ |FE1 IFNil J [ |FNi| |FQ| _ 

For the given points F, E and Q there are four combinations possible; (FE, FQ), 
(FE, -FQ), (-FE, FQ) and (-FE, -FQ). So there are four solutions to unit normal 
vector FNj possible. Let the vectors FEi and FQi be in the direction opposite to 
PE and FQ respf.ctively. Then the four solutions to unit normal vector are FNia, 
FNib, FNi, and FNid as shown in Fig. 5.6. FNia and FNxc are bisectors of angles 
QFE and EFQi respectively. Therefore FNi* and FNie are perpendicular to each 


(5.34) 

(5.35) 

(5.36) 

(5.37) 
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Figure 5.6; Four Possible Solutions to Unit Normal Vector on the Plane of Reflection 

other. Similarly FNib and FNia are perpendicular to each other. Also the vectors 
FNia and FNic are collinear with vectors FNib and FNm respectively. 

The vector FNja is the desired solution. This is because the incident ray is to be 
reflected downwards, towards the resin surface which is on the AiVx-plane. For such 
a reflection, the unit normal to the mirror plane is given by 

[fnl] = [fnla] = [ A B C 0 ] 

where 

C < 0 (5 38) 

A < 0 

As mentioned earlier, all the derivations and position vectors are in the base fra 
XoYoZq. The coordinate B is given by Eq. (5.30). It follows that 

for Qy > h 
for Qy < hi 


B > 0 
< 0 


(5.39) 
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So there are two equations, Eqs. (5.16) and (5.36) in two unknowns C and as 
follows: 


Mfx - 02) + B{h - bi) + C{Ci - C2) = 0 

{(q* - fx)^ + {qy - bi)^ + {q^ - ci)^} - {A{qx - fx) 

+ B{qy-bi) + C{qx-ci)}^ = 0 

where B is given by Eq. (5.30), i.e.. 


(5.40) 


B = kC 


Since yl < 0, therefore, from Eq. (5.32) it follows that 

A = -^1 - (1 + A:2)C2 (5.41) 


Since A is a real number, therefore, Eq. (5.41) gives the upper limit of \C\. Moreover, 
using Eq. (5.38) it follows that 


1 

vTT^ 


< ^<0 


(5.42) 


The set of equations, Eq. (5.40) above are non-linear in C and fx- Newton-Raphson 
method, an iterative technique, is used to solve them [97]. Let Fi and F 2 represent 
the left hand sides of Eqs. (5.16) and (5.36) above, i.e.. 


F2 


A{fx ~ 02 ) + B{bx — 62 ) + C{cx C 2 ) 

(5.43) 

{fe - + (iv - hf + (S. - Cl)'} 


- {A(q^ - U) + - 6,) + C(q. - Cl)}' 

(5.44) 

of B in Eq. (5.30) with respect to C gives 



(5.45) 

dc 


of A in Eq. (5.32) with respect to C gives 


M = -(1 + k'‘)CIA 

(5.46) 


The Jacobian matrix [J| of the set of equation [F] 


[J] 


Fx 

F 2 


is given by 


- dFx 

m: 

dC 

dfx 

dF2 

dF2 

. dC 

dU - 


(5.47) 
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The partial derivatives of Fi and are as follows: 


dFi 

dC 

m 

dfx 

dF2 

dC 


dFi 

dfx 


- ifx - + k{bi - 62 ) + (ci - C 2 ) ^5 

(5.49) 

= -2(1 + e)C {(9, - + (q, - 6.)^ + (g, - cO^} 

-2 {A{qx — fx) + B{qy - bi) + C{q^ — ci)} 

+ H% - h) + {qz - Ci)| (5 50j 

= -2A^{qx - fx) + 2A {A{qx - fx) + B{qy - bi) + C{qx - ci)} (5.51) 


The initial value of C for Newton-Raphson iteration is taken such that it lies in the 
range given by Eq. ( 5 . 42 ). 


5.6 Orientation of the Frame X2Y2Z2 Using Direc- 
tion Cosines of X2 and Z2 axes 

The frame X2Y2Z2 is obtained by a rotation a about Xo-axis, a rotation jd about 
Fo-axis and a translation {a2,b2,C2) in the base frame XqYqZq successively. For these 
two successive rotations, the X2-axis lies on a plane parallel to J^O'^o-plane for any 
values of a and /?. So once the mirror plane is obtained, the X2-axis is parallel to the 
line of intersection of the mirror plane and the Xo^o-plane. In Fig. 5 . 7 , RST is the 
mirror plane and the line RS is the line of intersection. 

The point F lies on the j!(r2k2-plane. Let the point F in the frame X2Y2Z2 be given 
by 

(/2|=[/2. 0 1] ( 5 . 52 ) 

The position vector [/2] can be obtained from [f] as follows: 

1/2) = (/] [l?] ( 5 - 53 ) 

where [T^] is as given in Eq. ( 5 . 9 ). 

The directions of X2, Y2 and Z2 axes axe chosen such that both fix and fiy are 
positive. This ensures that the ray physically falls on the mirror (Fig. 5 . 8 ). The 
angles a and for these X2, Y2 and Z2 axes are obtained. 

The mirror plane given by Eq. ( 5 . 15 ) can be rewritten as 


Ax + By A Cz = D 


( 5 . 54 ) 
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where D is a constant given by 

D — A.Q2 "h Bb2 Cc2 (5-55) 


The coefficients A, B and C have already been obtained in the previous sections. 
Let the line of intersection RS makes an angle $ from Xo-axis as shown in Fig. 5.7. 
For the plane given by Eq. (5.54), the angle ^ is given by 


tan^ = 


dz 

dx 


Or, 

tan« = (5-56) 

Since the coefficients A and C are negative for any point on the resin surface, 
therefore, the slope of line RS 

tan^<0 (5-^'^) 

The choice of X 2 and Y 2 axes depends upon the coordinates of the points E and 
G. The points E and G are chosen such that Qy = fy> 9y — fz < 9z- Then for 

F to lie in the first quadrant of X 2 Y 2 plane, X 2 -axis should be directed downwards 
and l 2 -axis should lie away from the line RS, as shown in Fig. 5.8. The directions of 
X 2 and F 2 are diosen in this manner. 
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Figure 5.8: Directions of and Y 2 Axes 


Since reflection is taking place in the first octant of the frame XoYqZq, therefore 
there are points on the miror lying in the first octant. Moreover, since both the 
coefficients A and C are negative, therefore, the mirror plane has positive Xo and Zq 
intercepts as shown in Fig. 5.8. Therefore from Eq. (5.57) it follows that 


90 '’ < ^ < 180 " 

Therefore, 

and 

• . ^ 

+ (J2 

Hence the direction cosines of X 2 -axis are given by 


(5.58) 

(5.59) 

(5.60) 


[® 2 ) 


C 


(5.61) 


For the chosen X 2 and Y 2 axes, the Z 2 ~ 8 xis is opposite to [fnl]. Hence the direction 
cosines of .^ 2 - 0 x 18 are given by 


[^2] 


■B 


-c 0 ] 


(5.62) 
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Consider the transformation [r|] given in Eq. (5.7). The first, second and third 
rows of [To^] give direction cosines of X 2 , Y 2 and Z 2 axes in the base frame XqYoZo. 
The fourth row of [Tq^] gives the coordinates of origin G of X 2 Y 2 Z 2 frame. 
Comparing Eq. (5.61) and the first row of [Tq] gives 


sin/? = 


A • 


and 


cos/3 


C 


+ (72 

Since both A and C are negative, therefore, sin/3 and cos/3 both are positive. 
Hence (3 lies in the first quadrant and is given by 


(5 — sin 


-1 


+ C2 

Comparing Eqs. (5.62) and the third row of [Tq] gives 


(5.63) 


sin a = B 


and 


cos a sin /3 = —A 


Substituting sin (3 in the above gives 

cosQ! = \l (A 2 + (72) 

Therefore cos a is always positive. And using Eq. (5.39) it follows that sin a is 
positive for qy > bi and negative for Qy <bi. Hence the angle a is given by 

a = sin“^ B (5.64) 

Prom Eq. (5.53) the coordinates /21 and f 2 y are given by 

f2a> = (A-02)cos/3+(c2-Ci)sin/3 (5.65) 

f 2 y = ( 61 - 52 ) sec a 

With both A and C negative the mirror plane is oriented such that fx > 02 - 
Therefore from Eq. (5.65) it follows that f 2 x is positive. Ftom Eq. (5.66) it follows 
that f 2 y is positive. So / 2 ® and f 2 y are always positive as desired. 
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5.7 Expression for f/, the Value of the Parameter 
t at the Point F 

The ray EF in the frame X 2 Y 2 Z 2 is given by 

h(<)] = [r(t)] [Tj®] (5.67) 

where [r(t)] is given by Eq. (5.5). Let the coordinates of [r 2 (i)] be given by 

{r2{t)] = [ r2x{t) r2y{t) r2z{t) 1 

The ray EF intersects the plane of mirror GHIJ in the X 2 F 2 -plane. Therefore the 
point of intersection is given by 

r2z(t) = 0 (5.68) 

Let [t02_3] be the third column of [T^]. That 


[i02_3] 


' IS 

cos a sin /? 

- sina 
cos Oi cos (5 

-02 cos Q! sin /5 + 62 sin a — C 2 cos a cos /? 


(5.69) 


Therefore Eq. (5.68) yields 

[e] + tf[ef]][t02^]=0 (5.70) 

where tf is the value of the parameter t at the point F. On solving it follows that 


- [e/][t02.3] 


(5.71) 


Let the position vectors [T' 2 (t)] and [r(t)] at t = tf be given by [r2Jf] and [rJ/] 
respectively. Therefore, from Eq, (5.5), 


[rJf] = [e] + tf [ef] 


and from Eq. (5.67), 


[r2Jf] = [r-tf] [r® ' 
The point F in the frame X 2 Y 2 Z 2 is given by 

[f 2 ] = [r2-tf] 


(5.72) 


(5.73) 


(5.74) 
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Figure 5.9: The Reflected Ray FQ 

5.8 The Reflected Ray FQ 


The ray EF is reflected from the mirror as the ray FQ (Fig. 5.9). 
on the ray FQ such that |EF| = |FP|. 

The point E in the frame is given by 

Let P be a point 

h\ = (el [I?] 

Let the coordinates of [ 62 ] be given by 

(5.75) 

[62] = [ 62 ® e2y 622 1 ] 

(5.76) 

The point F in the frame X 2 Y 2 Z 2 , [/z] is given by Eq. (5.52). 
FP in the frame X 2 Y 2 Z 2 is given by 

Hence the vector 

[.^ 2 ] = [ / 2 ®-e 2 x f2y-e2y +^ 22 : 0 ] 

The reflected ray FQ is, therefore, given by 

(5.77) 

K] = [/ 2 ] + sM 0 <S<sq 

(5.78) 

where sq is the value of the parameter s at the point Q. The ray FQ in the frame 
is given by 

[mi] = [m 2 ] [t?] 

where 

(5.79) 

11 

(5.80) 
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The transformation [T^] is given by Eq. (5.7). [To^] is given by 


frM 

[■^oJ 


10 0 0 

0 10 0 

0 0 10 

^>3 C3 1 


(5.81) 


Therefore [Tf], the inverse of [To^], is given by 


[rf] 


1 

0 

0 

— 03 


0 

1 

0 

-bz 


0 

0 

1 


0 

0 

0 


-C3 1 


(5.82) 


Substituting Eqs. (5.7) and (5.82) into Eq. (5.80), it follows that 


[Z’i 


cos/3 0 —sin/3 0 
sin a sin (3 cos a. sin a cos /3 0 

cos a sin /3 —sin a cosocosyd 0 
a% — az ^2 ~ bz C2 — C3 1 


(5.83) 


Let the coordinates of [uq] in Eq. (5.79) be given by 


H] 


Wl® VJly Wu 


1 ] 


(5.84) 


5.9 The Point Q on the Contour of a slice Drawn 
by the Laser on the Resin Surface 

The point Q drawn by the laser on the resin surface is the intersection of the reflected 
ray FQ with the resin surface, i.e., the XiYi-plane. Rrona Eqs. (5.79) and (5.83), the 
point Q in the frame XiY\Zi is given by 


wu = 0 

Let [t20| be the third column of [Tf]. That is 

-sin/3 
sin a cos /3 
cos Q! cos ^3 
C2 - C3 

Therefore Rq. (5.85) yields 


(5.85) 


{i21_3] 


(5.86) 


(/a] + 


[t21-3]=0 


(5.87) 



nhapter 5. Analysis of Mechanical Error in ST. 


113 


whGrG sQ is the of the parameter at the point Q. It follows that 


[/2][t21_3] 

[^2][t21_3] 


(5.88) 


Substituting the value of sq in Eq. (5.79) gives [wi] at sq . Let [w^] at sq be given 
by [wl-sq]. Let tulxsq and wly^sq be the x and y coordinates of [u;l_sg]. Therefore, 


[ml_.sq] = 

wlx-sq wly^q 0 1 

(5.89) 

Let the coordinates of the point Q in the frame XiYiZi be given by 

fe] = 

qiy qiz 1 ] 

(5.90) 

Therefore, 

qix = wlx^q 
qiy = wly^q 
qu = 0 

The point Q in the frame XqYqZq is given by 

(5.91) 

II 

(5.92) 

This yields 

= [ wlx^q+a^ wlysq+bz cz 1 
Let the coordinates of point [q] be given by 

(5.93) 

l9l = 

= [ a; y z 1 

(5.94) 


5.10 The Means and the Variances of the Depen- 
dent Variables x, y and z 

FVom Eqs. (5.93) and (5.94) the coordinates of a point on the contour of a shce drawn 
by the laser beam are given by 

X = 

y = 

z = 


xiVi,V2,...,Vn) 
y(VuV2,...,Vn) 
l^ii) 
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Using the approach given in Chapter 2, the means and the variances of the de- 
pendent variables x, y and z are as follows: 


m[x] = x{m[Vi\, i= 1,2, 11) 

m\y\ = i = 1,2, . . 11) 

m{z] = z{m[V^, i = 1,2, 11) 



m . g (*)>Ki 

The partial derivatives are evaluated at the mean values of the variables Vi as 
denoted by the subscript m. The mean values m[T^] of the variables Vi are as given 
in Sectta6.2. Tte necessary expressions for the influence coefficients 
and can be derived as in the following sections. 


2 2 

5.11 The influence coefficients and 


The point Q on the contour of a slice drawn on the resin surface is given by Eq. (5.93). 
It gives [q] in the base frame XqYqZo in terms of the nominal values of the variables. 
Writing [q] in terms of the variables Vi yields 

(5.97) 


[ql = [ wlx^q+az wly^q+h Vn 1 
Therefore the partial derivatives of [q] with respect to the variables Vi are given 


by 


A 

dVi 

d 

dVi 


11 


[g] 

[g] 


dwlx^q dvly-s^ ^ p 
dVi dVi 

dwlx^q dwl y-sq ^ ^ 
dVn dVii 


i = 1 to 10 


(5.98) 


5.11.1 Partial Derivatives of [wl_sq] 

The position vector [mi] at sq on the resin surface is given by Eq. (5.89). Therefore, 


d 

Wi 


[ml-sg] = 


dwlx^q dwlysq ^ ^ 

dVi dVi 


Taking partial derivative of Eq. (5.79) gives 
d 


dVi 


[wljsq] 


dVi dVY^^ I ^ dVi \ 


+ 


[/2] + ^g [fp^] 


A 

J dVi 


(5.99) 


[T^ 


i = 1 to 11 
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5.11.2 Partial Derivatives of sq 

The value of the parameter s at the point Q, sq is given by Eq. (5.87). It can be 
rewritten as 

sq{fp2] [t21-3] = -[/s] [f21_3] ■ 

Therefore the partial derivative of sq with respect to the variables Vi are given 
by 


dsq 

dVi 


[#2][t21_3| 


film..] . 1/.,^ , 


dVi 
i = 1 to 11 


dVi 

(5.100) 


5.11.3 Partial Derivatives of [/p2], [/ 2 ], [r2_^, [rJfi and tf 

The vector [fp2] given by Eq. (5.77) has its partial derivatives as follows: 


d 

dVi 


[m] = 


df2x de2x df2y_d^ 
dVi dVi dVi dVi dVi 


i = 1 to 11 (5.101) 


Taking partial derivatives of Eqs. (5.74), (5.73) (5.72) and (5.70) it follows that 


A 

dVi 

dVi 
_d_ 

dVi 


[/ 2 ] 
[r2-tf] 
[r-tf] 


\ df2x df2y 


L dVi dVi 

[T2°] + M/l 


0 0 


d 

dVi 


9[r-tf](rr0^ , 


dVi 
d[e] dtf 


dVi 

dW] 


dVi dVi dVi 


[r2_t/] i = 1 to 11 (5.102) 

i = 1 to 11 (5.103) 

i = 1 to 11 (5.104) 


dtf 

dVi 


1 


d[e\ 


a[f02_3] , d[ef] 


[e/l[t02.3] [ 




[t02-3] 


dVi 


i = 1 to 11 (5.105) 


5.11.4 Partial Derivatives of [ef\, [ej and [e] 

Equation (5.4) gives [ef] in terms of the nominal values of the variables. Writing [ef] 


in terms of the variables Vi yields 


[ef] 


sin 14 cos 14 sin 14 sin 14 — cos 14 0 ] 


(5.106) 


Therefore the partial derivatives of [ef] with respect to the variables are given 



Chapter 5. Analysis of Mechanical Ermr in SL 


116 


= [ cos Vi cos ^5 cos Vi sin 1/5 sin Vi o] 

= [ — sinVisinVs sin Vi cos Vi 0 o] (5.107) 

= [ 0000 ] i = 1 to 3 , 6 to 11 

Taking the partial derivative of Eq. (5.75) with respect to the variables Vi gives 

i = ltoll (5.108) 

Equation (5.6) gives the point [e] in terms of the nominal values of the variables. 
Writing [e] in terms of the variables Vi yields 




d 

ev. 

dVi 


W] 

W] 

[ef] 


[e] = Vi Vi Vi 1 


(5.109) 


Therefore the partial derivatives of [e] with respect to the variables Vi are given 


by 


d 

avi 

a 

avi 

a 

avi' 


[e] = [ 1 0 0 0 ] 

[e] = [ 0 1 0 0 ] 

[e] = [ 0 0 1 0 ] 

0000] 


(5.110) 


i = 4 to 11 


5.11.5 Partial Derivatives of the Elements of [T^] and [t21_3] 

Equation (5.83) gives the transformation [If] in terms of the nominal values of the 
variables. Writing [Tf ] in terms of the variables Vi yields 


T2 


cosVio 0 -sinVio 0 

sinVgsinVio cos Vi sinVgCosVio 0 
cosVgsinVit) -sinVg cosVicosVio 0 
Ve — 03 Vi — 63 Is — 1^11 b 


(5.111) 


Therefore the partial derivatives of [Tf ] with respect to the variables Vi are given 
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by 


d 

dVi 


[Tn 


_^rr2 
dVe 


d r^2 


dVr 


m = 


d rr^2 


dVs 


m 


9 Tr^2 


dVg 


d 

dVxo 


[Tn 


[Th 


d 

dVi 


11 


■[T^] 


0 


i = 1 to 5 


[Tt] = 


0 0 0 0 
0 0 0 0 
0 0 0 0 
10 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 10 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 10 


(5.112) 


0 


cos Vg sin Fio 
- sin Vg sin Vio 
0 

— sin Vio 0 
sin Vg cos Vio 0 
cos Vg cos Vio 0 
0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0-10 


0 

-sinl^ 
cosV^ 

0 

- cos ViQ 0 

— sin Vg sin Vio 0 

- cos Vg sin 1^10 0 

0 0 


0 0 

cos Vg cos Vio 0 
sin Vg cos Vio 0 
0 0 


The partial derivative of [t21_3] with respect to a variable Vi is the third column 
of the partial derivative of [T^^] with respect to that variable Vi, i.e., 

d 


d 

dVi 


[t21_3] = third column of 


i = 1 to 11 


(5.113) 
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5.11.6 Partial Derivatives of the Elements of [T^] and [^02_3] 


Equation (5.8) gives the transformation in terms of the nominal values of the 
variables. The transformation [Tf ] in terms of the variables Vi is given by 


T? 


COS Vio 

sin Vg sin Vio 

cos Vg sin Vio 

0 ' 

0 

cosVg 

— sin Vg 

0 

— sin Vio 

sin Vg cos Vio 

cos Vg cos Vio 

0 

— Ve cos Vio 

— Vg sin Vg sin Vio — V? cos Vg 

— V^ cos Vg sin Vio + V 7 sin Vg 

1 

-f sin Vio 

— Vg sin Vg cos Vio 

-V^ cos cos Vio 
/»»■ -* 

-t 4 \ 


(5.114) 

Therefore the partial derivatives of the elements of [T§] with respect to the vari- 
ables Vi is given by 

d 


dVi 


d 

dVe 


d 

dV 


[T^] 


[r«] 


0 




d 

dVs 


[T^] 


d 

dVg 


[T^] 


d 


dVio 




0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


i = 1 to 5, 11 

0 0 0 

0 0 0 

0 0 0 

cosVio — sinVgsinVio -cosVgsinVio 0 

0 0 0 

0 0 0 

0 0 0 

■ cos Vq sin Vg 0 

0 “ 0 
0 0 

0 0 


0 

0 

0 


sinVio -sinHcosVio -cosFgCosVio 0 


0 


0 


0 


0 


cos Vg sin V\Q 


■ sin 


cos Vg COS Vio 


-Vg cos Vg sin Vio + I 7 sin Vg 


■ sin H sin lio 


— cos 

■ sin Vg cos Vio 


Vg sin Vg sin Vio + V 7 cos Vg 
-(-Vg sin Vg cos Vio 


. 1 

— sin Vio 

sin Vg cos Vio 

J — 

cos Vg COS Vio 

0 ' 

0 

0 

0 

0 

— cos Vio 

— sin Vg sin Vio 

— cos Vg sin Vio 

0 

Ve sin Vio 
-l-Vg cos Vio 

— V^ sin V^ cos Vio 
-t-V^ sin Vg sin Vio 

— Ve cos V^ cos Vio 
cos H sin Vio i 

0 


(5.115) 


The partial derivative of 1*02^1 with respect to a variable V, is the third colunm 
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fi 

— [t02_3] = third column of -577 [To] 
dVi dVi 


of the partial derivatives of the elements of [T°] with respect to that variable V^, i.e, 

i = 1 to 11 (5.116) 

f, 

5.12 Numerical Results 

Numerical results are obtained with the following input values (see Figs. 5.1 and 5.2): 

e = 90.0° 4> = 0.0° 


ai 

= 0.0 m 

61 

= 0.15 m 

Cl 

= 1.0 m 

02 

= 0.25 m 

62 

= 0.05 m 

C2 

= 1.02 m 

03 

= 0.0 m 

bz 

= 0.0 m 

C3 

= 0.5 m 


Length of the platform: 

Width of the platform: 

Absolute tolerances on x-coordinates: 
Absolute tolerances on y-coordinates: 
Absolute tolerances on z-coordinates: €3 
Absolute tolerances on angles: €4 = €5 


0.30 m 
0.30 m 

€1 = eg = 0.0001 m 
62 = 67 = 0.0009 m 
^8 = Cii = 0.0008 m 
69 = 610 — 0.03° 


As mentioned in Section 5.6, the coordinates 62 and C2 are chosen such that 62 < h 
and C2 > Cl- 

For the laser beam to be horizontal and along Xi direction, the angles 6 and ^ are 
taken to be 90° and 0° respectively. For these values of 6 and (j)^ the plane of reflection 
EFQ is normal to the Fb^^o-plane as discussed in Section 5.5. This is a degenerate 
case and should be taken care of separately. The mirror plane is obtained for these 
nominal values and they cannot be varied in the present formulation. 

The standard deviation cr of error at a point on the work surface is given by [110] 


sjm + m + DN 

say. 


where 


Hence, 


VD^, 


DSUM = D[x] + D[y] + D[z] 


DSUM = cr2 

i.e., DSUM is the square of the standard deviation of error at a point. 

The variance D[x}, D[y] and D[z] are obtained at several points Q on the resin 
surface, i.e., the XiFrplane. The x and y coordinates of Q, which lie in the range 

0 < 9® < -30 and 0 < % < .30, are incremented by .03 each. Table 5.1 hsts the 

variances D[x], D[y] and D[z], and their sum DSUM for these values of q^ ajod gy. 



lYI 

qy 

m 

III 

0.00 

0.00 

0.03 

0.00 

0.06 

0.00 

0.09 

0.00 

0.12 

0.00 

0.15 

0.00 

0.18 

0.00 

0.21 

0.00 

0.24 

0.00 

0.27 

0.00 

0.30 

0.03 

0.00 

0.03 

0.03 

0.03 

0.06 

0.03 

0.09 

0.03 

0.12 

0.03 

0.15 

0.03 

0.18 

0.03 

0.21 

0.03 

0.24 

0.03 

0.27 

0.03 

0.30 

0.06 

0.00 

0.06 

0.03 

0.06 

0.06 

0.06 

0.09 

0.06 

0.12 

0.06 

0.15 

0.06 

0.18 

0.06 

0.21 

0.06 

0.24 

0.06 

0.27 

0.06 

0.30 

0.09 

0.00 

0.09 

0.03 

0.09 

0.06 

0.09 

0.09 

0.09 

0.12 

0.09 

0.15 

0.09 

0.18 

0.09 

0.21 

0.09 

0.24 

0.09 

0.27 

0.09 

0.30 

0.12 

0.00 

0.12 

0.03 

0.12 

0.06 

0.12 

0.09 

0.12 

0.12 

0.12 

0.15 

0.12 

0.18 

0.12 

0.21 

0.12 

0.24 

0.12 

0.27 

0.12 

0.30 


Variances 


and T heir I 

dFI 


Variances ana JL ueir aum 


0.15 

0.15 

0.15 

0.15 

0.15 


0.00 

0.03 

0.06 

0.09 

0.12 


2.1056 
2.1071 
2.1131 
2.1238 
2.1392 
2.1592 
2.1837 
2.2127 
2.2460 
2.2834 
2.3249 
1.9988 
1.9980 
2.0019 
2.0104 
2.0236 
2.0415 
2.0640 
2.0909 
2.1223 
2.1578 
2.1973 
1.9170 
1.9138 
1.9162 
1.9213 
1.9322 
1.9478 
1.9680 
1.9929 
2.0221 
2.0656 
2.0931 
1.8600 
1.8540 
1.8527 
1.8662 
1.8646 
1.8776 
1.8954 
1.9180 
1.9460 
1.9763 
2.0117 
1.8278 
1.8188 
1.8145 
1.8150 
1.8204 
1.8307 
1.8459 
1.8659 
1.8905 
1.9195 
1.9527 
1.8206 
1.8082 
1.8006 
1.7977 
1.7998 


10 “" 

1.5646 

1.4983 

1.4480 

1.4137 

1.3951 

1.3925 

1.4057 

1.4349 

1.4804 

1.5424 

1.6212 

1.5354 

1.4699 

1.4203 

1.3864 

1.3681 

1.3655 

1.3785 

1.4074 

1.4523 

1.5134 

1.5912 

1.5087 

1.4441 

1.3951 

1.3616 

1.3436 

1.3410 

1.3639 

1.3823 

1.4267 

1.4871 

1.6639 

1.4848 

1.4209 

1.3724 

1.3393 

1.3215 

1.3189 

1.3316 

1.3597 

1.4035 

1.4632 

1.5391 

1.4634 

1.4001 

1.3622 

1.3194 

1.3017 
1.2991 
1.3117 
1.3395 
1.3828 
1.4418 
1.5168 
1.4445 
1.3818 
1.3343 

1.3018 
1.2843 


0.71111 


0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 

0.71111 


DSUM 

4.3814 
4.3165 
4.2723 
4.2486 
4.2454 
4.2627 
4.3005 
4.3587 
4.4375 
4.5369 
4.6572 
4.2452 
4.1790 
4.1333 
4.1079 
4.1029 
4.1181 
4.1536 
4.2094 
4.2856 
4.3823 
4.4996 
4.1368 
4.0690 
4.0214 
3.9940 
3.9869 
3.9999 
4.0330 
4.0863 
4.1599 
4.2637 
4.3680 
4.0559 
3.9860 
3.9363 
3.9066 
3.8971 
3.9076 
3.9381 
3.9888 
4.0596 
4.1506 
4.2619 
4.0023 
3.9300 
3.8778 
3.8455 
3.8332 
3.8409 
3.8687 
3.9165 
3.9844 
4.0724 
4.1806 
3.9762 
3.9011 
3.8459 
3.8106 
3.7952 


0.15 
0.15 
0.15 
0.15 
0.15 
0.15 
0.18 
0.18 
0.18 
0.18 
0.18 
0.18 
0.18 
0.18 
0.18 
0.18 
0.21 
0.21 
0.21 
0.21 
0.21 
0.21 
0.21 
0.21 
0.21 
0.21 
0.21 
0.24 
0.24 
0.24 
0.24 
0.24 
0.24 
0.24 
0.24 
0.24 
0.24 
0.24 
0.27 
0.27 
0.27 
0.27 
0.27 
0.27 
0.27 
0.27 
0.27 
0.27 
0.27 
0.30 
0.30 
0.30 
0.30 
0.30 
0.30 
0.30 
0.30 
0.30 
0.30 
0.30 


0.15 
0.18 
0.21 
0.24 
0.27 
0.30 
0.00 
0.03 
0.06 
0.09 
0.15 
0.18 
0.21 
0.24 
0.27 
0.30 
0.00 
0.03 
0.06 
0.09 
0.12 
0.15 
0.18 
0.21 
0.24 
0.27 
0.30 
0.00 
0.03 
0.06 
0.09 
0.12 
0.15 
0.18 
0.21 
0.24 
0.27 
0.30 
0.00 
0.03 
0.06 
0.09 
0.12 
0.15 
0.18 
0.21 
0.24 
0.27 
0.30 
0.00 
0.03 
0.06 
0.09 
0.12 
0.15 
0.18 
0.21 
0.24 
0.27 
0.30 


O-'^m 
1.8071 
1.8194 
1.8366 
1.8585 
1.8850 
1.91.58 
1.8387 
1.8225 
1.8110 
1.8045 
1.8068 
1.8158 
1.8299 
1.8489 
1.8726 
1.9008 
1.8827 
1.8623 
1.8466 
1.8358 
1.8303 
1.8302 
1.8355 
1.8461 
1.8618 
1.8824 
1.9077 
1.9533 
1.9282 
1.9077 
1.8922 
1.8821 
1.8776 
1.8787 
1.8853 
1.8973 
1.9145 
1.9366 
2.0512 
2.0209 
1.9951 
1.9743 
1.9590 
1.9495 
1.9458 
1.9480 
1.9558 
1.9691 
1.9876 
2.1775 
2.1413 
2.1096 
2.0829 
2.0618 

2.0466 
2.0376 
2.0347 
2.0378 

2.0467 
2.0610 


‘mr 
1.2817 
1.2940 
1.3215 
1.3643 
1.4227 
1.4970 
1.4280 
1.3658 
1.3187 
1.2865 
1.2664 
1.2786 
1.3058 
1.3481 
1.4059 
1.4795 
1.4138 
1.3521 
1.3053 
1.2732 
1.2559 
1.2532 
1.2653 
1.2922 
1.3341 
1.3914 
1.4643 
1.4018 
1.3405 
1.2940 
1.2621 
1.2448 
1.2421 
1.2540 
1.2806 
1.3222 
1.3789 
1.4512 
1.3920 
1.3310 
1.2847 
1.2529 
1.2357 
1.2329 
1.2447 
1.2710 
1.3123 
1.3685 
1.4402 
1.3841 
1.3234 
1.2773 
1.2457 
1.2284 
1.2256 
1.2372 
1.2633 
1.3042 
1.3601 
1.4313 


DM 

10 - 

0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 
0.71111 


DSUM 


3.7999 
3.8245 
3.8692 
3.9340 
4.0188 
4.1239 
3.9778 
3.8995 
3.8408 
3.8021 
3.7843 
3.8055 
3.8468 
3.9082 
3.9897 
4.0914 
4.0076 
3.9255 
3.8630 
3.8202 
3.7973 
3.7945 
3.8119 
3.8494 
3.9070 
3.9849 
4.0831 
4.0662 
3.9797 
3.9127 
3.8654 
3.8381 
3.8308 
3.8438 
3.8770 
3.9306 
4.0046 
4.0989 
4.1542 
4.0629 
3.9908 
3.9384 
3.9058 
3.8935 
3.9016 
3.9301 
3.9792 
4.0488 
4.1389 
4.2727 
4.1758 
4.0980 
4.0397 
4.0013 
3.9833 
3.9859 
4.0091 
4.0531 
4.1179 
4.2034 


Table 5.1: Variances and Their Sum at Several Poi 
by varying Qx and % 


the. Resin Surface obtained 
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Prom the table it is observed that the variances D[x], D[y] and D[z] are of the 
same order. The variance D[z] doesn’t vary over the resin surface. This is because 
the influence coefficient is non-zero for i=ll only, and is constant over 

the resin surface. This can be verified from Eq. (5.93). Table 5.2 fists the influence 
coefficients at {qx = 0.0,% = 0.30). The variances D[x] and D[y] vary over the resin 
surface causing DSUM to vary over it. The maximum value of the sum of variances 
DSUM occurs at (% = 0.0, % = 0.30). 


The Influence CoelBEicients 

i 


7Jil\ 


1 

1.687e-09 

-4.470e-10 

3.338e-16 

2 

1.723e-01 

l.OOOe+00 

6.085e-17 

3 

-1.132e+00 

3.000e-01 

-1.876e-16 

4 

-9.852e-01 

2.611e-01 

-1.036e-16 

5 

1.500e-01 

8.705e-01 

2.936e-18 

6 

l.OOOe-00 

4.470e-10 

1.094e-16 

7 

-1.723e-01 

-7.703e-ll 

-4.865e-17 

8 

5.744e-0l 

2.568e-10 

O.OOOe-fOO 

9 

2.947e-01 

1.026e-f00 

4.163e-17 

10 

-1.362e-f00 

2.611e-01 

5.551e-17 

11 

5.574e-01 

-3.000e-01 

l.OOOe-bOO 


Table 5.2: The Influence Coefficients at % = 0.0 m and % = 0.3 m in SL 

The three-sigma bands of error in tracing several curves on the resin surface is 
obtained. The three-sigma band is the band of 3cr and — 3cr values of error at a 
point. Figure 5.10 shows the three-sigma band for an inclined line from (.01, .28) 
to (.28, .1) in the XiTl-plane. The three-sigma band decreases and then increases 
with the parameter s along the fine. Figure 5.11 shows the three-sigma band for an 
ellipse with center at (.15, .15), semimajor axis .14 and semiminor axis .10 m the 
XiTi-plane. The inclination angle of ellipse from Xj direction is 0°. The three-sigma 
band varies slightly with the parameter psi along the ellipse. From Table 5.1 it is 
observed that the three-sigma band varies appreciably over the resin surface for this 

set of input values. 
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Chapter 6 


OPTIMAL ALLOCATION OF 
TOLERANCES AND 
CLEARANCES 

The stochastic model of an RP process is discussed in Section 2.2.5. The coordi- 
nates of a point on the work surface traced by the laser beam or the tip of the extruder 
head is expressed as a function of random variables involved in the process. If there 
are n random variables Vi, V 2 , • • ■ , Ki involved in the process then the coordinates 
of the point are given by x{Vi, V 2 , 14), y{Vi, V 2 , 14) and z{Vi, V 2 , 

14 ). The section gives the methodology to associate any given range of a dependent 
variable x, y or z with its corresponding probability. For normally distributed ran- 
dom variables, the dependent variable may as well be taken as normal for more than 
five random variables. For a given probability, the range of a dependent variable is 
evaluated from the variance of the dependent variable. Equations (2.22) and (2.24) 
give expressions for the variances of the dependent variables x, y and z in terms of 
the influence coefficients (^) j (^) {dvi) • 

Chapters 3, 4 and 5 derive expressions for the influence coefficients in the Lami- 
nated Object Manufacturing (LOM), the Fused Deposition Modeling (FDM) and the 
Stereolithography (SL) process. The variances D[x], D[y] and D[z] are tabulated at 
the mean values of the random variables at several points on the work surface. The 
three-sigma bands of mechanical error in tracing several curves on the work surface 
for given tolerances and clearances are plotted. This forms the analysis part of the 
problem. The synthesis part is the inverse of the above problem. In synthesis, the 
designer has to decide about the levels of tolerances on ki s for certain allowable tol 
erance limits on x, y and z. One very obvious solution is to keep very strict tolerances 
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on Ws so that the output errors on x, y and z will be less than the specified lim- 
its. However, such a solution is impractical since it is very expensive to manufacture 
any component or assembly with allowances close to zero. As a matter of fact, it is 
desirable to give as much of tolerance as possible to keep the manufacturing costs 
low. 

Consider an RP process with n random variables Vi involved in its stochastic 
model. The random variables K have variances D[Vi]. The problem of optimal allo- 
cation of tolerances and clearances can be stated as follows: 

n 2 

minimize C — (6.1) 

i=i pi 


where 

II 

t) 


subject to 


n 

J2aipi = D[x] 

n 

(6.2) 


Ei'iPi = -Db) 

i=l 

n 

(6.3) 


= D[z] 

(6.4) 


where 


ai 




Ci 



The minimization problem is stated at that point P(x,y,z) on the work surface 
where the errors in x, y and z have been observed to be critical while conducting 
the error analysis. In other words, some study in analyzing mechanical errors should 
be made prior to undertaking the synthesis. On the basis of this study values of 
variances D[x], Djy] and D[z] should be specified at a particular point P(x,y,z) on the 
work surface. Moreover, in the statement of the problem above, it has been assumed 
that the production cost C is inversely proportional to the variances D[T^]. 


6.1 Optimization Using Lagrange Multiplier Tech- 
nique 


The coefficients in the above minimization problem are evaluated at the mean 

value of Vi. The expressions for the mean values of the random variables m[Vi] 
discussed earlier do not depend upon Ci or and hence on /?». So pt and * 
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independent. Hence the constraints are linear functions of design variables for the 
optimization of pi’s. We can take derivatives easily. So the optimization problem is 
solved using Lagrange multiplier technique. 

Using Lagrange multiplier technique, the modified cost function can be written as 


M 


C Xi — D[. 

+A3 [Y^Cipi - D[z] 


X 


where 


U=:l 

n 1 

i=l rt 


and Ai, A2 and A3 are Lagrange multipliers [100, 101]. 
The optimality conditions for minimum M are 

dM 


dpi 

dM 

dXi 


0 

dM 


i = 1,2, 


dM 


0 


dXz dXz 

Applying the optimality conditions to the above problem gives 


AjOj + Xzbi + A3Ci 

i= 1 , 2 ,..., n 


YaiPi-D[x] = 0 

fsssl 

YciPi-D[z] = 0 


i=l 


Let 


(Tl 

A2 

“ Ai' 


<7-2 

Kx 

D[x] 

and 

K 2 

“ D\yY 






Ax’ 

£M. 

D[z] 

then 0-1 and <73 can be obtained by solving the following two equations: 

f{(ri,(T2) = Ki 

and 

g{ai, 0-2) = JTj 


(6.5) 


( 6 . 6 ) 
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where 


Or, 


Similarly 


f {(^1,0-2) 




O/i 


_ 

/ Mpi 


E 


hi 

T' + X2bi + XzCi 




9i<^i,cr2) 


El 

i 


\j Ci + Cibi + C 2 Ci 


1 ... 

/ Oi + cri 6 i + C 2 Ci 


1 

di + cri 6 i + (T 2 Ci 

El, 

1 


( 6 . 7 ) 


( 6 . 8 ) 


For Cl and crj to be real, it is necessary that 


Oi + cibi + C2Ci >0, i = 1,2,. . . ,n 


( 6 . 9 ) 


The region given by the expression ( 6 . 9 ) is plotted in Fig. 6.1 for Oj < 0 and aj > 0. 
Feasible region of crj and C2 is where the expression ( 6 . 9 ) is satisfied for all i. 

The feasible region is bounded by two convex hulls (Fig. 6.2). One convex hull is 
formed by the straight lines Oi + (ri6i + cr2Ci = 0 , for those i for which Uj > 0 . The 
other convex hull is formed by a* < 0. 

The function f{ci , 0^2) and g{ci , cr2) can be plotted in the feasible region of (ui , C2 ) . 
This will give allowable range of Ki and K^. Then for given values of Ki and K2 in 
the allowable range, Ci and cr2 can be found by solving Bqs. ( 6 . 7 ) and (6.8). 

There is a feasible region of ci and <72. Therefore arbitrary choice of Ki and K2 is 
not allowed. A real solution of (Ai, A2, A3) cannot be obtained for Ki and K2 chosen 
in the infeasible region. But if a real solution is obtained then the expression ( 6 . 9 ) is 
verified. The solution satisfies the first of equations (6.6). This gives 

Oj + Cibi + C2Ci = ^ (6.10) 

^IPi 

So the expression ( 6 . 9 ) is verified if Ai obtained in the solution is positive. This 
verification avoids plotting f io'i,C2) and g(ci,C2) to find the allowable range of Ai 
and K2. 

The first of equations (6.6) can be rewritten as 

^ i = U...,n (6.11) 

-y/AiOi 4- X2bt + AaCi 
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Figure 6 . 1 : The region ai + cribi + cr^Ci > 0 

Substituting the above expression for pi in the last three equations of (6.6) gives 
three non-linear equations in three unknowns Aj, A2 and A3 as follows: 

^i(Ai,A 2, A3) = ]^ai(Aiai -f A26i + AaCj) - D[x] = 0 

i^sl 

h,(Xi.h,M} = i;6i(A,a, + A 26 , + A3cr‘''^-i)|!/l = 0 (6.12) 

= ^Ci(Aiai-l-A26i + A3Ci) - D[z] — 0 

i=l 

The above set of non-linear equations can be solved by Newton-Raphson method. 
The elements of the Jacobian [J] can be derived as follows. 

?h = y'a,(-f)(A,a,H-A2(., + A3Ci)-=/^ 

^Ai S ^ 

Sh = Vai(-^)(Aiai + A 26 i + A3Ci)-^/^ 

tti 2 

= y;a,(-|)(A,a. + A2M + Aj,Ci)-='* 

. dA$ ■ i=s:l ; 
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Similarly the expressions for |^, i = 1,2,3 and |^, i - 1,2,3 can be obtained. 
Then the Jacobian [JJ can be written as follows: 


i=i ^ 


h%CLi 

Cj Qj% 


Ojjhx 

b] 

Cib% 


CLiCi 

biCi 

c? 


(6.13) 


So analytical expression for Jacobian [J] is available and it is not needed to find 
the numerical derivatives of h*, i = 1,2,3. On solving Eq. (6.12) by Newton Raphson 
method the values of Xi , A 2 and A 3 can be obtained. Then p, s can be obtained from 

expression ( 6 . 11 ). 

6.2 Optimization Using Genetic Algorithms 

Genetic algorithms (G As) are computerized search and optimization algorithms based 
on the mechanics of natural genetics and natural selection. GAs minuc the suxv,^- 
ol-the-flttest principle of nature to make a search process [25). The variabte m 
are coded in some string structures. Mostly binary coded string are used. The length 
of the string is usually determined according to the desired so ^ 

In general, a fitness function Ax) ^ derived from the objective function 

X' +irtn Fnr maximization problem, the fitness function 

used in successive genetic operation. Form ^ n _ fu\ Fnr 

can be conddered to be the same as the object!™ fimetion or H-) - /W- 
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minimization problem the following function is often used. 




1 

1 + f{x) 


The fitness function value of a string is known as the string’s fitness. 

The operation of GAs begins with a population of random strings representing 
design or decision variables. Thereafter, each string is evaluated to find the fitness 
value. The population is then operated by three main operators — reproduction, 
crossover and mutation — to create a new population of points. The new population is 
further tested for termination. If the termination criterion is not met, the population 
is iteratively operated by the above three operators and evaluated. This procedure is 
continued until the termination criterion is met. One cycle of these operations and 
the subsequent evaluation procedure is known as a generation in GA’s terminology. 

As GAs use a coding of variables, they work with a discrete search space. Even 
though the underlying objective function is a continuous function, GAs convert the 
search space into a discrete set of points. GAs have also been developed to work 
directly with continuous variables (instead of discrete variables). In those GAs binary 
strings are not used. Instead, the variables are directly used. These are called real- 
coded GAs [25, 26|. 

A code developed by Deb et al [26] for real-coded GA is used in the present work. 


6.3 Optimum Tolerances in LOM 

Chapter 3 finds the position vector q(u, w) = [ x y z 1 ] of a point on the contour 
of a slice traced by the laser beam in LOM. There are nineteen random variables 
involved in the stochastic model of the LOM process. The influence coefficients Oi, bi 
and Ci and the partial derivatives required for them are derived in Section 3.2, pp. 50- 
66. The variance of the dependent v«iriables x, y and z and their sum at several 
points on the work surface is tabulated in Table 3.1, Section 3.3. Sections 3.3 also 
gives the plots of three-sigma bands of error in tracing several lines on the bilinear 
surface q(u,w). Erom the error analysis it is found that the maximum value of the 

sum of variances DSUM occurs at (u=1.0, w=1.0). 

Numerical results for the optimal allocation of tolerances are obtained with the 
mean values of variables and other input values mentioned in Section 3.3. Optimal 
allocation is based on the point where DSUM is maximum, that is, at (u=1.0, w=1.0) 
in uw parametric space on the bilinear surface ABDC. Assuming D[x] — 1.6076 x 
10“^’'m^, D[y| = 1.1111 X and D[z] = 2.4994 x the optimal values 

of tolerances obtained are listed in Table 6.1. The optimal allocation is done using 
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Lagrange multiplier techniques. Since a real solution for the Lagrange multipliers 
(Ai,A 2,A3) is obtained with positive Ai, therefore, the ratio of variances, Ki and 
Ki are in the allowable range and the solution for (Ai,A2,A3) exists. The solution 
obtained for the Lagrange multipliers are as follows: 

Ai = 8.4607x 10^^ A2 = 8.0997x 10^^ A3 = 1.6493 x 10^^ 


Optimal Allocation of Tolerances 


Tolerance 



Units 

Absolute tolerance on coordinate ai 

€1 

= 2.152e- 

-04 

m 

Absolute tolerance on coordinate bi 

€2 

= 5.382e- 

-06 

m 

Absolute tolerance on coordinate Cj 

^3 

= 7.111e- 

-07 

m 

Absolute tolerance on angle ai 

€4 

= 2.537e- 

-05 

rad 

Absolute tolerance on angle /3i 

€5 

= 2.133e- 

-03 

rad 

Tolerance per unit length on radius ri 


= 2.370e- 

-05 

m/m 

Tolerance per unit length on length k 

e? 

= 4.337e- 

-02 

m/m 

Absolute tolerance on coordinate 02 

^8 

= 1.271e- 

-07 

m 

Absolute tolerance on coordinate 62 

eg 

= 5.382e- 

-06 

m 

Absolute tolerance on coordinate C2 

eio 

= 7.111e- 

-07 

m 

Absolute tolerance on angle 02 

eii 

= 7.500e- 

-06 

rad 

Absolute tolerance on angle 02 

ei2 

= 5.551e- 

-06 

rad 

Tolerance per unit length on radius r2 

ei3 

= 4.963e- 

-07 

m/m 

Tolerance per unit length on length l 2 

ei4 

= 7.712e- 

-04 

m/m 

Absolute tolerance on coordinate 03 

ei6 

= 9.892e- 

-09 

m 

Absolute tolerance on coordinate 63 

eie 

= l.OOOe- 

-08 

m 

Absolute tolerance on coordinate C3 

ei7 

= 2.168e- 

-02 

m 

Absolute tolerance on angle 0 

€18 

= 1.430e- 

-08 

rad 

Absolute tolerance on angle ^ 

ei9 

= 2.168e- 

-02 

rad 


Table 6.1: Optimal Allocation of Tolerances in LOM 


Prom Table 6.1 it is observed that much stricter tolerances are aUocated on the 
variables ^ 5 , Vi^ and compared to other variables. This is because the corre- 
spondiEg influeiice coefficients are much larger. Stricter tolerances are also allocated 

on the variables Vs, V'lo and 1 ^ 3 , whereas large tolerances are demanded on the 
variables V„ Ur and Vi,. Therefore, the error is most sensitive to V^s, Vis an is, 
and least sensitive to 14, Vi 7 and Vi 9 - 
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6.4 Optimum Tolerances and Clearances in FDM 


Chapter 4 finds t.ho c;oordinates of a point on the contour of a slice traced by the nozzle 
tip Q- There an; fifteen random variables Vi, V 2 , ^15 involved in the stochastic 

model of the FDM process. The expressions for the influence coefficients Oj, bi and Cj 
and the partial derivatives needed for them are given in Sections 4.7, 4.8 and 4.9, (pp. 
86-88). The variances of the dependent variable x, y and z and their sum at several 
points on the work surface is tabulated in Table 4.1, Section 4.10. Section 4.10 also 
gives the plots of three-sigma bands of error in tracing several curves by the nozzle tip. 
From the error analysis it is found that the maximum value of the sum of variances 
DSUM occurs at h = 0.25 m and k = 0-25 m. • 

The optimal allocation of tolerances and clearances is done at the mean values 
of variables and other input values mentioned in Section 4.10. Optimal allocation is 
based on the location of nozzle tip where DSUM is maximum, that is, at h = 0.25m 
and k = 0.25 m. Assuming D[x] = 6.5 x lO'^m^, D[y] = 3.0 x lO-^m^ and D[z] = 19.5 
X 10"®m2, the optimal values of variances of variables Vi are obtained using Lagrange 
multiplier technique. Sinc.e a real solution for the Lagrange multipliers (Ai,A 2 ,A 3 ) 
is obtained with positive Ai, therefore, the ratio of variances, Ki and K 2 are in the 
allowable range and the solution for (Ai, A 2 , A 3 ) exists. The solution obtained for the 
Lagrange multipliers are as follows: 

Ai = 1.0396 xlO^'^ A2 = 1.2301 xl0^« A 3 = 9.6734x 10^^ 


The optimal values of variances a of variables Ui obtained from Lagrange mul- 
tiplier technique are listed in Table 6.2. The derivatives and , - , and 

hence » are identical for (i - 2,7). They are also identical for the indrees 1 in eaA 
of tho sot. (i = 3,6), (i = 4,8) and (i = 5,9,10,11). Therefore, from Eq. (&6) the 
optimum variances p. are identical for the indices i in each of these sets as obtained 

by Lagrange multiplier technique in Table 6.2. , . • i 

The optimum values of variances are also obtained using real coded genetic algo- 
rithm optimization method. The input values taken for GA runs are as o ows. 


Population size: 
Number of variables: 
Bounds on variables: 
Crossover probability: 
Mutation probability: 
Random seed number: 


150 

15 

rigid 

0.90 
0.06 
0.15 


Random seed numoer. nf runt? The 

Optimization is done by varying number of generation and ^ ^ 

op Jnm values of variances p. on variables Vi, obtained for 600 generations and 
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single run, are listed in lable 6.2 alongwith those obtained from Lagrange multiplier 
technique. 


Optimum Variances from Lagrange 

Optimum Variances from Real GA 

i 

Variance pi 

i 

— '' " - — 1 — 

Variance pi 

1 

9.018e-10 

1 

8.533e-10 

2 

3.2150-08 

2 

2.953e-08 

3 

3.101e-08 

3 

3.078e-08 

4 

9.()18e-10 

4 

7.374e-10 

5 

2.730C-08 

5 

2.873e-08 

6 

3.1010-08 

6 

3.123e-08 

7 

3.215e-08 

7 

2;912e-08 

8 

9.018e-10 

8 

7.405e-10 

9 

2.730e~08 

9 

2.920e-08 

10 

2.730(1-08 

10 

2.809e-08 

11 

2.730e-08 

11 

2.799e-08 

12 

6.531(1-09 

12 

1.445e-08 

13 

8.058(1-08 

13 

7.328e-08 

14 

6.005C-09 

14 

1.545e-08 

15 

6.940e-08 

15 

7.611e-08 


Table 6.2: Optimum Values of Variances pi in FDM from Lagrange Multiplier Tech- 
nique and from Real Coded GA. 


The optimum values obtained by real GA compare quite closely with those ob- 
tained by Lagrange method. The optimum values of pu and pu obtained by real 
GA differ from those obtained by Lagrange method. This needs to be investigated 
further. It has been observed that when the optimum values differ by several orders, 
then some of the optimum variances obtained by GA differ from those obtained by 
Lagrange methofl 

The tolerances and clearances are computed from the optimal values of variances 
Pi of variables Vi obtained from Lagrange method. Both pe and pr give the clearance 
Ci 2 in pair 12. The smaller of the two values is chosen for the clearance ci 2 because 
narrower clearance will ensure that the variances on the dependent variables doesn t 
exceed the specified limit. Similarly smaller of the two values of clearances obtained 
from the v?irianceH ps and ps is chosen for the clearances C 34 on pair 34. The optimal 

values of tolerances and clearances are listed in Table 6.3. 

Converting the tolerance coefficients into absolute tolerances on the respective 
links and comparing the tolerances on links we find that stricter tolerances are 
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Optimal Allocation of Tolerances 


Tolerance 

Units 

Absolute tolerance in positioning member 2 along link 1 

ei = 9.009e-05 

m 

Tolerance per unit length on hnk 2 

62 = 7.685e-03 

m/m 

Absolute tolerance in positioning member 4 along link 3 

63 — 5.283e-04 

m 

Tolerance per unit length on hnk 4 

64 = 9.009e-03 

m/m 

Tolerance per unit length on hnk 5 

65 = 3.305e-03 

m/m 

Radial Clearance in pair 12 

C 12 = 3.5226 — 04 

m 

Radial Clearance in pair 34 

C34 = 6.006e-05 

m 

Absolute tolerance in Z 2 direction at the head attachment 

610 = 4.957e-04 

m 

Absolute tolerance in Z 2 direction at the nozzle tip attachment 

611 = 4.957e— 04 

m 

Absolute tolerance on angle cki 

612 = 2.424e-04 

rad 

Absolute tolerance on angle pi 

613 = 8.516e— 04 

rad 

Absolute tolerance on angle 0:2 

614 = 2.325e— 04 

rad 

Absolute tolerance on angle ^2 

615 = 7.903e — 04 

rad 


Table 6.3: Optimal Allocation of Tolerances in FDM 


manded on links 1 and 4. It may be noted that these two Hnks are the only hnks 
in lo-direction. Stricter clearance is demanded in pair 34. Among angles, stricter 
tolerances are demanded on angles ai and a 2 , i.e., the angles about Ao-direction. 
Therefore, the error is most sensitive to Vi, V 4 , Vs, Vg, Via and Vu. 


6.5 Optimum Tolerances in SL 


Chapter 5 finds the position vector q(u, w) — [ x y z 1 ] of the point Q on the 
contour of a slice in SL. The laser beam in SL draws the contour on the resin surface. 
There are eleven random variables Vi, Va, • . • , Vii involved in the stochastic model of 
the process. The expressions for the influence coefficients Oj, bi and Ci and the partial 
derivatives needed for them are given in Section 5.11, pp. 114-118. The variances of 
the dependent variable x, y and z and their sum is tabulated in Table 5.1, Section 5.12. 
Section 5.12 also gives the plots of three-sigma bands of error in tracing several curves 
on the resin surface. The maximum value of the sum of variances DSUM occurs at 


(?a: — — •3)’ 

The optimal allocation of tolerances is done at the mean values of variables and 
the other input values mentioned in Section 12. The optimal allocaticn is based on the 
point at which the sum of variances DSUM is maximum, that is, at {q^ = 0, % = .3) 
on the resin surface. Assuming D[x] = 8.0 x 10“W, D[y] = 12.0 x 10 w and D[z] 
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_ 5.0 X the optimal values of variances of the variables are obtained using 

Lagrange multiplier technique as discussed in Section 6 . 2 . Since a real solution for 
the Lagrange multipliers (Ai, A2, A3) is obtained with positive Ai, therefore, the ratio 
of variances, Ki and K 2 are in the allowable range and the solution for Ai, A 2 and 
A 3 exists. No real solution was found when K 2 5.0. The solution obtained for the 
Lagrange multipliers are as follows: 

Ai = 7.740x 10^® A 2 = 3.164x10^'^ A 3 = -2.034x 10 ^® 

The optimal values of variances pi are listed in Table 6 . 4 . 

The optimum values of variances are also obtained using real coded genetic al- 
gorithm optimization method as discussed in Section 6 . 2 . Since the error at the 
point Q doesn’t depend upon the variable Vi, therefore, only the other ten variable 
are considered for optimization using GA. The influence coeSicients with respect to 
variable Vi are smaller by many orders compared to those for other variables. Prom 
Eq. ( 6 . 6 ) it can be seen that the terms corresponding to Vi contribute negligibly to 
the constraints. The input values taken for GA runs are as follows: 


Population size: 

100 

Ntimber of variables: 

10 

Bounds on variables: 

rigid 

Crossover probability: 

0.85 

Mutation probability: 

0.05 

Random seed number: 

0.51 

Optimization is done by varying number of generation and number of runs. The 


optimum values of variances pi on variables Vi, obtained for 600 generations and 
single run, are listed in Table 6.4 alongwith those obtained from Lagrange multiplier 
technique. 

The optimum values of variances obtained by real GA compare very closely to 
those obtained by Lagrange multiplier technique. This confirms that the optimum 
values obtained are global optimum. 

The tolerances on coordinates and angles, obtained from the optimal values of 
variances pi, are listed in Table 6.5. The tolerance allocated on the variable Vi is 
many orders higher compared to those on other variables. This only means that the 
error at the point Q doesn’t depend upon the variable Vi. This can be easily seen 
that changing the x-coordinate of the laser beam doesn’t affect the coordinates of the 
point Q. The absolute tolerance on other variables are of the same order for this set 
of input values. Among them, stricter tolerance are demanded on the z-coordinate of 
the source E and the x-coordinate of mirror. Stricter tolerance is demanded on the 

angle /?. Therefore the error is most sensitive to V 3 , Ve and Pio- 
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Optimum Variances from Lagrange 

Optimum Variances from Real GA 

i 

Variance pi 

i 

- -'V ' — 1 

Variance pi 

1 

6.730e+00 


^ - 

2 

4.279e-08 

2 

4.463e-08 

3 

1.003e-08 

3 

1.147e-08 

4 

1.152e-08 

4 

CO 

o 

0 

CD 

1 

o 

oo 

5 

4.915e-08 

5 

4.807e-08 

6 

1.137e-08 

6 

1.300e-08 

7 

6.596e-08 

7 

7.647e-08 

8 

1.979e-08 

8 

2.280e-08 

9 

3.153e-08 

9 

3.253e-08 

10 

8.341e-09 

10 

9.546e-09 

11 

5.000e-08 

11 

2.239e-08 


Table 6.4: Optimum Values of Variances pi in SL from Lagrange Multiplier Technique 
and from Real Coded GA. 


Optimal Allocation of Tolerances 



Tolerance 

Units 

Absolute tolerance on coordinate ai 


= 7.782e+00 

m 

Absolute tolerance on coordinate bi 

^2 

= 6.205e- 

-04 

m 

Absolute tolerance on coordinate Cj 

€3 

= 3.004e- 

-04 

m 

Absolute tolerance on angle oi 

€4 

= 3.220e- 

-04 

rad 

Absolute tolerance on angle pi 

€5 

= 6.651e- 

-04 

rad 

Absolute tolerance on coordinate 

66 

= 3.198e- 

-04 

m 

Absolute tolerance on coordinate 62 

67 

= 7.705e- 

-04 

m 

Absolute tolerance on coordinate C 2 

^8 

= 4.220e- 

-04 

m 

Absolute tolerance on angle 02 

^9 

= 5.327e- 

-04 

rad 

Absolute tolerance on angle P 2 

610 

= 2.740e- 

-04 

rad 

Absolute tolerance on coordinate C 3 

^11 

= 6.708e- 

-04 

m 


Table 6.5: Optimal Allocation of Tolerances in SL 
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Th<i synt hesis procrediin; in the present work enables the designer to allocate 
tolerances and ch'araiKH's optimally. The synthesis procedure clearly shows that 
the allowances on diffen'iit- variables are different depending upon the influence of 
a particnilar variable. 



Chapter 7 

CONCLUSIONS 


7.1 Conclusions 


In the present work the mechanical error in rapid prototyping processes is studied 
using stochastic approach. Three RP processes, LOM, FDM and SL are considered. 
The coordinates of a point on the work surface traced by the laser beam or the tip 
of the extruder head is expressed as a function of random variables involved in the 
process. The methodology to associate any given range of a dependent variable x, y 
or z with its corresponding probability is given. The variances of the dependent vari- 
ables X, y and z are expressed in terms of the influence coefficients (^) , (|^) and 


f— 

In LOM, the veriances DM, D[yl and DN are of the same order. The r^riance DW 
wies significantly over the work surface compared to DM] and D[y|. The variance 
D[yl does not vary over the work surface and is equal to the variance of the random 
variable Vie- The variances and their sum are listed in a table to show their variation 
across the work surface. The thiecsigma bands of error in tracing few lines on the 
work surface is plotted. Phom the list of variances and the three-sigma bands it is 

, vxi xi-v. ..„„,.oUprrnr does not vary much over the work surface and a contour 

observed that the overall error aoebuu j 

, ■. the work surface without much difference m accuracy, ine 

can be drawn anywhere on tne woriv ^ nn 

1 • • T rvM iR done when the laser beam traces the contour of a shce on 

error analysis in LOM is aone 

the sheet held between idle rollers. . nM nn 

The variance DM and Dly] in FDM are of the same cider. The variance DM is an 
order higher. This is because the influence coefficients (^) have larger magmtudes 

J i. The three-sigma bands in tracing few example curves 

compared to and ^ 

by the nozzle tip IS plotted, -the variai 
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show their variation across the work surface. The overall error is found to be varying 
appreciably across the work surface. The error is the least at the front-left end of the 
work surface and maximum at the rear-right end. The mechanical error in FDM was 
studied assuming that the material solidifies as it comes out of the nozzle tip. 

The variances D[x], D[y] and D[z] in SL are of the same order. The variance D[z] 
does not vary over the work surface. The variances and their sum are listed in a table 
to show their variation across the work surface. The three-sigma bands of error in 
tracing a line and an ellipse on the resin surface is plotted. 

A methodology for optimal allocation of tolerances and clearances is presented. 
Since the constraints are linear functions of design variables, therefore, the optimiza- 
tion is done using Lagrange multiplier technique. The optimal allocation is done 
at that point on the work surface where the sum of variances DSUM is found to 
be maximum while conducting error analysis. Optimization is also done using real 
coded Genetic Algorithm {real coded GA) for FDM and SL. The optimum values of 
variances for LOM differ from each other by several orders. The real GA code does 
not give good results for this case. Therefore, optimization in the case of LOM is not 
done using GA. 

The synthesis procedure in the present work enables the designer to allocate tol- 
erances and clearances optimally. The synthesis procedmre clearly shows that the 
allowances on different variables are different depending upon the influence of that 
variable. In LOM, stricter tolerances are demanded in the x and y-coordinates of the 
source of the laser beam and the angle 9 of the laser beam. In FDM, stricter toler- 
ances are demanded on links 1 and 4, i.e., the hnks in Fo-dfrection and on angles ax 
and 0 ! 2 , i.e., the angles about Ao-direction. In SL, stricter tolerances are demanded 
on the z-coordinate of the source of the laser beam, the x-coordinate of the mirror and 
the angle of of the mirror. It has been found that the error in SL is not sensitive 
to the location of the source of the laser beam along the direction of the laser beam. 

7.2 Scope for Further Work 

Geometric modeling of the LOM process can be improved by considering the mecha- 
nisms which cause the movement of the optics and the platform. The tolerances and 
clearances in these mechanisms can be considered in the error studies. The tolerances 
in the belt-driven mechanism for XY-positioning table can be considered. 

There is a positive displacement belt drive for the motion of extruder head in 
FDM. The Z-stage platen moves on the power screws and there are two cylindrical 
rods for maintaining its alignment. The geometric model of FDM can be improved by 
considering these mechanisms. The manufacturing errors in them can be consider 
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in the error studies. 

The geometric model of SL has considered the reflection through the mirror. It 
has considered the error in location and orientation of the mirror. In this way, it is an 
improvement over the geometric model of LOM. The geometric model can be further 
improved by considering the mechanisms which cause the movement of the optics and 
the elevator platform. 

The study of mechanical error in FDM is carried out assuming the material solid- 
ifles as it comes out of the nozzle tip. The material flow characteristics and solidifi- 
cation properties can be considered in further work. 

In the present work the mechanical error has been studied in drawing the slice 
of a part. The cumulative error in FDM, as the part is built progressively, depends 
upon the solidification and flow characteristics of modehng material and the error in 
the platform. 

The cumulative error in LOM depends more on the error in the location and orien- 
tation of the platform and the heated roller. Geometric model is to be further refined 
taking them into consideration. The error at a particular layer is still stochastic as 
the thickness of the LOM sheet is sensed and the solid model is sliced on the fly. 

The cumulative error in SL depends upon the error in positioning the platform 
and the error in sweeping by the recoater blade. Since the resin level is sensed after 
every layer, therefore, the error at any layer can be found out by stochastic approach. 
Further work is needed to find the cumulative error in RP processes as the part is 
built progressively. 

In the present work synthesis has been carried out at a particular position on the 
work surface. The choice of the position is based on the trial analysis of the RP 
process. This approach seems to yield satisfactory results. However, a methodology 
can be developed for finding that critical point. For instance, the critical point can 
be chosen by repeating the synthesis procedure throughout the work surface. The 
critical point would be the point for which the maximum value of the sum of variances 
over the work surface is minimum. The optimal values of tolerances and clearances 
correspond to this point. The function /(u-i.crz) and g(cri,(T 2 ) can be plotted in 3D 
in the feasible region of (cri, cr 2 ) to find the allowable range of Ki and in Lagrange 
multiplier technique. 

Optimization using real coded GA for the cases of FDM and SL gives optimal 
values comparable to those obtained by Lagrange multiplier, technique. Results of 
real coded GA can be further refined. Optimization using real coded GA in the case 
of LOM needs to be further investigated. The optimum values of variances in LOM 
differ from each other by several orders. The code for real coded GA needs to be 

improved to take into account this case. 
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